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Abstract 

In this review we discuss the important progress made in recent years towards under- 
standing the experimental data on cosmic rays with energies <^ 10 19 eV. We begin with 
a brief survey of the available data, including a description of the energy spectrum, 
mass composition, and arrival directions. At this point we also give a short overview 
of experimental techniques. After that, we introduce the fundamentals of acceleration 
and propagation in order to discuss the conjectured nearby cosmic ray sources. We 
then turn to theoretical notions of physics beyond the Standard Model where we con- 
sider both exotic primaries and exotic physical laws. Particular attention is given to 
the role that TeV-scale gravity could play in addressing the origin of the highest energy 
cosmic rays. In the final part of the review we discuss the potential of future cosmic 
ray experiments for the discovery of tiny black holes that should be produced in the 
Earth's atmosphere if TeV-scale gravity is realized in Nature. 
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I. APPETIZER 



Cosmic ray (CR) and accelerator based particle physics share common roots, and in 
fact many of the key discoveries early in the history of particle physics came from the 
study of CRs. After a period of divergence between the two fields, both in methodology 
and in the key areas of interest, a confluence is now underway, driven in no small part by 
the mystery of the highest energy CRs. Much of the interest in these CRs is rooted in 
recent developments in both experiment and theory. The most recent data reported at 
the 27th International Cosmic Ray Conference (ICRC) suggest that some CRs arrive 
at the outer limits of the Earth's atmosphere with energies above 10 20 eV. This is 
in itself remarkable, as it seems difficult to explain how such high energies can be 
attained without invoking new physics. Moreover, the center-of-mass (cm.) energies 
achieved when these CRs impinge on a stationary nucleon of the air molecules are 
higher than any that can be reached in present-day terrestrial experiments, or indeed 
any experiment in the foreseeable future. Adding to the puzzle, the HiRes and AGASA 
experiments continue to revise their estimates of the fluxes and energies, making the 
subject more confusing. 

Although current theoretical and experimental uncertainties make it impossible to 
determine the origin and nature of the highest energy events with any certainty, the 
mere presence of such events suggests the exciting possibility that CRs will, once again, 
provide us with glimpses of physics beyond our present theories. Speculations ranging 
from supersymmetry to magnetic monopoles, large extra dimensions, and captivating 
tiny black holes produced in CR collisions fill the Los Alamos preprint archives at an 
ever-increasing rate. There is clearly great hope that we are at the verge of learning 
something revolutionary. 

In this spirit then, it seems opportune to review our present understanding of the 
highest energy CRs: what has been measured and with what confidence, what do we 
know about the composition and possible acceleration mechanisms for these particles, 
and how much do our measurements depend on uncertainties in Standard Model cal- 
culations? Throughout the first part of this review we focus tightly on the interplay 
between experiment and phenomenology assuming no physics beyond the Standard 
Model. This provides something of a sanity check on where we are now and where we 
think we are going. Afterwards, we relax our grip on sanity and summarize some of 
the new, and sometimes fantastic, ideas that may be probed by the next generation of 
experiments. To paraphrase Neils Bohr, the hope is that none of these ideas is crazy 
enough, and that in a few years we can look forward to genuine surprises. 
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II. OBSERVATION OF THE HIGHEST ENERGY COSMIC RAYS 



A. The energy spectrum 



In 1912 Victor Hess carried out a series of pioneering balloon nights during which he 
measured the levels of ionizing radiation as high as 5 km above the Earth's surface [TJ. 
His discovery of increased radiation at high altitude revealed that we are bombarded by 
ionizing particles from above. These CR particles are now known to consist primarily 
of protons, helium, carbon, nitrogen and other heavy ions up to iron. 

Below 10 14 eV the flux of particles is sufficiently large that individual nuclei can be 
studied by detectors carried aloft in balloons or satellites. From such direct experiments 
we know the relative abundances and the energy spectra of a variety of atomic nuclei, 
protons, electrons and positrons as well as the intensity, energy and spatial distribution 
of X-rays and 7-rays. Measurements of energy and isotropy showed conclusively that 
one obvious source, the Sun, is not the main source. Only below 100 MeV kinetic 
energy or so, where the solar wind shields protons coming from outside the solar sys- 
tem, does the Sun dominate the observed proton flux. Spacecraft missions far out into 
the solar system, well away from the confusing effects of the Earth's atmosphere and 
magnetosphere, confirm that the abundances around 1 GeV are strikingly similar to 
those found in the ordinary material of the solar system. Exceptions are the overabun- 
dance of elements like lithium, beryllium, and boron, originating from the spallation 
of heavier nuclei in the interstellar medium. 

Above 10 14 eV, the flux becomes so low that only ground-based experiments with 
large apertures and long exposure times can hope to acquire a significant number of 
events. Such experiments exploit the atmosphere as a giant calorimeter. The incident 
cosmic radiation interacts with the atomic nuclei of air molecules and produces exten- 
sive air showers (EASs) which spread out over large areas. Already in 1938, Pierre 
Auger concluded from the size of EASs that the spectrum extends up to and perhaps 
beyond 10 15 eV @, HI- Nowadays substantial progress has been made in measuring the 
extraordinarily low flux (~ 1 event km~ 2 yr -1 ) above 10 19 eV. Continuously running 
experiments using both arrays of particle detectors on the ground and fluorescence de- 
tectors which track the cascade through the atmosphere, have detected events with pri- 
mary particle energies higher than 10 20 eV [(§ §, §, 0, §, |, 0, [IT], 0, [13], [IJ, |IJ, f§ |I7| , 
with no evidence that the highest energy recorded thus far is Nature's upper limit. 

In contrast to the irregular shape of the isotropic electromagnetic background spec- 
trum from, say, 10 8 — 10 20 Hz, the CR energy spectrum above 10 9 eV can be described 
by a series of power laws, with the flux falling about 3 orders of magnitude for each 
decade increase in energy (see Fig. |l]). In the decade centered at ~ 10 15 ' 5 eV (the 
knee) the spectrum steepens from E~ 2 ' 7 to E~ 3 ' . This feature, discovered around 40 
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FIG. 1: Compilation of measurements of the differential energy spectrum of CRs. The dotted 
line shows an E~ s power-law for comparison. Approximate integral fluxes (per steradian) 
are also shown fiUM. 



years ago [IS , is still not consistently explained. The spectrum steepens further to 
E~ 3:i above ~ IO 17 ' 7 eV (the dip) and then flattens to E~ 2 - 7 at ~ IO 18 ' 5 eV (the ankle). 
Within the statistical uncertainty of the data collected by AGASA which is large 
above 10 20 eV, the tail of the spectrum is consistent with a simple extrapolation at 
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that slope to the highest energies, possibly with a hint of a slight accumulation around 
10 195 eV. A very widely held interpretation of the ankle is that above 10 18 5 eV a 
new population of CRs with extragalactic origin begins to dominate the more steeply 
falling Galactic population [^TJ. This hypothesis is supported by the data recorded at 
AGASA, which shows that around 10 18 eV the angular distribution correlates with the 
Galactic center (anisotropy ~ 4%) and is consistent with a Galactic origin, whereas at 
higher energy the anisotropy disappears [T2^, [23[ . 

At the very high end of the spectrum the flux now appears more uncertain than 
was thought as recently as a year ago. At that time the event rate at 10 19 eV re- 
ported by different experiments was in agreement at the 10 - 15% level. In particular, 
preliminary data reported by the HiRes group [ 16j showed 7 events above 10 20 eV, in 
good accord with the number anticipated from the flux observed by the AGASA ex- 
periment pO| . Very recently, the AGASA Collaboration reported a total of 17 events 
above 10 20 eV |T7[ , quite consistent with their previous work PD| . However, one of 
the HiRes cameras, which has an exposure slightly greater than that of the AGASA 
experiment, recorded only 2 events above 10 20 eV compared to the 20 expected for a 
spectrum similar to that reported earlier by AGASA ||24j| . The HiRes group also re- 
ported data from their stereo system, which has 20% of the monocular exposure. They 
observed one event with an energy estimated to be close to 10 20 ' 5 eV. These unexpected 
discrepancies are not yet understood. The "disappearance" of the events reported as 
being above 10 20 eV in 1999 [16| is attributed to a better understanding of the atmo- 
sphere, which is now claimed to be clearer than had been previously supposed. Adding 
to the puzzle, the Haverah Park energy estimates have been re-assessed [f25|, resulting 
in a steeper reconstructed spectrum which shows differences of up to 30% compared to 
the one derived by fitting Akeno, AGASA and the "old" Haverah Park data PEfl. 1 The 
average energy of the 4 events observed by the Haverah Park experiment that were pre- 
viously above 10 20 eV is shifted 30% downwards to energies below 10 20 eV. The energy 
weighted cosmic ray flux corresponding to these 4 events is JE 3 = 8.3^4q x 10 24 eV 2 
gr -i m -2 g -i a £ an ener gy Q f 7_g2 x 10 19 eV (see Fig. |113| ). The highest energy Haverah 
Park event is now estimated to have an energy of 8.3 x 10 19 eV. In the "old" Haverah 
Park data sample, this event had a reported energy of 1.2 x 10 20 eV. In Fig. ^|we show 
a compilation of the updated data from Haverah Park |2"5| , AGASA [^(]], monocular 
HiRes (HiRes 1 and HiRes 2) |4j, HiRes-MIA |§ and Fly's Eye stereo 0. For 



comparison, the figure also shows 4 recent parameterizations of the energy spectrum 



It is noteworthy that AGASA energies have been estimated under the assumption that the primaries 
are protons above 10 18 eV. Though there is no evidence as to what mass species is dominant at 
the highest energies, the results from AGASA decrease by only about 20% when the hypothesized 
chemical composition is changed. 
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including the fit to the Akeno, AGASA and Haverah Park data sets by Nagano and 
Watson f26j, the one obtained by Szabelski, Wibig and Wolfendale f29j, the param- 
eterization given by the HiRes Collaboration [33] and a parameterization based on a 
recent analysis of Haverah Park data p5 |. The recalibrated Haverah Park spectrum 
is in very good agreement with that obtained by the fluorescence detectors of Fly's 
Eye and HiRes. Reconciling the observations of the various experiments depends on a 
detailed understanding of detector behavior and atmospheric properties as well as reli- 
able models to predict the evolution of cosmic air showers. These issues are discussed 
in the following sections. 



B. Properties of extensive air showers and mass composition 

We begin this section with a brief discussion of the principal observables characteriz- 
ing an EAS. Next, we discuss how the shower development differs for different primary 
species and how these differences are manifest in the shower observables. 

The incidence of a single high energy particle on the upper atmosphere gives rise 
to a roughly conical cascade of particles which reaches the Earth in the form of a 
giant "saucer" traveling at nearly the speed of light. As the cascade develops in the 
atmosphere, the number of particles in the shower increases until the energy of the 
secondary particles is degraded to the level where ionization losses dominate. At this 
point the density of particles starts to decline. The number of particles as a function 
of the amount of atmosphere penetrated by the cascade (in g cm -2 ) is known as the 
"longitudinal profile" shown schematically in Fig. [5]. The atmospheric depth at which 
the number of particles in the showers reaches its maximum, A max , is often regarded as 
the most useful observable of the shower, as it strongly depends on the primary energy 
and composition. For example, A max increases with primary energy since more cascade 
generations are required in the cooling of secondary products. The way the average 
depth of maximum (X max ) changes with energy depends on the primary composition 
and particle interactions according to 

(X max )=D e \n(^-) , (1) 



where D e is the so-called "elongation rate" and E is a characteristic energy that 



depends on the primary composition [[30]. Therefore, since (X max ) and D e can be 



determined directly from the longitudinal shower profile measured with a fluorescence 
detector, E and thus the composition, can be extracted after estimating E from the 
total fluorescence yield. Indeed, the parameter often measured is D w , the rate of 
change of (A max ) per decade of energy. We note that one can discern changes in the 
primary composition from breaks in the elongation rate, and such breaks are relatively 
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FIG. 2: A composite energy spectrum including recently reanalysed Haverah Park data 
assuming proton and iron primaries (the parameter A measures the attenuation length of the 
density of charged particles at 600 m from the shower core), stereo Fly's Eye data, monocular 
HiRes data from both eyes up to 60°, and hybrid HiRes-MIA data. The measurements are 
compared to spectrum parameterizations given by different authors as described in the text. 
Published in Ref. (27 ]. 



insensitive to certain systematic uncertainties. 

If the primary particle is a nucleon or a nucleus the shower begins with a hadronic 
interaction. The number of hadrons increases through subsequent generations of parti- 
cle interactions. However, in each generation about 30% of the energy is transferred to 
an electromagnetic cascade by the prompt decay of neutral pions. Ultimately, the elec- 
tromagnetic cascade dissipates around 90% of the primary particle's energy, and hence 
the total number of electromagnetic particles is very nearly proportional to the shower 
energy. The remaining energy is carried by muons and neutrinos from 7r =l= decays. The 
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FIG. 3: Particles interacting near the top of the atmosphere initiate an electromagnetic and 
hadronic cascade. Its profile is shown on the right. The different detection methods are illus- 
trated. Mirrors collect the Cerenkov and nitrogen fluorescent light, arrays of detectors sample 
the shower reaching the ground, and underground detectors identify the muon component of 
the shower. 

number of muons does not increase linearly with energy, since at higher energy more 
generations are required to cool the pions to the point where they are likely to decay 
before interaction. Production of extra generations results in a larger fraction of the 
energy being lost to the electromagnetic cascade, and hence a smaller fraction of the 
original energy being delivered to the 7r ± . With this in mind, one can evaluate the 
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muon production in a heavy nucleus shower relative to that in a proton shower by 
viewing a nucleus with atomic number A and energy E as a collection of individual 
nucleons each with energy E/A. The muon production in a proton shower increases 
with energy as E ' 85 ||31||, and thus the total number of muons produced by a nucleus of 



mass A is N A oc A(E/A) - 85 , or, comparing to proton showers, Nw~ = A°' 15 NJ^. There- 
fore, an iron nucleus produces a shower with around 80% more muons than a proton 
shower of the same energy. Around 10 20 eV, the hadronic mean free path is roughly 
40 g/cm 2 . 

On the other hand, if the primary particle is a 7-ray the basic interactions that cool 
the cascade development are pair production, Bremsstrahlung, ionization losses, and 
Compton scattering. Above 10 19 eV the Landau- Pomeranchuk-Migdal (LPM) effect 
leads to significant suppression of the Bethe-Heitler cross section for pair production 
and Bremsstrahlung, and so the growth of with energy gradually disappears |]32[ . 
However, the muon suppression in electromagnetic showers caused by the LPM ef- 
fect become largely compensated by the interactions of 7-rays with the geomagnetic 
field [Q. Above the ankle (~ 10 18,5 eV), a 7-shower produces less than 20% as many 
muons as a proton shower of the same energy. The energy dependence of this effect has 
a rather complicated (non power-law) form |34 ]. Gamma-rays of energy ~ 3 x 10 20 eV 
have interaction lengths of approximately 60 g/cm 2 with a standard deviation of 80 
g/cm 2 if the LPM effect is taken into account, but these figures reduce to 46 g/cm 2 
and 45 g/cm 2 , respectively, if the LPM effect is ignored f32|| . 

The experimental observables for inclined showers (incident zenith angle 9 > 60°) 
are complimentary to those for vertical showers, since inclined showers are mostly 
comprised of muons which were produced far from the detection zone [ 35| . In contrast 
to vertical showers the shower front is extremely flat (radius of curvature > 100 km) and 
the particle time spread is very narrow (At < 50 ns). For 9 > 70° the electromagnetic 
component arising from the hadronic channel through 7r° decays is totally extinguished, 
as more than 3 equivalent vertical atmospheres are traversed. There remains only a 
residual electromagnetic component in the shower front which is produced by the muons 
themselves, mostly through muon decay. 

Extracting information on the nature of the primaries from D e and (A max ) , or from 
the variation of the muon content has proved to be exceedingly difficult for a number 
of reasons. The most fundamental drawback is that the first few cascade steps are 
subject to large inherent fluctuations and consequently this limits the event-by-event 
mass resolution of the experiments. In addition, the cm. energy of the first interactions 
is well beyond any reached in collider experiments. Thus, one needs to rely on hadronic 
interaction models that attempt to extrapolate, using different mixtures of theory and 
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phenomenology, our understanding of particle physics. 2 

The analysis of the elongation rate and the spread in X max at a given energy reported 
by the Fly's Eye Collaboration suggests a change from an iron dominated composition 
at 10 17 ' 5 eV to a proton dominated composition near 10 19 eV fll0 |. Such behavior of D e 
is in agreement with an earlier analysis from Haverah Park [J36|J . However, the variation 
of the density of muons with energy reported by the Akeno Collaboration favours a 
composition that remains mixed over the 10 18 — 10 19 eV decade [[57]. At this stage, it 
is worthwhile to recall that any inferred mass composition would be model-dependent 



(see Appendix A). Recently, Wibig and Wolfendale |38[ reanalyzed the Fly's Eye data 
considering not only proton and iron components (as in [JTOj] ) but a larger number of 
atomic mass hypotheses. Additionally, they adopted a different hadronic model that 
shifts the prediction of A max for primary protons of 10 18 eV from 730 g cm~ 2 |10 



to 751 g cm -2 . The difference, although apparently small, has a significant effect on 
the mass composition inferred from the data. The study indicates that at the highest 
energies (10 18,5 — 10 19 eV and somewhat above) there is a significant fraction of nuclei 
with charge greater than unity. This result is more in accord with the conclusions of 
the Akeno group than those of the Fly's Eye group. Deepening the problem, there 
is a difference between the conclusion reached from data on A max from the Fly's Eye 



experiment 10] and those reached from the HiRes prototype, operated with the MIA 



detector in the range 10 17 to 10 18 eV [[39 



Nowadays, the most powerful tool for extracting the primary mass spectrum relies 
on comparing the flux of vertical showers to that of inclined showers. For example, if all 
primaries were photons, one would expect to observe a great asymmetry in the number 
of vertical and inclined events. Ave et al. {4(| f|T] have used the observed Fly's Eye 
spectrum to predict the rate of inclined events. Comparing the predicted rate to the 
rate observed by Haverah Park for showers in the range 60° < 9 < 80°, they conclude 
that above 10 19 eV, less than 48% of the primary cosmic rays can be photons and above 
4 x 10 19 eV less than 50% can be photons. Both of these statements are made at the 
95% CL. This analysis is consistent with the search for showers which have significantly 
fewer muons than normal as reported by the AGASA Collaboration |42j| . In addition, 
according to the predictions of QGSJET for high energy interactions (see Appendix 
A), the difference in abundance of muons with respect to photons and electrons in the 
Haverah Park inclined shower measurements seems to favour a light composition above 
10 19 eV [I3|. A more robust statement awaits a better understanding of the sensitivity 



of this kind of analysis to different models of particle interactions. 

Finally, if the primary is a neutrino the first interaction occurs deep in the atmo- 
sphere, triggering showers in the volume of air immediately above the detector. The 



A brief summary of the current status of hadronic models is provided in Appendix A. 
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shower would thus present a very curved front (radius of curvature of a few km), with 
particles arriving over O(fis). If the primaries are u e and as expected from tc^ 
decays, one expects a different type of shower for each neutrino species: "mixed" (with 
full energy) for v e or "pure hadronic" (with reduced energy) for v^. In the charged 
current interaction of a v e , an ultrahigh energy electron is produced which initiates a 
large electromagnetic cascade parallel to the hadronic cascade. In contrast, the charged 
current interaction of a produces an ultrarelativistic muon that is not detectable 
by the experiments. In the presence of maximal v il jv T mixing, v T showers must also 
be considered. However, since the r mean flight distance is ~ 50E km/EeV, only r's 
with energy < 8 x 10 17 eV will decay. Thus, v T showers above this energy will be 
indistinguishable from showers. The characteristics of these anomalous showers can 
be easily identified by fluorescence eyes and surface arrays. 



C. The arrival direction distribution 

The investigation of anisotropy, when taken together with analyses of the spectral 
shape and particle species, can yield very important clues to reveal the CR origin. 
For example, if CRs mostly come from the direction of the Galactic plane and are 
largely protonic in composition, one expects to see a dipole anisotropy favoring the 
direction of the Galactic center. Moreover, since magnetic rigidity increases with en- 
ergy, the angular width of the Galactic plane as seen in CRs would shrink slowly with 
rising energy. Of course, a heavy composition dominated by iron nuclei would show 
much smaller anisotropy at any given energy because of the smaller Larmor radius. 
The deviation from isotropy in galactic latitude is generally expressed in terms of the 
Wdowczyk-Wolfendale plane enhancement function, 



= C 1 - + L437 f* ex P [- fe2 ] > ( 2 ) 



where b is the galactic latitude in radians, J b s and I exp are observed and expected 
(for isotropy) intensities at latitude b, and is an energy dependent galactic latitude 
enhancement factor. A galactic origin for most of the particles would be expected to 
result in a positive value of Je that increases steadily with energy. A negative shows 
depression around the plane and = indicates the arrival direction distribution 
is isotropic. Much of the data suggests that the magnitude of the Galactic plane 
enhancement increases systematically with energy until a little more than 10 19 eV, 
above which it disappears and, indeed, there is evidence for a deficit in the direction 
of the Galactic plane above this energy The most significant Galactic plane 

enhancement factor reported by Fly's Eye, J'e = 0.104 ± 0.036, is in the energy range 
(0.4 — 1.0) x 10 18 eV ||46|| . In a similar energy range, the AGASA Collaboration reported 
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a strong anisotropy with first harmonic amplitude of 



probability of 0.2% after taking the number of independent trials into account fl22 



corresponding to a chance 

A 



recent analysis of SUGAR data confirms the existence of an excess flux of CRs from a 
direction near the Galactic center |47 |. The signal is consistent with that from a point 
source H8|. 




FIG. 4: Arrival directions of cosmic rays detected by the AGASA and Akeno (A20) experi- 
ments in equatorial coordinates. Open circles and open squares represent cosmic rays with 
energies (4 - 10) x 10 19 eV, and > 10 20 eV, respectively. The Galactic and super-Galactic 
planes are shown. Large shaded circles indicate event clusters within 2.5°. The shaded re- 
gions indicate the celestial regions excluded by a zenith angle cut of < 45°. Published in 
Ref. POL 



In a much higher energy range, the analysis of Haverah Park data seems to indicate 
that events with primary energy > 4 x 10 19 eV reach the Earth preferentially from the 
direction of the super-Galactic plane, a swath in the sky along which radio galaxies 
are clustered []49[| . The magnitude of the observed excess is found to be 2.5 - 2.8 a in 
terms of Gaussian probabilities. However, such large-scale correlation with the super- 
Galactic plane was not observed in the data sets of the AGASA [ 50 ] , SUGAR [51, |52 
and Fly's Eye [^BJ experiments. 

The arrival direction of the events with energy > 4 x 10 19 eV registered by the 
AGASA experiment is shown in Fig. |J The data show no significant super-Galactic 
plane enhancement. However, some excess of particles from the general direction of 
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the super-Galactic plane is observed in the energy bin log(i?[eV]) = [19.1, 19.2], with 
/I = 0-36 ± 0.15 |53|]. 3 The highest energy events are distributed widely over the 
sky, without any apparent counterparts, such as sources in the Galactic plane or in 
the local supercluster. Moreover, the data are consistent with an isotropic distribu- 
tion of sources, in sharp contrast to the anisotropic distribution of light in the local 
supercluster [f)4| . 

Adding to the puzzle, the AGASA experiment draws attention to the fact that 
the pairing of events on the celestial sky could be occurring at higher than chance 
coincidence f2(| . Specifically, when showers with separation angle less than the angular 
resolution # min = 2.5° are paired up, AGASA finds five doublets and one triplet among 
the 58 events reported with mean energy above 4 x 10 19 eV (see Fig. |J). The chance 
probability of observing such topology in an isotropic distribution can be estimated 
by means of the Goldberg- Weiler formalism |55| . To do so, consider the solid angle 



O ~ 4.8 sr on the celestial sphere covered by AGASA to be divided into N equal 
angular bins, each with solid angle u ~ tt9 2 . Then, by tossing n events randomly into 



N ~ -— = 1045 - I —j- J (3) 
tt6 z 1 sr \1 J 

bins, one is left with a random distribution. Now, identify each event distribution by 
specifying the partition of the total sample of AGASA into a number m of empty bins, 
a number mi of single hits, a number ni2 of double hits, etc., among the N angular bins 
that constitute the whole exposure. The probability to obtain a given event topology 
is 



55 



A> n\ yr (m-r ; (4) 



N N n'° AA mi 

n—n J 



where 

m~ = N ( 

Using Stirling's approximation for the factorials with the further assumption N ^> n ^> 
1, Eq. (fD can be re- written in a quasi- Poisson form 



n 



e 

mA 
b=2 J 



(6) 



3 The super-Galactic plane enhancement factor is given by I b s (b)/I eX p (b) = (1 — /^ G ) + 
1.402 exp [— 6 2 ] . Compare this with the Fly's Eye parameterization for the Galactic plane 
given in Eq. (||). 
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FIG. 5: Inclusive probabilities for five doublets and one triplet in a 58-event sample at 
AGASA. Solid (exact), dashed (Poisson). From Ref. @. 



where r = (N — mn)/n ~ 1, and the prefactor V is given by 



V = e 



-(n-mi) 



7? 



mi 



(7) 



For "sparse events" , where N ^> n, one expects the number of singlets nil to approxi- 
mate the number of events n. In such a case the prefactor is near unity. Figure |5| shows 
the inclusive probabilities for observing 5 doublets and one triplet given 58 events at 
AGASA as a function of the angular resolution. 4 The inclusive probability is extremely 
sensitive to the angular binning. The chance probability within the experimental an- 
gular resolution of AGASA is less than 10~ 3 , in agreement with the estimates of the 
AGASA Collaboration @. 

The "world" data set has also been studied p8[. Six doublets and two triplets out 



The specified number of j-plets plus any other cluster, counts as all the j-plets + extra-clusters. 
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of 92 events with energies > 4 x 10 19 eV were found, with the chance probability 
being less than 1% in the restricted region within ±10° of the super-Galactic plane. 
Very recently, the angular two-point correlation function of a combined data sample 
of AGASA (E > 4.8 x 10 19 eV) and Yakutsk (E > 2.4 x 10 19 eV) was analyzed 
For a uniform distribution of sources, the probability of chance clustering is reported 
to be as small as 4 x 10~ 6 . The implication of such clusters would be profound, but 
the case for them is not yet proven. To calculate a meaningful statistical significance 
in such an analysis, it is important to define the search procedure a priori in order to 
ensure it is not (inadvertently) devised especially to suite the particular data set after 
having studied it In the above mentioned analysis, for instance, the angular bin 
size was not defined ahead of time. 



D. Measurement techniques 

There are several techniques which can be employed in detecting ultrahigh energy 
CRs, ranging from direct sampling of particles in the shower to measurements of associ- 
ated fluorescence, Cerenkov or radio emissions, or possibly even radar detection of the 
air shower. Direct detection of shower particles is the most commonly used method, 
and involves constructing an array of sensors spread over a large area to sample par- 
ticle densities as the shower arrives at the Earth's surface. Another well-established 
method involves measurement of the longitudinal development of the EAS by sensing 
the fluorescence light produced via interactions of the charged particles in the atmo- 
sphere. A more recently proposed technique uses radar echos from the column of 
ionized air produced by the shower. In the rest of this section, we give an overview of 
the main features of these experimental techniques, and afterwards outline the status 
of existing and pending ultrahigh energy CR experiments. More detailed and rigorous 
treatments of the current experimental situation of ultrahigh energy CRs are given in 
several review articles [^6], [61], [62], |63| . 



1. Surface arrays 

A surface array is comprised of particle detectors, such as plastic scintillators or 
Cerenkov radiators, distributed with approximately regular spacing. Such detectors 
measure the energy deposited by particles in the EAS as a function of time, and from 
energy density measured at the ground and the relative timing of hits in the different 
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FIG. 6: Average lateral distribution of simulated showers (with the program 
aires/qgsjet [ |64|| ) compared to average measurements from Volcano Ranch |p5|| . Here 
r/ry[ is the distance to the shower axis in units of the Moliere radius (at the Volcano Ranch 
elevation rjyi ~ 100 m) and S is the lateral distribution of particles at ground in units of 
minimum ionizing particles per square meter (mips/m 2 ). The estimated size of the showers 
is 10 s particles, which corresponds to an energy of about 10 18 eV. Taken fron Ref. 66]. 



detectors one can estimate the energy and direction of the primary CR. 5 

Reconstruction of air showers involves fitting the lateral distribution (LD) function 
of particle densities at the ground (see for example Fig. |5p. The exact form of the 
LD depends on how the experimental apparatus responds to the shower particles. 
The AGASA experiment, for example, employs 5 cm thick plastic scintillators. For 
energies which are typical of particles at about 1 km from the shower axis (usually 



5 An assumption of axial symmetry is generally made; this assumption is only valid for zenith angles 
9 < 60°, because at larger angles the low energy secondaries are deflected by geomagnetic fields. 
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FIG. 7: Simulation of a muon entering a tank of water from above and radiating Cerenkov 
light. Housings for three photomultiplier tubes are visible at the top of the tank. This is a 
simulation of a Pierre Auger Observatory water tank [E7[]. 



called the "core" ) these scintillators generate similar signals for the electron and muon 
components of the shower, while producing much smaller signals for the photons. The 
Haverah Park and Auger experiments, on the other hand, employ tanks of water about 
a meter deep and a meter and a half in radius, as illustrated in Fig. [?]. In this case, 
Cerenkov radiation produced by charged particles passing through the water bounces 
off the reflective tank walls and is picked up by photomultiplier tubes mounted inside. 
The signal size is thus proportional to the track length of the charged particle in the 
water. For energies typical of particles about 1 km from the core, muons tend to 
penetrate the full depth of the tank and thus produce large signals, whereas electrons 
and photons are quickly absorbed and produce signals an order of magnitude smaller. 

Furthermore, the fraction of the energy deposited by the muon and electromagnetic 
components changes with distance from the core and has been determined experimen- 
tally. For the Haverah Park experiment, for instance, the lateral density distribution 
of the water Cerenkov signal in units of fully penetrating vertical equivalent muons per 
m 2 , can be fitted with a modified power- law [|68[ 



where the distance to the shower core, r, is in meters, k is a normalization parameter, 
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(3 = 1.03 ± 0.05, and rj is given by 



77 = 3.49 - 1.29 sec 9 + 0.165 log 



E 



10 17 eV/ ' 



(9) 



with 9 the shower zenith angle, and E ^ 10 19 eV the shower energy. 

Since the response of plastic scintillator slabs to electrons, muons and photons is 
somewhat different from that of deep water Cerenkov radiators, a different form of the 
LD function is appropriate. For instance, in the case of AGASA, the empirical formula 



for the LD of charged particles with sec# < 1.7 reads [^9 
S{r 



C 



r 



1 + 



-77+a r 



1 + 



r \ - 



1 km 



) 



(10) 



where r is the distance in meters from the core, tm is the Moliere radius at two radiation 
lengths above the observation level (91.6 m at Akeno), C is a normalization factor, a 
and S are 1.2 and 0.6, respectively, 77 is the parameter which indicates the slope of the 
LD 6 for r > function of the zenith angle of the arrival direction, and is given 

by 

77 = 3.97 -1.79 (sec 9-1.0). (11) 
The primary CR energy is determined using a formula derived from simulations, 



E = 2.0 x HrS (600) eV 



(12) 



where S*o(600) is the charged particle density in 1/m 2 at 600 m from the core for a 
vertical equivalent shower. So (600) is evaluated from the observed local density at 
600 m, 

.^( sec fl_l)_^( sec £_l)2l 

Ai A 2 



5" (600) = S(600) exp 



(13) 



5941S g/cm 2 . 



The first order 



where X = 920 g/cm 2 , A x = 500 g/cm 2 , and A 2 
approximation agrees with measurements for sec^ < 1.4. The error quoted in the 
second term is at the 68% CL. 

As discussed previously, the muon content of the EAS at ground level depends on the 
composition of the primary CR. Thus, surface arrays with some ability to distinguish 
muons from electrons and photons can determine something about the primary CR 
species. Water Cerenkov detectors and scintillator-absorber combinations have such 
distinguishing power, for example. It is also possible to gauge the muon content of 
an EAS from the signal rise time, as the muonic component tends to be compressed 



6 The selection criterion of deeply developing showers is 77 of 4 to 6, which indicates a steeper lateral 
distribution than that of the average EAS (more on this below) |7Q]. 
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in time compared to the electromagnetic component. Typically experiments measure 
something like the time interval between the arrival of the 10% and 50% points of the 
integrated signal (£10/^50) • 

2. Fluorescence eyes 

As an EAS develops in the atmosphere it dissipates much of its energy by exciting and 
ionizing air molecules along its path. Excited nitrogen molecules fluoresce producing 
ultraviolet radiation. The shower development appears as a rapidly moving spot of light 
describing a great circle path across a night-sky background of starlight, atmospheric 
airglow, and man-made light pollution. The angular motion of the spot depends on 
both the distance and the orientation of the shower axis. The apparent brightness 
of the spot depends on the instantaneous number of charged particles present in the 
shower, but it is also affected by Cerenkov contamination and atmospheric scattering. 

The fluorescence trail is emitted isotropically with an intensity that is proportional 
to the number of charged particles in the shower, N e . 7 The ratio of the energy emitted 
as fluorescence light to the total energy deposited is less than 1%, hence low energy 
(< 10 17 eV) showers can hardly be observed. Furthermore, observations can only be 
done on clear moon-less nights, resulting in an average 10% duty cycle. The emitted 
light is typically in the 300 - 400 nm ultraviolet range to which the atmosphere is quite 
transparent. Under favorable atmospheric conditions EASs can be detected at distances 
as large as 20 km, about 2 attenuation lengths in a standard desert atmosphere at 
ground level. 

A fluorescence eye consists of several large light collectors (or telescopes) which im- 
age regions of the sky onto clusters of light sensing and amplification devices. The 
basic elements of a telescope are the diaphragm, which defines the telescope aper- 
ture, the spherical mirror that must be dimensioned to collect all the light entering 
the diaphragm in the acceptance angular range, and the camera which consists of an 
array of photomultiplier tubes (PMTs) positioned approximately on the mirror focal 
surface [[H]]. 8 The PMTs effectively pixelize the region of the sky covered by the tele- 
scope. The shower development is detected as a long, rather narrow sequence of hit 
PMTs. As the point-like image of the shower proceeds through an individual PMT, 
the signal rises, levels off, and falls again. The collection of vectors pointing from the 



7 We stress that N e denotes the charged multiplicity, but in practice it is often used to express the 
number of electrons and positrons since they completely dominate the total number of charged 
particles. 

8 A potential future development could involve the use of avalanche photodiodes (APDs) instead of 
PMTs. APDs have several advantages @. 
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hit PMTs defines the shower detector plane (SDP). 

The sensitivity of the detector depends primarily on the signal (S) to noise (JV) 
ratio. The signal is proportional to the diaphragm area, whereas the background (B) 
is proportional to the pixel (the area read out by a single PMT) solid angle times the 
diaphragm area, thus 

S S^d^ 

N y/B «pix 

where cfaph is the diaphragm diameter, and a P i x is the angular diameter of a pixel. 

Shower reconstruction involves first determining the geometry of the shower, and 
then reconstructing the shower's longitudinal profile. The first step of the geometrical 
reconstruction is the SDP determination from the PMT hit pattern. Next, the tube 
hit times are used to find the shower impact parameter and incident angle in the SDP; 
the resolution obtained for these parameters depends on the length of the observed 
track [f7ql . Once the geometry is known, the track can be chopped into angular bins 



and the longitudinal profile can be determined. Extracting the profile usually involves 



3-parameter fits to the Gaisser-Hillas function [74 



V V \ [(-^max-^o)/^] | V \ ' 

A — -A- ft \ vVma.x — A 



N e (X) = iV ejmax exp , X > X , (15) 

where A^ e max is the size at the maximum, X is the depth of the first observed in- 
teraction, and A =70 g/cm 2 . The integral of the longitudinal profile is a calorimetric 
measure of the total electromagnetic shower energy. A charged particle in the cas- 
cade deposits an average of 2.2 MeV into the atmosphere in each depth interval of 1 
g/cm 2 |j75|| , so the total electromagnetic energy (in MeV) is given by 

E em = 2.2 J N e (X) dX . (16) 

The largest cosmic ray air shower ever recorded has an estimated energy of (3.2 ±0.9) x 
10 20 eV, reaching the maximum size near a depth of 815 g/cm 2 [12"]. The longitudinal 



development is shown in Fig. The size at maximum is greater than 200 billion 
particles! 



3. Radio detection techniques 

In addition to the ultraviolet emission exploited by fluorescence eyes, it may also be 
possible to call into service the radio frequency (RF) energy generated by air-showers. 
CR showers may induce radio pulses through several mechanisms, though it is thought 
that from about 20-100 MHz, the dominant process can be described as coherent 
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FIG. 8: The 3-parameter best fit shower profile for the highest energy event observed by the 
Fly's Eye experiment. Published in |T^ |. 



synchrotron emission by the electron and positron pairs propagating in the Earth's 
magnetic field W%. 



RF pulses coincident with EASs were first measured in the 1960's [pf] , and these early 
results were followed up with additional work over another decade or so. Subsequently 
the technique was mostly abandoned in light of the promising results from surface 
arrays and fluorescence eyes, though quite recently a few groups have undertaken a 
new examination of the method. For example, a prototype RF detection system using 
a single antenna in conjunction with the buried muon detectors of the CASA-MIA 
array has recently been evaluated, with the goal of demonstrating the feasibility of 
incorporating such an RF receiver into an existing experiment like the Pierre Auger 
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Observatory |79|] . Though no significant EAS-induced pulses were observed, a number 
of technical issues relevant to future development were uncovered. 

More ambitious possibilities have been recently described |8(| in the context of next- 
generation digital telescopes. In particular, the planned Low-Frequency Array (LO- 
FAR) 9 , which could begin construction as early as 2004, would be well-suited to study 
RF emissions from EASs. This device will comprise around 100 stations of 100 dipoles 
distributed over a region some 400 km in radius. A cluster of computers will correlate 
data from the antennae to synthesize "virtual telescopes" which look in any desired di- 
rection. The CR energy range for which one achieves both a reasonable signal-to-noise 
ratio and a reasonable event rate depends on the number of dipoles actually employed. 
A single LOFAR-type station of about 100 dipoles, for example, would be useful for 
measurements in the range 10 15 -10 17 eV; indeed just such a prototype station, known 
as LOPES |yj, is scheduled to be operational in 2003 at the KASCADE experiment. 
Falcke and Gorham |50j estimate that LOFAR should be able to observe events up to 
10 20 eV at a rate of order 1 per year. 

In addition to such passive radio techniques, it may also be possible to detect radar 
reflections off the ionization columns produced by EASs. This idea was suggested al- 
ready in 1940 [82], and recently the issue has been re-explored [83, 84] as either an 
independent method to study air showers, or as a compliment to existing fluorescence 
and surface detectors. For example, one could use a fluorescence detector to trigger a 
radar system to interrogate in the approximate direction of the shower. It has been 
estimated that showers from primaries of order 10 18 eV should be detectable f83|, and 



that range estimates with precision around 20 m may be possible. The technique would 
work best for distant or horizontal showers, the latter being of particular interest in 
searches for deeply penetrating showers induced by neutrino primaries. One important 
uncertainty is the lifetime of the ionization trail, which decays away due to electron 
recombination and attachment and can thus limit the time available for radar interro- 
gation; these effects depend on altitude and atmospheric conditions. A proposal has 
recently been put forth to evaluate the method using the Jicamarca radar system near 



Lima, Peru [85 



E. CR-observatories: past and present 

The first measurements of ultrahigh energy CRs were carried out by Linsley at 
Volcano Ranch (35°09'N, 106°47'W) in the late 1950's f8|, 0. The detector comprised 
an array of 19 scintillation counters each of 3.3 m 2 , spaced on a 884 m hexagonal grid. 



9 tittp: //www. lofar.org 
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More recently, experiments have been performed at Haverah Park in England (53°58'N, 
1°38'W), Yakutsk in Russia (62°N, 130°E), Sydney in Australia (30°32'S, 149°43'E), 
Dugway in the Utah desert (40°N, 112°W), and near the village of Akeno, about 100 km 
west of Tokyo (38°47'N, 138°30'E). 

The Haverah Park detector || ceased operation in July 1987 after a life of 19 years. 
The detector array consisted of deep (1.2 m) water Cerenkov detectors with stations 
ranging in area from 1 to 54 m 2 (typically 34 m 2 ) distributed over an area of approxi- 
mately 12 km 2 . Restrictions on land access made it impossible to position the detectors 
on a uniform grid. Therefore, a 4 detector array with 500 m spacing was deployed at 
at the center and 6 sub-arrays of 50 m and 150 m spacing were located at about 2 km 
from the center. In the ongoing analyses of the data, primary energies are derived from 
the water Cerenkov detector density at 600 m from the shower axis, p(600). Random 



errors in the determination of p(600) are less than 30% ||88|| . The energy threshold 
was approximately 6 x 10 16 eV. The angular resolution of the array was measured to 
be: for zenith angle (9), 2.5° x sec# with 9 e (0°, 75°), whereas for azimuth angle ((f)), 
2.5 x cosec# with 9 G (15°, 90°) [g|]. 

The Yakutsk array, in operation since 1972, was expanded in 1974 for sensitivity to 
ultrahigh energy CRs, and spread over a ground area of approximately 18 km 2 . Most 
of the detectors are plastic scintillators arranged on a triangular grid. The spacing of 
the detectors was around 1 km, with the center of the array filled with 7 detectors on 



a 500 m triangular grid ||90|| . Underground stations measured the muon flux, with an 
energy threshold of 1 GeV. The primary energy is estimated from the particle density 
at 600 m from the shower core. In 1995, the array was re-arranged into a smaller area 
(~ 10 km 2 ) in order to investigate in detail the spectrum around 10 19 eV ||9 1|| . 

The Sydney University Giant Air-shower Recorder (SUGAR) [0 operated between 
1968 and 1979, and is the largest array to date in the Southern hemisphere. 47 detector 
stations were arranged to cover an area of approximately 70 km 2 . Each of these stations 
consisted of two 6 m 2 conical liquid scintillator tanks separated by 50 m. The angular 
resolution of the array is quoted [|7| to be 3° x sec 9 for showers detected by more 
than three stations. However, most of the events were only viewed by 3 stations, and 
consequently the resolution is closer to 6° for vertical showers f92fl . 

The Akeno Giant Air Shower Array (AGASA) |93| , |9~4|j consists of 111 scintillation 
detectors each of area 2.2 m 2 , spread over an area of 100 km 2 with 1 km spacing. This 
gives an acceptance of 125 km 2 sr for CRs above 10 19 eV. 10 The array detectors are 



The array has been divided into 4 branches for the purpose of trigger and data acquisition. These 
branches are called the "Akeno Branch" , the "Sudama Branch" , the "Takane Branch" , and the 
"Nagasaka Branch". The 20 km 2 array f95j, operational since the end of 1984, forms a part of the 
"Akeno Branch" . 
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connected and controlled through a sophisticated optical fiber network. The array also 
contains a number of shielded scintillation detectors which provide information about 
the muon content of the showers. The full AGASA experiment has been running since 
1992. 

The fluorescence technique has so far been implemented only in the Dugway desert by 
the Fly's Eye group from the University of Utah. Following a successful trial at Volcano 
Ranch |9B| the Utah group built a device containing two separated Fly's Eyes made 
up of 880 and 460 PMTs, 3.3 km apart [57, 55]. The two-eye configuration monitored 
the sky from 1986 until 1993. The first detector, Fly's Eye I, had full operation since 
1981. It was made of 67 telescopes of dd p h = 1-5 rn and spherical curvature, each with 
12 or 14 PMTs at the focus. The mirrors were arranged so that the entire night sky 
was imaged, with each PMT viewing a hexagonal region of the sky with a pix = 5°. 
Fly's Eye II was made of of 36 mirrors of the same design. This detector viewed only 
half of the night sky. 

As an up-scaled version of Fly's Eye, the High Resolution (HiRes) Fly's Eye detec- 
tor has recently begun operations J5?|. It uses 14 spherical telescopes of d^ p h = 2.0 m 
to collect the light from a 0.95 sr portion of the sky. The image plane of each tele- 
scope is populated with an array of 256 hexagonal PMTs, yielding a P i x = 1°. In 
monocular mode, the effective acceptance of this instrument is ~ 350(1000) km 2 sr at 
10 19 (10 20 ) eV, on average about 6 times the Fly's Eye acceptance, and the threshold 
energy is 10 17 eV. This takes into account a duty cycle of about 10% typical of the 
fluorescence technique. The field of view of the telescopes is centered on the Chicago 
Air Shower Array (CASA) and the Michigan Muon Array (MIA) ||100|| , situated 3.3 km 
to the northeast. The combination of these instruments has been used as a prototype 
hybrid detector in which HiRes records the development profile, CASA records the 
ground particle density and MIA detects the muonic component of a common shower. 



F. The potential of the Pierre Auger Observatory 



The Pierre Auger Observatory (PAO) [|101| , |102|| is designed to work in a hybrid mode 



employing fluorescence detectors overlooking a ground array of deep water Cerenkov 
radiators. During clear, dark nights, events will be simultaneously observed by fluores- 
cence light and particle detectors at ground level. The PAO is expected to measure the 
energy, arrival direction and primary species with unprecedented statistical precision. 
It will eventually consist of two sites, one in the Northern hemisphere and one in the 
Southern, each covering an area of 3000 km 2 and consisting of 1600 particle detectors 
overlooked by 4 fluorescence detectors. The overall acceptance (2 sites) is 14000 km 2 
sr. A prototype engineering array, l/40th-scale, is (at the time of writing) already 
taking data in Malargiie, Argentina (35°12'S, 69°12'W) fllO^ . 
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The surface array stations are cylindrical water Cerenkov detectors with area 10 m 2 , 
spaced 1.5 km from each other in a hexagonal grid. Cerenkov radiation emitted by 
charged particles penetrating the detector is read out by 3 PMTs. The output signal 
is digitised by flash ADCs, with the aim of separating the electromagnetic signal (low 
energy electrons and photons) from the muons crossing the tank. Event timing is made 
possible via global positioning system (GPS) satellites with a precision of a few tens 
of ns. Communication between the stations is achieved using radio signals by methods 
similar to cellular telephone techniques. The stations are powered by solar panels and 
batteries which allow autonomous operation. 

The configuration of the fluorescence detectors is arranged to optimize the hybrid 
detector performance. The number of eyes and their location on the site are chosen 
so that all showers of energy > 10 19 eV that hit the surface array are seen by at least 
one eye. A further constraint comes from the need to limit the systematic error in the 
measurements deriving from uncertainty in the attenuation length of the atmosphere 
traversed by the light in its path from shower to detector. The optimal configuration 
was determined by Monte Carlo simulation, guided by the orographic constraints on the 
site. The base-line design of the detector includes 4 fluorescence eyes, each comprised 



of six telescopes with dd p h = 1-7 m and a P j x = 1.5° ||104 . 

The angular and energy resolution of the ground array (without coincident fluo- 
rescence data) are typically less than 1.5° and less than 20%, respectively. "Golden 
events," events detected by both methods simultaneously, will have a directional re- 
construction resolution of about 0.3° for energies near 10 20 eV. If an event trigger is 
assumed to require 5 detectors above threshold, the array is fully efficient at 10 19 eV. 
In three years of running, the surface arrays in both hemispheres, operating 24 hours 
per day, will collect more than 1000 showers above 4 x 10 19 eV with approximately 
uniform sky exposure. This will enable a straightforward search for correlations with 
discrete sources and also a sensitive large scale anisotropy analysis. 

In addition, PAO offers a window for neutrino astronomy above 10 17 eV. For standard 
neutrino interactions in the atmosphere, each site of PAO reaches ~ 15 km 3 w.e. sr 
of target mass around 10 19 eV [|105|| , which is comparable to other neutrino detectors 
being planned. 11 An even greater acceptance | 106|| should be achievable for the case of 
Earth-skimming neutrinos which produce a r, as we discuss in Sec.II-I. 



w.e. = water equivalent. 
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G. The Telescope Array Project 



The Telescope Array will comprise a collection of fluorescence detectors dispersed 
over a large area near Salt Lake City, Utah ||107|| . Ten observational stations separated 
from each other by 30 — 40 km are planned, each station containing 40 telescopes of 
ddph = 3 m. An imaging camera with 256 PMTs (a P i x = 1°) will be installed on the 
focal plane of each telescope. The effective aperture of the array will be approximately 
5000 km 2 sr for 10 20 eV particles assuming a 10% observation duty factor. This aperture 
is around 30 times larger than the existing AGASA ground array. The energy, arrival 
direction, and shower maximum will be determined with an accuracy of 6%, 0.6°, 
and 20 g/cm 2 , respectively. In addition, the Telescope Array will observe high energy 
gamma rays from point sources in the sub-TeV energy region. Gamma rays will be 
distinguished from the large hadronic background using the imaging patterns observed 
with many telescopes. 



H. 200?: A space odyssey 

Recently, NASA initiated a concept study for space-based detectors which will stereo- 
scopically image, from equatorial orbit, the nitrogen fluorescence light generated by 
EASs induced by ultrahigh energy (> few x 10 19 eV) CRs. The Orbiting Wide-angle 
Light-collectors (OWL) |108|1 mission will involve photodetectors mounted on 2 satel- 



lites orbiting at 640 km above the Earth's surface. The eyes of the OWL will monitor 
a large atmospheric volume and record the ultraviolet fluorescent trails with a time 
resolution of ~ 1 //s or less in segmented focal plane arrays. The segmentation of the 
arrays will be such as to sample the atmosphere near the Earth's surface in ~ 1 km 2 
pixels. The use of a space-based platform enables an extremely large event acceptance, 
allowing a high statistics measurement. OWL is set for possible implementation after 
2007. 

Two different baseline instruments have been proposed to achieve wide-field of views. 
The first is a refractive design using two Fresnel lenses which focus onto a large focal 
plane array ||109 |. The second design uses Fresnel correcting optics which focus onto 



a spherical reflector in a Maksutov design that in turn focuses the light onto a focal 



plane array [110, [1 llj. For a satellite separation of 500 (2000) km, the instantaneous 
acceptance is 1.5 x 10 6 (3.75 x 10 6 ) km 2 sr. Using a 10% duty factor, the 500 (2000) km 
configuration leads to an effective acceptance of 1.5 x 10 5 (3.75 x 10 5 ) km 2 sr. There- 
fore, assuming a continuation of the CR-spectrum oc E~ 2 ' 75 , one expects rates of 1500 
events/yr (500 km) and 3750 events/yr (2000 km) for E > 10 20 eV JTE§. The effective 



acceptance for the Maksutov baseline assuming 1000 km orbits and 500 (2000) km is 
2.0 x 10 5 (4.0 x 10 5 ) km 2 sr which leads to an event rate of 2000 (4000) events per year 
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for E > 10 20 , again assuming an E 2,75 spectrum continuation ||112j . 

For the near future (~ 2006), the European Space Agency is studying the feasibil- 
ity of placing a single eye on the International Space Station, which will serve as a 
pathfinder mission to develop the required technology to observe the fluorescent trails 
of EASs. The smaller viewing volume of the Extreme Universe Space Observatory 
(EUSO) ||113|| will result in a smaller event rate by a factor of ~ 5 compared to OWL. 



I. Ultrahigh energy neutrino experiments 

Up to now, we have discussed experiments designed primarily to study air showers 
initiated by hadrons or photons. As pointed out in Sec.II-B, however, neutrinos may 
also induce extensive air showers, so current and future air shower experiments might 
also function as neutrino detectors. In this section, we briefly review dedicated neutrino 
detectors as well as neutrino searches carried out by existing air shower experiments. 

The traditional technique to observe high energy neutrinos involves detecting the 
optical Cerenkov light emitted by muons produced in charged current interactions of 
neutrinos with nucleons either in ice or water. The largest pilot experiments (~ 0.1 km 
in size) are: the now defunct DUMAND (Deep Underwater Muon and Neutrino Detec- 
tor) experiment ||114|| in the deep sea near Hawaii, the underwater experiment in Lake 



Baikal ||115|| , and AMANDA (Antarctic Muon And Neutrino Detector Array ) |lT6| in 



the South Pole ice. Next generation neutrino telescopes aim towards an active volume 
in the range of 1 km 3 of water. Projects under construction or in the proposal stage are: 
two deep sea experiments in the Mediterranean, the French ANTARES (Astronomy 
with a Neutrino telescope Abyss environment RESearch) [ |1171| and NESTOR (Neutrino 
Experiment SouthwesT Of GReece |118| ), and ICECUBE ||119|1 , a scaled up version of 



the AMANDA detector. 

Above 10 17 eV the neutrino interaction length is below 2000 km w.e. in rock, and so 
upward going neutrinos are typically blocked by the Earth. This shadowing severely 
restricts the high energy event rates in underground detectors. Current limits on high 
energy (10 15 — 10 17 eV) neutrino fluxes come from measurements of the Extensive Air 
Shower array on the TOP of the underground Gran Sasso Laboratory in central Italy 
(EAS-TOP) fng, [m|, and the Frejus detector JT^, g, located in an underground 



laboratory near the middle of the road tunnel connecting Modane (France) and Bar- 
donecchia (Italy) in the Alps. 

Large ground arrays or fluorescence eyes, for which the interaction medium is not the 
Earth but the atmosphere, are complementary to neutrino telescopes in the analysis 
of the energy spectrum. Specifically, the neutrino cross section at ultrahigh energies 
is non-negligible (about 10~ 32 cm 2 at 10 18 eV ||1 24|| ) , and so quasi- horizontal (9 > 
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75°) neutrinos that traverse an atmospheric depth of up to 360 m w.e. can initiate 
a shower. As discussed earlier, hadronic showers at large zenith angles have their 
electromagnetic component extinguished, and only high energy muons created in the 
first stages of the shower development survive. Therefore, the shape of the shower 
front is relatively flat. The Fly's Eye Collaboration searched for deeply developing air 
showers (DDASs) and upward-moving showers above 10 17 eV and no unusual events 
were found in 10 6 s of running time 125|| . Additionally, the AG AS A group searched for 
giant deeply developing air showers [70|. During the observation live time of 9.7 x 10 7 s 
no candidate DDAS was found. From this, the AGASA Collaboration concluded that 
the 90% CL upper limit on the flux of DDASs with energy > 10 19 5 eV is 1.9 x 10" 16 m" 2 
s -i sr — i ? a factor of 10 lower than the flux of cosmic rays above 10 19,5 eV. The current 
upper bound (E u ~ 10 19 5 eV) on the total cosmic neutrino flux, obtained by combining 
the exposures of the AGASA and Fly's Eye experiments (the latter integrated over all 
its operating epochs), is found to be 5.1 x 10 -12 m -2 s _1 sr -1 , at the 95% CL ||126 |. 
The sensitivity of surface arrays and fluorescence eyes could be significantly enhanced 
by triggering on neutrinos that skim the Earth, traveling at low angles along chords 
with lengths of order their interaction length [ 127 , |128| , [129| , |130| , |131| , |132| , |133|| . Some 
of these Earth-skimming neutrinos may be converted into tau leptons in the Earth's 
crust. Unlike electrons, which do not escape from rocks, or muons, that do not produce 
any visible signal in the atmosphere, taus can escape even from inside the rock and 
produce a clear signal if they decay above the detector, increasing the ^-event rate. 

The study of radio pulses from electromagnetic showers created by neutrino in- 
teractions in ice would provide an increase in the effective area up to 10 4 km 2 . A 
prototype of this technique is the Radio Ice Cerenkov Experiment (RICE) |[L34|| . Sim- 
ilar concepts are used by the Goldstone Lunar ultrahigh energy neutrino Experiment 
(GLUE) to set an upper bound on the ultrahigh energy neutrino flux ||135| , |136|| . In this 
experiment, the non-observation of microwave Cerenkov pulses from electromagnetic 
showers induced by neutrinos interacting in the moon's rim leads to flux upper limits: 
log(£ 2 J) = -3.03, -2.66, and -2.30 GeV cm" 2 s" 1 sr" 1 , for 1.0 x 10 22 , 3.0 x 10 22 , 
and 1.0 x 10 23 eV, respectively. 

For a more extensive discussion on the status of ultrahigh energy neutrino experi- 
ments the reader is referred to yi37[|. 



III. COSMIC RAY ACCELERATION 

It is most likely that the bulk of the cosmic radiation is a result of some very general 
magneto-hydrodynamic (MHD) phenomenon in space which transfers kinetic or mag- 
netic energy into cosmic ray energy. The details of the acceleration process and the 
maximum attainable energy depend on the particular physical situation under consider- 
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ation. There are basically two types of mechanism that one might invoke for bottom-up 
CR production. 12 The first type assumes the particles are accelerated directly to high 



energy by an extended electric field ||138|| . This idea can be traced back to the early 
1930's when Swann ||139|| pointed out that betatron acceleration may take place in the 
increasing magnetic field of a sunspot. These so-called "one-shot" mechanisms have 
been worked out in greatest detail, and the electric field in question is now generally 
associated with the rapid rotation of small, highly magnetized objects such as neutron 
stars (pulsars) or active galactic nuclei (AGN). Electric field acceleration has the ad- 
vantage of being fast, but suffers from the circumstance that the acceleration occurs 
in astrophysical sites of very high energy density, where new opportunities for energy 
loss exist. Moreover, it is usually not obvious how to obtain the observed power law 
spectrum in a natural way, and so this kind of mechanism is not widely favored these 
days. The second type of process is often referred to as statistical acceleration, because 
particles gain energy gradually by numerous encounters with moving magnetized plas- 
mas. These kinds of models were mostly pioneered by Fermi ||140 1 . In this case the E~ 2 



spectrum very convincingly emerges. 13 However, the process of acceleration is slow, 
and it is hard to keep the particles confined within the Fermi engine. In this section we 
first provide a summary of statistical acceleration based on the simplified version given 



in Ref |HJ (See also [133]). For a more detailed and rigorous discussion the reader is 
referred to ||143|| . After reviewing statistical acceleration, we turn to the issue of the 
maximum achievable energy within diffuse shock acceleration and explore the viability 
of some proposed ultrahigh energy CR-sources. 



A. The Fermi mechanism 



In his original analysis, Fermi ||140|| considered the scattering of CRs on moving 
magnetized clouds. In this CR entering into a single cloud with energy Ei and 

incident angle 9i with the cloud's direction undergoes diffuse scattering on the irreg- 
ularities in the magnetic field. After diffusing inside the cloud, the particle's average 
motion coincides with that of the gas cloud. The energy gain by the particle, which 
emerges at an angle Of with energy Ef, can be obtained by applying Lorentz transfor- 
mations between the laboratory frame (unprimed) and the cloud frame (primed). In 



12 In bottom-up models the CR particles start with low energy and are accelerated, as opposed to 
top-down models where usually exotic particles start initially with very high energy and cascade 
decay to the CR particles. 

13 A point worth noting at this juncture: A power law spectrum does not necessarily point to the 
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the rest frame of the moving cloud, the CR particle has a total initial energy 



£[ = T £7, (1-/3 cos 0i), (17) 

where T and (3 = V/c are the Lorentz factor and velocity of the cloud in units of the 
speed of light, respectively. In the frame of the cloud we expect no change in energy 
[E[ = E'j), because all the scatterings inside the cloud are due only to motion in the 
magnetic field (so-called collisionless scattering). There is elastic scattering between the 
CR and the cloud as a whole, which is much more massive than the CR. Transforming 
to the laboratory frame we find that the energy of the particle after its encounter with 
the cloud is 

E f = TE' f (l + (3 cos9 f ). (18) 
The fractional energy change in the laboratory frame is then 



AE Ef-Ei 1-/9 cos Oi + P cos 6' } - (3 2 cos 0* cos 8' f 
~E~ = ~E~ = 1-/3 2 



1. (19) 



Inside the cloud the CR direction becomes randomized and so (cos 6^) = 0. The average 
value of cos 0j depends on the relative velocity between the cloud and the particle. The 
probability P per unit solid angle Q of having a collision at angle 0j is proportional to 
(v — V cos 9i), where v is the CR speed. In the ultrarelativistic limit, i.e., v ~ c (as 
seen in the laboratory frame), 



so 



dP 

' oc (1-/3 cos 0<)> (20) 



(cos0 t ) = -^. (21) 



Now, inserting Eq. fl2T|) into Eq.flH]) one obtains for f3 <C 1, 

(A£) = l + / 373 4 

e 1-/32 ~ 3 p • y ll > 

Note that (AE)/E cx /3 2 , so even though the average magnetic field may vanish, there 
can still be a net transfer of the macroscopic kinetic energy from the moving cloud to 
the particle. However, the average energy gain is very small, because (3 2 1. 

A version of Fermi's mechanism which is first order in /3, and thus a more efficient 
accelerator, is realized for CR encounters with plane shock fronts ||144j| . In this case, a 
large shock wave propagates with velocity —u[, as depicted in Fig. [5| Relative to the 
shock front, the downstream shocked gas is receding with velocity 112, where | w 2 1 < 
and thus in the laboratory frame it is moving in the direction of the front with velocity 
V = —u\ + u 2 . In order to find the energy gain per shock crossing, we identify the 
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upstream 



downstream 



FIG. 9: CR acceleration at a shock front. A planar shock wave is moving with velocity —u\. 
A CR particle is repeatedly crossing the front and scattering in magnetic irregularities. 



magnetic irregularities on either side of the shock as the clouds of magnetized plasma 
in the Fermi mechanism we discussed previously. By considering the rate at which CRs 
cross the shock from downstream to upstream, and upstream to downstream, one finds 
(cos6>,,) = —2/3 and (cosO'j) = 2/3. Hence, Eq. ( [L9|) can be re-written as 

m * p = l^i, (23) 

h 3 3c 

Note this is first order in f3 — V/c, and is therefore more efficient than Fermi's original 
mechanism. This is because of the converging flow - whichever side of the shock you 
are on, if you are moving with the plasma, the plasma on the other side of the shock 
is approaching you. 

The rate at which CRs cross from upstream to downstream is given by the projection 
of the isotropic CR flux onto the plane shock front 

1 a( a\ aa HcrV a UcrV (oa\ 
d(coso) / dcp cosy ~ , (24) 

o Jo 4?r 4 
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where ncR is the number density of particles undergoing acceleration. On the other 
hand, the rate of convection downstream away from the shock front is 

n oss = n C RU 2 , (25) 

so the probability of crossing the shock once and then escaping from the shock (being 
lost downstream) is given by 

Prob. (escape) = loss ~ 4 — . (26) 

Across f 

The probability of returning to the shock after crossing from upstream to downstream 
is 

Prob. (return) = 1 — Prob. (escape) , (27) 

and so the probability of returning to the shock n times and also of crossing the shock 
at least n times is 

Prob. (cross > n) = [1 — Prob. (escape)]™. (28) 
Therefore, the energy after n shock crossings is 

where E is the overall initial energy. With this in mind, the number of encounters 
needed to reach an energy E reads 

n - HE/Eo) (30) 
U ~ ln(l + (AE)/E) ' m 

and the number of particles accelerated to energies greater than E is 

Q(> E) a ±[1 - P r ob.(e S cape)]- = l - Z^^l . (31 ) 

Prob. (escape) 

m=n 

Substitution of Eq. (|30D into Eq. (|3lD leads to 

Q(> E) oc - - / ( ^] 7 (32) 

^ V ; Prob. (escape) \E J K ' 

with 

ln[l — Prob. (escape)] _1 
7= H1 + (AE)/E) ■ {dS} 
All in all, Fermi's mechanism yields the desired power law spectrum for CR acceleration. 
Note that in the first order mechanism the spectral index, 7, is independent of the 
absolute magnitude of the velocity of the plasma, and depends only on the ratio of the 
upstream and downstream velocities. 
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B. Desperately seeking Zevatrons 



1. Plausible sources 



Now we focus attention on the direct identification of candidate CR-sources. For 
well- researched reviews on this topic see ||145| , |146| |. A variety of astrophysical objects 
have been proposed to account for the origin of high energy CRs. The options include: 



Supernovae explosions ||147| , |148|| . 



Large scale Galactic wind termination shocks [|149 . 



Pulsars (neutron stars) ||150|| . 



Active galactic nuclei (AGNs) ||151 



BL Lacertae (BL Lac) - a sub-class of AGN [152, 153 



Spinning supermassive black holes associated with presently inactive quasar rem- 
nants PI [T55H. 14 



Large scale motions and the related shock waves resulting from structure forma- 
tion in the Universe | |157|| such as accretion flow onto galaxy clusters and cluster 
mergers [|158| , |159|| . 15 

Relativistic jets and "hot-spots" produced by powerful radiogalaxies. Relativistic 
jets occurring in very powerful radiogalaxies carry large amounts of energy up to 
the radio hot-spots situated far (~ 100 kpc) from the central engine. These hot- 
spots are believed to harbor strong, mildly relativistic shocks dissipating the jet 
bulk kinetic energy into heating plasma, thus generating magnetic fields which 

Additionally, a 



efficiently accelerate CRs to ultrahigh energies ||161| , |162| , IToc 
velocity shear layer at the relativistic jet side boundary can play an active role 
in the acceleration process ||164| , |165| . 



The electrostatic polarization fields that arise in plasmoids produced in planetoid 
impacts onto neutron star magnetospheres ||166| . 



Magnetars - pulsars with dipole magnetic fields approaching 
169|| - appear also as serious candidates ||170| , |171| . 



10 15 G HTB7L H6&1 



14 An unexpected apparent correlation between the highest energy events reported by the AGASA 
Collaboration and dead quasars was recently reported [|156| . 
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The famous AGASA triplet points towards the merger galaxies Arp 299 (NGC 3690 + IC 694), at 



a distance of 70 Mpc [16C] 
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Starburst galaxies ||172|, |173|, \L74\ . 



Magnetohydrodynamic (MHD) winds of newly formed strongly magnetized neu- 
tron stars [IT7 



Gamma ray burst (GRB) fireballs ||176| , 177 , 178 , |179 . 

Strangelets, stable lumps of quark matter, accelerated in astrophysical environ- 



ments 180 



Hostile aliens with a big CR gun Hi 8 HI . 



In general, the maximum attainable energy of Fermi's mechanism is determined by 
the time scale over which particles are able to interact with the plasma. Sometimes 
the acceleration region itself only exists for a limited period of time; for example, 

In such a case, 



supernovae shock waves dissipate after about 10 yr [182 



Eq. 



would have to be modifed accordingly. Otherwise, if the plasma disturbances persist for 
much longer periods, the maximum energy may be limited by an increased likelihood 
of escape from the region. The latter case is relevant to ultrahigh energies, because 
when the Larmor radius of the particle (of charge Ze), 



' L 



110 



R 



20 



ZB 



kpc 



(34) 



approaches the accelerator size it becomes very difficult to confine the CR magnetically 
to the acceleration region. Here, B^q is the magnetic field in units of fiG, and R20 = 
E/10 20 eV. If one includes the effect of the characteristic velocity (3c of the magnetic 



scattering centers, the above argument leads to the general condition ||138| , 



2f3cZeBr L 



(35) 



for the maximum energy acquired by a particle traveling in a medium with magnetic 
field B. This is sometimes called the "Hillas criterion." The upper limit on the energy 
of one-shot acceleration scenarii turns out to be quite similar to the shock acceleration 
case of Eq. (|35f) . For instance, a dimensional analysis suggests that the maximum 
energy that can be obtained from a pulsar is ||138| 



R 



= -ZeB r 2 

max ^-"^ J -'s 1 us 5 

C 



(36) 



where u is the pulsar angular velocity, B s the surface magnetic field and r ns the neutron 
star radius. Therefore, if B s ~ 10 12 G, r ns ~ 10 km, and uj ~ 607r s _1 (as for the Crab 
pulsar), a circuit connected between pole and equator would see an emf ~ 10 18 V for 
an aligned or oblique dipole. When realistic models of acceleration are constructed, 
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FIG. 10: The Hillas diagram showing size and magnetic field strengths of possible sites of 
particle acceleration. Objects below the diagonal lines (from top to bottom), derived from 
Eq. (35) assuming the extreme value f3 = 1, cannot accelerate protons above 10 21 eV, above 
10 20 eV and iron nuclei above 10 20 eV, respectively. (This version of the picture is courtesy 
of Murat Boratav). 
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however, this ideal dimensional limit is not realized, because the large potential drop 
along the magnetic field lines is significantly short-circuited by electron and positrons 
moving in the opposite directions along the field lines ||150||. 



The dimensional arguments of Eqs. (^) and (pq) are usually summarized in the form 
of the popular "Hillas diagram" ||138|| shown in Fig. [IUJ. Clearly, very few sites appear 
able to generate particle energies > 10 20 eV; either this occurs on highly condensed 
objects with huge B or enormously extended objects. 16 From a glance at Fig. ITU] it 
seems that the structures associated with active galaxies, neutron stars and GRBs have 
sufficient size and field strength to be considered as potential sources. In subsequent 
sections we discuss these potential sources in more detail. 



2. Radiogalaxies 

Fanaroff-Riley II (FRII) galaxies |184| are the largest known dissipative objects (non- 



thermal sources) in the Universe. Localized regions of intense synchrotron emission, 
known as "hot spots" , are observed within their lobes. These regions are presumably 
produced when the bulk kinetic energy of the jets ejected by a central active nucleus 
(supermassive black hole + accretion disk) is reconverted into relativistic particles 



and turbulent fields at a "working surface" in the head of the jets ||185|| . Specifically, 
the speed with which the head of a jet advances into the intergalactic medium of 
particle density n e can be obtained by balancing the momentum flux in the jet against 
the momentum flux of the surrounding medium. Measured in the frame comoving with 
the advancing head, ~ Vj [1 + (rie/nj) 1 / 2 ] -1 , where rij and v- } are the particle density 
and the velocity of the jet flow, respectively. Clearly, Vj > t>h for n e >n- 3 , in such a way 
that the jet will decelerate. The result is the formation of a strong collisionless shock, 



which is responsible for particle reacceleration and magnetic field amplification ||186 
The acceleration of particles up to ultrarelativistic energies in the hot spots is the result 
of repeated scattering back and forth across the shock front, similar to that discussed 
in Sec.III-A. The particle deflection in this mechanism is produced by Alfven waves in 
the turbulent magnetic field. This process has been studied in detail by Biermann and 



Strittmatter | 187|| . Dimensional arguments suggest that the energy density per unit of 



wave number of MHD turbulence is of the Kolmogorov type [|188]|, and so for strong 



shocks the acceleration time for protons is [189] 



40 1 I f E\ 



ir c j3h u \eB J 



1/3 



- 



R- 2/3 (37) 
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For a comprehensive discussion on the electrodynamical limitations of CR sources the reader is 



referred to the recent study in Ref. [183 
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where /3j e t is the jet velocity in units of c, u is the ratio of turbulent to ambient magnetic 
energy density in the region of the shock (of radius R), and B is the total magnetic 
field strength. The acceleration process will be efficient as long as the energy losses 
by synchrotron radiation and photon-proton interactions do not become dominant. 
The subtleties surrounding the conversion of a particle kinetic energy into radiation 
provide ample material for discussion ||187| , |19C1| , |191| , |192| , |193| , |194|| . The most pop- 



ular mechanism to date relates 7-ray emission to the development of electromagnetic 
cascades triggered by secondary photomeson products that cool instantaneously via 
synchrotron radiation [|1871 , |190| , |191| , |192| , |193|| . The synchrotron loss time for protons 
is given by ||195| 

Ts y n ~ 2-r P2 ' ( 38 ) 

ctt vn e 1 B A 

where m e , m p , o"t and T are the electron mass, proton mass, Thomson cross section, 
and Lorentz factor, respectively. The characteristic single photon energy in synchrotron 
radiation emitted by an electron is 

{3\ 1/2 heE 2 B 99 
E ^ (2 J ~ 5.4 x 10- B„ El TeV . (39) 

For a proton this number is (m p /m e ) 3 ~ 6 x 10 9 times smaller. Thus, it is evident 
that high energy gamma ray production through proton synchrotron radiation requires 
very large (0(100 G)) magnetic fields. Considering an average cross section a lp for the 
three dominant pion-producing interactions [|19(3|| . 

7P — >■ pre , (40) 

7P — > mr + , (41) 

7P — > p7T + TT~ , (42) 

the time scale of the energy losses, including synchrotron and photon interaction losses, 
reads |Lg/J 

6-7T mi c 3 r vvn 

T] ~ E E = y A3) 

loss " a T ml B* (1 + Aa) 1 + Aa> 1 ' 

where a stands for the ratio of photon to magnetic energy densities and A gives a 
measure of the relative strength of 7p interactions versus the synchrotron emission. 
Note that channel (f4*l"D involves the creation of ultrarelativistic neutrons (but T n ^ Y p ) 
with mean free path in the observer rest frame given by A„ = T n cr n , where r n ~ 900 s, 
is the neutron lifetime. Since X n > X p for T n < T p max , such neutrons can readily escape 
the system, thereby modifying the high end of the proton spectrum. Biermann and 
Strittmatter ||187|| have estimated that A pa 200, almost independently of the source 
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parameters. The most energetic protons injected in the intergalactic medium will have 



an energy that can be obtained by balancing the energy gains and losses ||197 



R 



20 



1.4 x 10 l 



jet 



■3/4 



(44) 



where i?k pc = -R/lkpc. 

As noted in the previous section, the canonical theory of diffuse shock acceleration 
not only assumes an infinite time, but also an infinite space for the particles to diffuse 
in the medium around the shock front. Therefore, in order to ascertain the capability 
of FRII to accelerate particles to ultrahigh energies, one also has to apply the Hillas cri- 



terion [|138|1 for localizing the Fermi engine in space, namely that the Larmor radius be 
less than the size of the magnetic region. For typical hot-spot conditions (B ~ 300 fiG, 
u ~ 0.5, and (5- iet ~ 0.3) and assuming that the magnetic field of the hot spot is limited 
to the observable region, one obtains E < 5 x 10 20 eV for a < 0.1 ||161||. Therefore, 



one concludes that protons can be accelerated to the highest observed energies in typ- 
ical FRII hot spots. Moreover, the shock structure in hot spots is likely to be much 
more extended than the visible region in the nonthermal radioemission, as suggested 



by magnetohydrodynamical modeling |161 



Particles can also attain ultrahigh energies (E > 10 20 eV) within the jets or the 
AGNs themselves. For instance, the knot A in the M87 jet, with a length scale Z§7 ~ 
2 x 10 20 cm, has a magnetic field strength ~ 300 /iG [ |198|| . Typical AGN sizes are 
^agn ~ 10 15 cm, and B AG ^ 
jets a < 1, whereas a 



1 G ||199| . Observational evidence suggests that in the 
1 for AGNs 



In Sec.V the ideas developed in this section will be applied to specific models. 



3. j-ray burst fireballs 

Gamma ray bursts (GRBs) are flashes of high energy radiation that can be brighter, 
during their brief existence, than any other gamma ray source in the sky. The bursts 
present an amazing variety of temporal profiles, spectra, and timescales that have puz- 
zled astrophysicists for almost three decades |200|| . In recent years, our observational 



insight of this phenomenon has been dramatically increased by the huge amount of 
data collected by the Burst and Transient Source Experiment (BATSE) on board the 
Compton Gamma Ray Observatory (CGRO), a satellite launched by NASA in 1991. 
In 9 years of operation, BATSE has accumulated a database of more than 2000 obser- 
vations. 

The temporal distribution of the bursts is one of the most striking signatures of the 
GRB phenomenon. There are at least four classes of distributions, from single-peaked 
bursts, including the fast rise and exponential decaying (FREDs) and their inverse 
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(anti-FREDs), to chaotic structures ||201| , [202|| . There are well separated episodes of 
emission, as well as bursts with extremely complex profiles. Most of the bursts are 
time asymmetric but some are symmetric. Burst timescales go through the 30 ms scale 
to hundreds of seconds. The measurement of these timescales is a rather complicated 
task, since it may depend on the intensity of both the background and the source. 
At high energies (> 100 MeV), some extremely long bursts have been detected. The 
angular distribution of these bursts is isotropic within the statistical limits, and the 
paucity of comparatively faint bursts implies that we are seeing to near the edge of the 
source population ||203|| . Both effects, isotropy and non- homogeneity in the distribu- 
tion, strongly suggest a cosmological origin of the phenomenon. Moreover, the recent 
detection of "afterglows" , delayed low energy (X-ray to radio) emission of GRBs, has 
confirmed the cosmological origin of the burst via a redshift determination of several 



GRB host-galaxies [204, 205 



If the sources are so distant, the energy necessary to produce the observed events by 
an intrinsic mechanism is astonishing: about 10 51 erg of gamma rays must be released 
in less than 1 second. The most popular interpretation of the GRB-phenomenology 
is that the observable effects are due to the dissipation of the kinetic energy of 
a relativistic expanding plasma wind, a "fireball", whose primal cause is not yet 
known [|06], |07], |08|, ^09], gK| gll], pT^| . The very short timescale observed in the 



time profiles indicate an extreme compactness that implies a source initially opaque 
(because of 77 pair creation) to 7-rays. The radiation pressure on the optically thick 
source drives relativistic expansion, converting internal energy into the kinetic energy 
of the inflating shell. Baryonic pollution in this expanding flow can trap the radia- 
tion until most of the initial energy has gone into bulk motion with Lorentz factors 



of T > 10 — 10 ||213|| . The kinetic energy, however, can be partially converted into 



heat when the shell collides with the interstellar medium or when shocks within the 
expanding source collide with one another. The randomized energy can then be ra- 
diated by synchrotron radiation and inverse Compton scattering yielding non-thermal 
bursts with timescales of seconds. 

We now consider Fermi acceleration in the fireball internal shocks. As usual, r L 
should be smaller than the largest scale Iqrb over which the magnetic field fluctuates, 
since otherwise Fermi acceleration may not be efficient. One may estimate Iqkb as 
follows. The comoving time, i.e., the time measured in the fireball rest frame, is 
t = r/Tc. Hence, the plasma wind properties fluctuate over comoving scale length 
up to Iqrb ~ r /r, because regions separated by a comoving distance larger than 
r/r are causally disconnected. Moreover, the internal energy is decreasing because 
of the expansion and thus it is available for proton acceleration (as well as for 7- 
ray production) only over a comoving time t. The typical acceleration time scale is 
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then |7§ 

GRB _ 

acc ' V > 

where j3c is the Alfven velocity. In the GRB scenario (3 ~ 1, so Eq. sets a 
lower limit on the required comoving magnetic field strength, and the Larmor radius 
r L = E' jeB = E/YeB, where E' — E/Y is the proton energy measured in the fireball 
frame. The dominant energy loss process in this case is synchrotron cooling. Therefore, 
the condition that the synchrotron loss time of Eq. fl38|) be smaller than the acceleration 
time sets the upper limit on the magnetic field strength. A dissipative ultra-relativistic 
wind, with luminosity and variability time implied by GRB observations, satisfies the 
constraints necessary to accelerate protons to energy > 10 20 eV, provided that Y > 100, 
and the magnetic field is close to equipartition with electrons. We stress that the latter 



must be satisfied to account for both 7-ray emission and afterglow observations 11213 ]. 
At this stage, it is worthwhile to point out that for the acceleration process at shocks 
with large Y the particle distributions are extremely anisotropic in shock, with the 
particle angular distribution opening angles ~ T _1 in the upstream plasma rest frame. 
Therefore, when transmitted downstream the shock particles have a limitted chance to 



be scattered efficiently to re-eneter the shock ||214|| . However, in this particular case, 
the energy gain by any "successful" CR can be comparable to its original energy, i.e., 
(AE)/E ~ 1. 

Now we comment on the data. It was pointed out that two of the highest energy 
CRs come from directions that are within the error boxes of two remarkable GRBs 



detected by BATSE with a delay of O(10) months after the bursts |215| |. However, 



a rigorous analysis shows no correlation between the arrival direction of ultrahigh 



energy CRs and GRBs from the third BATSE catalog ] p!6|] . Moreover, no correlations 
were found between a pre-CGRO burst catalog and the Haverah Park shower set that 
covered approximately the same period of time. These analyses, however, could have 
been distorted by the angular resolution (A6 ~ 3°) of the GRB measurements. A 
sensitive anisotropy analysis between ultrahigh energy CRs and GRBs will be possible 
in the near future with the facilities of PAO and the High Energy Transient Explorer 
(HETE). 17 

Finally, it is also interesting to note that if the GRBs are uniformly distributed 
(independent of redshift), and if in the past they emitted the same amount of energy 
in ultrahigh energy (~ 10 14 MeV) CRs as in ~ MeV photons, the energy input of these 
particles into the extragalactic space would be enough to account for the observed CR 



17 This satellite, designed to localize GRBs, is expected to detect around 50 events per year with a 
10 arc-minute accuracy in the medium energy X-ray band, and around 10 events per year with 10 



arc-second accuracy in the soft X-ray band, tittp : / / space . mit . edu/HETE/Welcome . html 
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flux ||176|| . However, recent afterglow studies indicate that their redshift distribution 
likely follows the average star formation rate of the Universe and that GRBs are more 
numerous at high redshift. If this is the case, ultrahigh energy CRs coming from GRBs 
would produce too sharp a spectral energy cutoff to be consistent with the AGASA 
data, because of energy degradation en route to Earth ||217| [21S|| . 18 



IV. INTERACTIONS EN ROUTE TO EARTH 



In this section we review the propagation of cosmic rays. We start with the hadronic 
component, continue with a discussion of electromagnetic cascades initiated by ultra- 
high energy photons in extragalactic space, and then discuss how propagation can 
be influenced by cosmic magnetic fields. For alternative discussions see, for instance, 
Refs. jT70l, ggg. 



The GZK cutoff 



Shortly after the discovery of the microwave echo of the big bang [|222|| , Greisen ||223 
and independently Zatsepin and Kuzmin [|224j1 (GZK) pointed out that this radia- 



tion would make the Universe opaque to protons of sufficiently high energy. At 
energies above a few 10 19 eV, the thermal photons are seen highly blue-shifted 
by the protons in their rest frames. The energy of the relic photons is sufficient 
to excite baryon resonances thus draining the energy of the proton via pion pro- 
duction and, coincident ally, producing a source of ultrahigh energy gamma rays 
and neutrinos. Since the implications of this process have been revisited in many 

forms |g2§ psq pan psq, p2q, pgq pij psq ggg , ggg , ggg , pgq , ggg , ggg , ggg , ggg , gg 



the concept of GZK sphere is by now a somewhat fuzzy notion. Naively, it is the radius 
of the sphere within which a source has to lie in order to provide us with protons of 
10 20 eV. In what follows we extensively discuss the phenomenology of proton-photon 
interactions and discuss in some detail the GZK radius. 

There are three sources of energy loss of ultrahigh energy protons. These are the 
adiabatic fractional energy loss due to the expansion of the Universe, pair production 
(P7 — > V e+e ~) an d photopion production (jry — > ttN), each successively dominating 
as the proton energy increases. The adiabatic fractional energy loss at the present 



18 It has been noted by Bahcall and Waxman [219] that if one excludes the AGASA results, the 
remaining data from Fly's Eye, HiRes, and Yakutsk CR experiments show evidence for a cutoff in 
the CR energy spectrum below 10 20 eV with la significance. In such a case, the door is open for a 



common origin between GRBs and ultrahigh energy CRs [ 220 1 
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cosmological epoch is given by 



1 fdE\ . , 

— ) =#o, (46) 

adiabatic 



E\dt J. 



where H ~ 100 h km s" 1 Mpc" 1 is the Hubble constant, with h ~ (0.71 ± 0.07) x 



1.15 
0.95 



the normalized Hubble expansion rate [|242|| (see Appendix B). The fractional energy 



loss due to interactions with the cosmic background radiation at a redshift z = is 
determined by the integral of the nucleon energy loss per collision multiplied by the 
probability per unit time for a nucleon collision in an isotropic gas of photons ||225 |. 
This integral can be explicitly written as follows, 

IdE c ^ p» , ^ . . f" , n(w) 
- -— = — ^ ^ dw r Kj a 3 (w r ) w r dw — (47) 

where w r is the photon energy in the rest frame of the nucleon, and Kj is the inelas- 
ticity, i.e., the average fraction of the energy lost by the photon to the nucleon in the 
laboratory frame for the jth reaction channel. 19 The sum is carried out over all chan- 
nels, n(w)dw stands for the number density of photons with energy between w and dw, 
<jj(w r ) is the total cross section of the jth interaction channel, T is the usual Lorentz 
factor of the nucleon, and w m is the maximum energy of the photon in the photon gas. 

Pair production and photopion production processes are only of importance for in- 
teractions with the 2.7 K blackbody background radiation. Collisions with optical and 
infrared photons give a negligible contribution. Therefore, for interactions with a black- 
body field of temperature T, the photon density is that of a Planck spectrum [ 243 , so 
the fractional energy loss is given by 



1 dE ckT 



POO 

/ dw r (Jj{w r ) Kj w r ln(l - e ~ Wr/2TkT ) (4c- 



E dt 2n 2 T 2 (ch) 3 

3 j 

where k and K are Boltzmann's and Planck's constants respectively, and Wq. is the 
threshold energy for the jth reaction in the rest frame of the nucleon. 

Berezinsky and Grigor'eva have examined in detail the pair production process [227|. 



At energies E m e m p /kT = 2.1 x 10 eV (i.e., w r /m e - 2 < 1), when the reaction 
takes place on the photons from the high energy tail of the Planck distribution, the 
fraction of energy lost in one collision and the cross section can be approximated by 
the threshold values 

K e+e - = 2 ^ , (49) 
m„ 



19 Here the laboratory frame is the one in which the cosmic microwave background is at a temperature 
w 2.7 K. 
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and 

- 3 

7T 2 / w r 



a e+e -(w r ) = —ar f — ^ - 2 J , (50) 

where a is the fine structure constant and ro is the classical radius of the electron. The 
fractional energy loss due to pair production for E < 10 18 eV is then, 

1 fdE\ 16c m e 2 fkT\ 3 [TkT\ 2 / m e \ 

~ ~E \~dt ) e+e _ = ~~n~ ~m~ v a r ° \~hc ) U7j eXP {-fkf)- (51) 

At higher energies (E > 10 19 eV) the characteristic time for the energy loss due to pair 
production is t ~ 5 x 10 9 yr ||244j| . In this energy regime, the photopion reactions p7 — > 
pir° and p7 — > n + n on the tail of the Planck distribution give the main contribution 
to proton energy loss. The cross sections of these reactions are well known and the 
kinematics is simple. 

Photopion production turns on at a photon energy in the proton rest frame of 
145 MeV with a strongly increasing cross section at the A(1232) resonance, which de- 
cays into the one pion channels n + n and 7T°p. With increasing energy, heavier baryon 
resonances occur and the proton might reappear only after successive decays of reso- 
nances. The most important channel of this kind is p~y — > A ++ 7r~ with intermediate 
A ++ states leading finally to A ++ — > p-K + . A ++ examples in this category are the 
A(1620) and A(1700) resonances. The cross section in this region can be described 
by either a sum or a product of Breit- Wigner distributions over the main resonances 



produced in A^7 collisions considering ttN, ttitN and KA (A — > Nir) final states ||245 



At high energies, 3.0 GeV < w r < 183 GeV, the CERN-HERA and COMPAS Groups 



have made a fit to the cross section [ 246 1 . The parameterization is 



a.K) = A + B In 2 (^) + C In (^) mb , (52) 

where A = 0.147 ± 0.001, B = 0.0022 ± 0.0001, and C = -0.0170 ± 0.0007. In this 
energy range, the cr total (n7) is to a good approximation identical to cr total (p7). 

We turn now to the kinematics of photon-nucleon interactions. Assuming that reac- 
tions mediated by baryon resonances have spherically symmetric decay angular distri- 
butions, the average energy loss of the nucleon after n resonant decays is given by 



K*(mRo) = 1 " ^TT 1 + Z2 > (53) 




where denotes the mass of the i th resonant system of the decay chain, uim the mass 
of the associated meson, = is the total energy of the reaction in the cm., and 
the mass of the nucleon. See Appendix C for details. For multi-pion production 
the case is much more complicated because of the non-trivial final state kinematics. 



45 




FIG. 11: Energy attenuation length of nucleons in the intergalactic medium. Note that after 
a distance of ~ 100 Mpc, or propagation time ~ 3 x 10 8 yr, the mean energy is essentially 
independent of the initial energy of the protons, with a critical energy around 10 20 eV. 



Published in Ref. [234]. 



However, it is well established experimentally ||247 that, at very high energies (\/s ^ 3 
GeV), the incoming particles lose only one-half their energy via pion photoproduction 
independently of the number of pions produced. This is the "leading particle effect" . 

For y/s < 2 GeV, the best maximum likelihood fit to Eq. fl48[) with the exponential 
behavior 

-^(§)r AeM - B/El (54) 
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derived from the values of cross section and fractional energy loss at threshold, 
gives ||234 



A = (3.66 ± 0.08) x 10 yr, B = (2.87 ± 0.03) x 10 20 eV . (55) 

The fractional energy loss due to production of multipion final states at higher cm. 
energies (-y/i > 3 GeV) is roughly a constant, 

- - ( I = C = ( 2 - 42 ± °- 03 ) x 10_8 • ( 56 ) 

From the values determined for the fractional energy loss, it is straightforward to 
compute the energy degradation of ultrahigh energy cosmic rays in terms of their flight 
time. This is given by, 

At - Ei (B/E) + Ei (B/E ) = 0, for 10 19 eV < E < 10 21 eV , (57) 

and 

E(t) = E exp[-Ct], for E > 10 21 eV , (58) 



where Ei is the exponential integral ||248|| . Figure |ll| shows the proton energy degrada- 
tion as a function of the mean flight distance. Notice that, independent of the initial 
energy of the nucleon, the mean energy values approach 10 20 eV after a distance of ~ 
100 Mpc. 

Nevertheless, the definition of the GZK radius is subject to possible magnetic field 
deflections. To obtain quantitative estimates, one generally defines the 50% horizon 
-R50 as the light propagation distance from the source at which 1/e of all injected 
protons have retained 50% or more of their energy, i.e., R 50 is achieved when 



/ — dE = N exp(-l), (59) 

JEo/2 UE 



where Nq is the number of particles injected with energy E$ ||237|| . For injection energies 
> 10 20 eV the protons are not affected much by magnetic field strengths in the range 
1 — 10 nG (more on this below) since their scattering angles are small, and so the horizon 



50 



energy dependence is similar to that of the energy loss distance. At E = 10 2U eV R 
is about 100 Mpc, whereas at 2 x 10 20 eV it decreases to 20 Mpc and becomes smaller 
than 10 Mpc for energies above 3 x 10 20 eV |237| . Below 10 20 eV the scattering in the 



magnetic field increases the propagation time and thus causes additional energy loss 
and an increase of the ratio between the mean energy loss length and -R50. In summary, 
the GZK-cutoff implies that if the primary ultrahigh energy CRs are protons energetic 
sources should be very close to the Earth, say within about 20 Mpc. 
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On a different track, the GZK-cutoff guarantees a cosmogenic ffux of neutrinos, 



caused by the decay of charged pions produced in the photon-nucleon interactions ||249 



250 



The resulting neutrino flux depends critically on the cosmological evolution of the 
cosmic ray sources and on their proton injection spectra ||230| , [232j , |252|| . Of course, the 
neutrino intensity will also depend on the homogeneity of sources. For example, semi- 
local objects, such as the Virgo cluster [|226j| , could countribute to the high energy tail 
of the neutrino spectrum. Additionally, there is a weak dependence on the details of the 
cosmological expansion of the Universe. For example, a small cosmological constant 



tends to increase the contribution to neutrino fluxes from higher redshifts [[252 

In Fig. [12| we show the cosmogenic neutrino flux obtained by Protheroe and John- 
son ||232| using Monte Carlo simulation techniques for a proton injection spectrum with 
-^cutoff = 3 x 10 21 eV. This study incorporates the cosmological evolution of the sources 
from estimates [|161|| of the power per comoving volume emitted as protons by powerful 
radiogalaxies, taking into account the radio luminosity functions given in Ref. ||254|| . 
The flux peaks around E « 2 x 10 17 eV, which is roughly the same energy suggested by 
other analyses following a source evolution scaling like (1 + z) 4 |230| , |252|| . The figure 
also shows the flux-estimates by Hill and Schramm |[226|| considering contributions from 
semi-local nucleon sources, and an earlier estimate by Stecker which does not take into 
account the source evolution ||250 |. 

In summary, cosmic microwave background (CMB) renders the Universe opaque to 
protons above about 5 x 10 19 eV. If proton sources are are at cosmological distances 
(> 100 Mpc), the observed proton spectrum should display the GZK cutoff. 



B. Photonuclear interactions 

The relevant mechanisms for the energy loss that extremely high energy nuclei suffer 
during their trip to Earth are: Compton interactions, pair production in the field of the 
nucleus, photodisintegration, and hadron photoproduction. The Compton interactions 
have no threshold energy. In the nucleus rest-frame, pair production has a threshold 
at ~ 1 MeV, photodisintegration is particularly important at the peak of the giant 
dipole resonance (15 to 25 MeV), and photomeson production has a threshold energy 
of ~ 145 MeV. 

Compton interactions result in only a negligibly small energy loss for the nucleus 



given by ||255 



dE Z 4 ( E X _ x 

~1I = T^\a^) eVs (60) 

where p 7 is the energy density of the ambient photon field in eV cm -3 , E is the total 
energy of the nucleus in eV, and Z and A are the atomic number and weight of the 
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FIG. 12: Cosmogenic + + z^ e fluxes from Protheroe and Johnson with energy cutoff of 
3 x 10 21 eV (solid) |251|| , Hill and Schramm (dashed) fl226|| , and previous estimate by Stecker 
without source evolution (dotted) |25C| ]. Published in Ref. [253]. 



nucleus. The energy loss rate due to photopair production is Z 2 /A times higher than 
for a proton of the same Lorentz factor | |256|| , whereas the energy loss rate due to 
photomeson production remains roughly the same. The latter is true because the 
cross section for photomeson production by nuclei is proportional to the mass number 



A P57 |, while the inelasticity is proportional to I /A. However, it is photodisintegration 
rather than photopair and photomeson production that determines the energetics of 
ultrahigh energy cosmic ray nuclei. During this process some fragments of the nuclei 
are released, mostly single neutrons and protons. Experimental data of photonuclear 
interactions are consistent with a two-step process: photoabsorption by the nucleus to 
form a compound state, followed by a statistical decay process involving the emission 
of one or more nucleons. 

Following the conventions of Eq. (f47|), the disintegration rate with production of i 
nucleons is given by ||258 | 



R 



Ai 



2r 2 



dw 



n[w) 
w 2 



2Tw 



dlV r W r <JAi(w r 



(61) 
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with OAi the cross section for the interaction. Using the expressions for the cross section 
fitted by Puget, Stecker and Bredekamp ||255|| , it is possible to work out an analytical 
solution for the nuclear disintegration rates (see Appendix D). Summing over all the 
possible channels for a given number of nucleons, one obtains the effective nucleon loss 
rate R = ^2 t iRAi- The effective nucleon loss rate for light elements, as well as for 



those in the carbon, silicon and iron groups can be scaled as [ 255 1 



dA 
~~dl 



dA 
~dt 



Fc 



56/ ' Fc V56 



with the photodisintegration rate parametrized by [|259| 

R(T)\ Fc = 3.25 x l(r 6 r- a643 exp(-2.15 x 10 10 /r) s" 1 
for T e [1.0 x 10 9 ,36.8 x 10 9 ], and 



mi 



1.59 x 10 



-12 -p-0.0698 -1 



(62) 



(63) 



(64) 



for T G [3.68 x 10 10 , 10.0 x 10 10 ]. 

For photodisintegration, the averaged fractional energy loss results equal the frac- 
tional loss in mass number of the nucleus, because the nucleon emission is isotropic 
in the rest frame of the nucleus. During the photodisintegration process the Lorentz 
factor of the nucleus is conserved, unlike the cases of pair production and photomeson 
production processes which involve the creation of new particles that carry off energy. 
The total fractional energy loss is then 



1 dE 
E~dT 



1 dT R 
f~dl + A' 



(65) 



For w r ^ 145 MeV the reduction in T comes from the nuclear energy loss due to pair 
production. The 7-ray momentum absorbed by the nucleus during the formation of the 
excited compound nuclear state that precedes nucleon emission is e>(10" 2 ) times the 
energy loss by nucleon emission ||260|1 . For F > 10 10 the energy loss due to photopair 
production is negligible, and thus 



E(t) ~ 938 A(t) T MeV ~ E exp 



56 



(66) 



Figure |TJ shows the energy of the heaviest surviving nuclear fragment as a function 
of the propagation time, for initial iron nuclei. The solid curves are obtained using 
Eq. fl66j)) whereas the dashed and dotted-dashed curves are obtained by means of Monte 
Carlo simulations [ 262| |. One can see that nuclei with Lorentz factors above 10 10 cannot 
survive for more than 10 Mpc. For these distances, the approximation given in Eq. 
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FIG. 13: The energy of the surviving fragment (Tq = 4 x 10 9 , To = 2 x 10 10 ) vs. prop- 
agation time obtained using Eq. (|66|) is indicated with a solid line. Also included is the 
energy attenuation length obtained from Monte Carlo simulations with (dashed) and with- 
out (dotted-dashed) pair creation production, for comparison. The region between the two 
dotted lines includes 95% of the simulations. This gives a clear idea of the range of values 
which can result from fluctuations from the average behaviour. It is important to keep in 
mind that a light propagation distance of 1.03 x 10 s corresponds to 1 Mpc. From Ref. [ p61| . 



always lies in the region which includes 95% of the Monte Carlo simulations. When the 
nucleus is emitted with a Lorentz factor T < 5 x 10 9 , pair production losses start to 
be relevant, significantly reducing the value of T as the nucleus propagates distances of 
(9(100 Mpc). The effect has a maximum for T ~ 4 x 10 9 but becomes small again for 
T < 10 9 , for which appreciable effects only appear for cosmological distances (> 1000 
Mpc), see for instance ||262| . 

Note that Eq. (|66|) imposes a strong constraint on the location of nucleus-sources: 
less than 1% of iron nuclei (or any surviving fragment of their spallations) can survive 
more than 3 x 10 14 s with an energy > 10 20,5 eV. For straight line propagation, this 
represents a maximum distance of ~ 3 Mpc. 

Because of the position of 56 Fe on the binding energy curve, it is generally considered 
to be a significant end product of stellar evolution, and indeed higher mass nuclei are 
found to be much rarer in the cosmic radiation. Specifically, the atomic abundance of 



51 



nuclei heavier than iron (in the cosmic radiation) is expected to be smaller by 3 or 5 
orders of magnitude relative to the lighter ones ||263|| . However, for super-heavy nuclei 
the GZK sphere is substantially increased, because for T < 10 9 the nucleus interacts 
only with the tail of the Planck spectrum and infrared backgrounds ||264| , |265|| . Indeed, 
because of the paucity of data, there is a great uncertainty in the chemical composition 
at the end of the spectrum, so the situation is ripe for speculation on super-heavy nuclei 
as primaries p66| . 



C. Relativistic specks of dust 

Dust is a very widespread component of diffuse matter in the galaxy and, apparently, 
in intergalactic space. In particular, specks expelled by radiation pressure from cool 



stars may be injected into the interstellar medium with speeds ~ 10 cm s ||267 
generally carrying a net electric charge. Therefore, these small solid particles can be 
re-accelerated very effectively to ultra-high energies at shock-wave fronts |[268| , |269| , |270 



However, the situation is unusual, because the Lorentz factors involved can be small 
even at very high energies. 

We now discuss the survival probability of such specks of dust. The break up of dust 
grains is mainly due to electrical stresses induced by the accumulation of charge as a 
result of the photoelectric effect in the galactic optical radiation and also as a result 
of ionization by the hydrogen atoms of interstellar gas. The critical electric charge on 
the surface of a dust grain of radius R, for which repulsion forces break up the particle 
is given by | |271| ] 

q c =^hR\ (67) 

where fo denotes the breaking stress. Subrelativistic (ti^3x 10 9 cm/s) dust grains 
disintegrate in collisions with the nuclei and electrons of the interstellar medium. Ion- 
ization takes place along the path of the nucleus in the dust grain, heating the matter 
up to 10 5 K in a narrow channel of radius ~ 10 -7 cm. Close to the surface of the 
channel the local pressure reaches 10 4 kg mm -2 , which results in microfractures. Be- 
cause of the small size of the grain, microfractures are unstable and increase until the 
complete break-up of the particle occurs. The charge accumulation rate on a small 
grain (R < R) is given by ||272| 



^ = 0.6 A- 1 v n H vr R 2 (£) * (pR) 0A1 (61 
and on a large one (R> R) 



'''' 2 A~ x v n H tt i? 2 a T N A p (£) * -^£= , (69) 
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with 

R=l.lxl0 2 p 1A4 f 1 - 22 t im. (70) 

Here, v is the velocity of the dust grain, p its density, riu is the density of hydrogen 
atoms in the surrounding medium, A is the mass number, Na is Avogadro's number, 
and <jt is the Thomson cross section. Significant accumulation of charge occurs during 
the trip through the galaxy, where the density of hydrogen is ~ 0.6 atoms/cm 3 . 
The path length up to the first break up in the galactic gas, 



turns out to be D < 10 15 cm, for graphite grains (fo ~ 1 kg mm -2 ), and D < 10 16 cm, 
for iron grains (fo ~ 10 kg mm -2 ) [|272|| . 

The threshold value of the Lorentz factor of a relativistic grain with energy E and 
density p (measured in g cm -3 ), under the influence of the photoelectric effect with 
photons of energy w, is given by [[271| 



, , — = 9.0 x 10 -2 J0 , GcV , (72) 
AnpJ \wcj y w%p J 

where E GcV = E/GeV, w cV = w/eV, and f is measured in kg mm -2 . Thus, a carbon 
dust grain with E ~ 10 20 eV and T > 25 is split up by starlight photons (w ^ 2 eV). 
Relativistic dust grain break up is particularly effective within the solar system. The 
distance the grain can approach the Sun without splitting up is given by [ 271] 



D = 2.3 x 10 -25 A 1_ ; cm, 73 

where L(w) is the number of photons with energy above threshold emitted by the Sun 
measured in s -1 , and x ~ 0.01 — 0.1 is the number of electrons stripped from the 
speck of dust by one photon. For carbon dust grains with r > 25 (uj ^ 2 eV and 
L(w) ~ 3.5 x 10 44 s -1 ), Eq. (^|) leads to D ~ 10 pc. Figure |l| shows the path length 
of a typical graphite grain in the field of the galactic optical radiation. 

Doubts have also been expressed about the prospects of surviving against the heating 
from photoionization within the solar radiation field [ [273|| . Consider a grain moving 



(radially) towards the Sun with constant velocity v = (3c and subtending an effective 
solid angle dQ, with a the angle between the position vector f of the grain and the 
positive x-axis defined by a line parallel to the direction of motion. In the rest system 
of the grain, the average energy of the photons is blue-shifted to Tw(l — (3 cos a). The 
Sun no longer appears to emit isotropically; namely, the number of incident photons 
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FIG. 14: The path length of a graphite grain of radius R = 1 /im up to the first break-up, in 
the field of the galactic optical radiation (GOR). Adapted from Ref. | |272| |. 

per second is T L{w) (1 — (3 cos a)cfO/47r. To estimate the surface temperature T of the 
grain we assume that it radiates black body, 



e a B T 4nR z = T L(w) w (1 - (5 cos a) 



An 



(74) 



where <jb is the Stefan-Boltzman constant and e the emissivity of the grain. The solid 
angle is determined by the cross section for absorption by the grain. For a particle 
traveling in the r-direction towards the Sun, 



dtt = 4 N a 7c R 3 r 



(75) 



where iV is the volume density of grain atoms and a the ionization cross section per 
atom. For iron, o ~ 3 x 10 _16 r _1 cm 2 [ 274]. Now, setting T equal to the melting 



point T c of the grain, Eqs. ( |74D and (75) combine to yield the following expression for 
the position r a at which the grain melts 



where 



[1 — (3 cos a) 



! T 2 RL(w)wNa 
127rea B T r 4 



(76) 
(77) 
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For an iron grain with T c ~ 1800 K, e ~ 0.5 and N ~ 8.6 x 10 22 atoms cm 3 , Eq. (j77|) 
leads to 

f a ~ 3 x 10 15 VTR cm. (78) 

Now, since 

E = -7iR 3 c i Tp. (79) 
3 

for an iron grain of E ~ 10 19 eV, f a = 1.0 x 10 13 T 1 / 3 cm. Hence, relativistic (logT < 2) 
iron grains would melt before reaching the Earth's orbit around the Sun, i.e., r a ~ 1 
au = 1.5 x 10 13 cm. There is no substantial change for an iron grain traveling non- 



radially, it would melt before reaching the Earth ||273|| . However, note that dust grains 
composed of carbon (T c ~ 4000 K) may reach the Earth's orbit with E ~ 10 20 eV if 
r ~ 20. 

In summary, subrelativistic dust grains have no chance to survive the trip to Earth. 
However, relativistic carbon specks could survive a trip of ~ 100 pc if one assumes 
the most optimistic parameters. Since the closest supernova remnant is about 200 pc 
away ||275|| , astrophysical models cannot rule out relativistic dust grains as ultrahigh 



energy CR primaries. 



D. Electromagnetic cascades 



The dominant absorption process for high energy 7-rays (energy E) traversing cosmic 
distances is pair creation through collisions with the various radiation fields (energy 
w) permeating the Universe pTl p77], E7[l 1279], EBO, [281], E82f. The fundamental 



process is well understood and its amplitude can be calculated accurately by general 
perturbation methods developed in quantum electrodynamics (QED) [ f283| , [284| , |285 |. 
The lowest-order single pair creation cross section ||276| 



a 



(3 - I3 A ) In 



1 + 
1-P 



2/3(2 - (3 2 ) 



50) 



peaks near the threshold E t h = m 2 e /w ~ 2.6 x 10 2 (w/eV) 1 GeV, and falls off asymp- 



totically as s — > 00 ||285 



4ira z 



a[s) ~ 



In I — 1 - 1 



Here, (3 = (1 — 1/s) 1 / 2 is the electron (positron) velocity in the cm. system and tq 
is the classical electron radius. Higher order QED processes with more than two final 
state particles become important with rising energy. For example, the fourth-order 



55 



two-pair creation cross section is a sharply rising function of s near the threshold that 
quickly approaches its asymptotic value [[285|| 20 



~ t 175 C( 3 ) " 38 ] ~ 6 - 45 ^ , ( 82 ) 

and becomes bigger than the single pair creation cross section when s <; 1 GeV 2 . With 
this in mind, the most efficient targets for 7-rays of energy E are background photons 
of energy w ~ m 2 e / E. This implies that for 7-ray energies > 10 20 eV, interactions with 
the radio background (w < 10" 6 eV ~ 100 MHz) become more important than those 
with the CMB. 

We have no direct knowledge of the cosmic radio background, mostly because it is 
difficult to disentangle the galactic and extragalactic component. For a temperature 
of the extragalactic component of 15 K at 178 MHz one finds 

n(w) ~ Kw~ 2 , (83) 

with K = 1.09 x 10 20 erg/cm 3 |[286|| . At frequencies somewhere below 1 MHz the in- 
tensity is expected to fall off exponentially; the principal absorption process would be 
inverse Bremsstrahlung or "free-free" absorption by the intergalactic plasma. Obser- 



vational estimates are consistent with a cutoff in the spectrum at w ^ 10~ 8 eV [|287 
However, that part of the spectrum where the Universe is optically thick to radio emis- 
sion depends on the abundance and clustering of electrons in the extragalactic medium 
and/or the radio source, and is uncertain between 0.1-2 MHz. A theoretical esti- 
mate ||288|| of the intensity down to kHz frequencies, based on the observed luminosity 



functions and radio spectra of normal galaxies and radio-galaxies, tends to give higher 
estimates than those of Ref. 



The probability per unit path length that a 7-ray is converted into e + e pairs in a 



collision with a low-energy photon is given by [ 276 ] 



— = 2 — f / ^-^-dw I sa(s)ds, (84) 



dx \ E J J w 

where sq = Ew/m e . The mean interaction length of pair production, showing the con- 
tribution of the various radiation fields, is given in Fig. For the radio background, 
we assume the conservative estimate of Ref. [|287|| . Using the radio spectrum derived by 



Protheroe and Biermann leads to a photon mean free path (at E ~ 10 eV) a factor 



of 3 - 10 smaller 288 



In the Klein-Nishina limit, s 3> m e) one of the pair-produced particles carries most of 
the 7-ray energy. In the absence of magnetic fields, this leading electron (positron) can 
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£(3) ~ 1.202 is the zeta function with argument 3. 
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FIG. 15: Mean pair production interaction length for 7-rays in the infrared and optical 
backgrounds (IR), radio background (R), and relic photons (CMB). Adapted from Ref. [{232(| , 

deliver a substantial fraction of the original 7-ray's energy to a new 7-ray via inverse 
Compton scattering (ICS), with cross section 

(85) 

The scattered photon, which is now the leading particle, can initiate a fresh cycle of pair 
production and inverse Compton scattering interactions, yielding an electromagnetic 
cascade. 
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For s <C m 2 , dies ~ ct, and the fractional energy loss is given by ||289 



1 / dE\ 4 „ wn(w) , . 

-— — ~-(T T cr 86 

E \ at J ICS 3 m e 

Corrections to the lowest order ICS cross section from processes involving additional 
photons in the final state are less than 10% in the ultrahigh energy range [290|| . Fur- 



thermore, various other possible processes, such as those involving production of one 
or more muons, taus, pion pairs, and Bethe-Heitler pairs (jX — > Xe + e~, where X 
stands for an atom, ion or free electron), can be safely neglected [[2851 . For example, 
the total cross section for the lowest-order single muon pair production is around an 
order of magnitude smaller than that for electron pair production. 

The development of electromagnetic cascades depends sensitively on the strength of 
the extragalactic magnetic field B, which is rather uncertain. For B > 10~ n G the 
electrons will lose much of their energy through synchrotron radiation before they have 
the opportunity to undergo ICS, and the electromagnetic cascade will be virtually 
terminated |p91 ]. Specifically, an electron of energy E moving in a field having a 



perpendicular component of magnitude B± radiates energy cit £l rcite 



dE _ _ ___ u /B ± Y (Ei _, 



_ = 3.3X10-"^ ^— 2 ) eVcm-. (87) 

Hence, if the radiated energy is less than say 10% of the electron energy in an interaction 
length A, the electron synchrotron cools before it can undergo ICS, and thus the cascade 
development stops. The above cooling time scale holds provided that 

B\ < 7.9 x 10 23 1 G 2 , (88) 

E e V Km 

where E aV = E/leV and A cm = A/1 cm. Then, for 3 x 10~ 12 G, 2 x 10- n G, and 
1 x 10 -10 G, the intergalactic space becomes opaque to the propagation of 7-rays of 
E = 10 21 eV, E = 10 20 eV, and E = 10 19 eV, respectively. 

Putting all this together, the GZK radius of the photon channel strongly depends on 
the strength of extragalactic magnetic fields. In principle, distant sources with z > 0.03 
can contribute to the observed rays above 5 x 10 19 eV if the extragalactic magnetic 
field does not exceed 10 -12 G ||292| . However, comparison of the expected and observed 
numbers of CRs from the direction of the nearby source Centaurus A allows one to 
place bounds on intergalactic magnetic fields of 3> 10~ n G [|197[|. In such a case, the 



survival probability for 3 x 10 20 eV 7-rays to distance D, 

P(> D) = exp[-L>/6.6 Mpc] , (89) 
becomes less than 10~ 4 after traversing a distance of 50 Mpc [|172|| . 



58 



E. Deflection and delay due to magnetic fields 



Cosmic rays en route to Earth suffer deflection and delay in magnetic fields, effects 
which can camouflage their origins. In this section, we discuss bending of cosmic ray 
orbits in both Galactic and extragalactic fields. Various combinations of field strength 
and coherence length are considered, and we discuss solutions to the diffusion equations 
as well as results from Monte Carlo simulation. 

The magnetic deflection of protons in the Galactic disk has been studied in detail 



by Stanev [293|. This analysis includes two extreme options for the behaviour of the 
field, reflecting the different symmetries with respect to field reversals in the r- and 
^-directions. The B-field has a 1/r behaviour, with deviations calculated out to 20 kpc 
from the Galactic center. The r.m.s. deviation averaged over arrival direction, for an 
energy of 2 x 10 19 eV (4 x 10 19 eV), varies from 17.7° to 23.7° (7.9° to 10.5°) in going 
between the two models. This deviation shows an approximately linear decrease with 
increasing energy. In the case of heavy nuclei, the arrival direction of the particles is 
strongly dependent on the coordinates of the source [294], p95| , |296|1 . The Larmor radius 
of a nucleus with E/Z = 10 19 eV in a uniform magnetic field of strength B ~ 3 fiG is 
slightly larger than 3 kpc. The typical large scale field intensity of the Milky Way is a 
few /iG, and is approximately uniform over scales of the order of a few kpc. Therefore, 
the propagation in the Galactic magnetic field of a nucleus with E/Z > 10 19 eV should 
in general be not very different from a quasi-rectilinear trajectory, with deflections away 
from the straight path becoming smaller with increasing energy. However, for 10 18 eV < 
E/Z < 10 19 eV, the deflection may be high enough (> 20°) to make any identification 
of the sources extremely difficult. 21 Moreover, lensing effects yield (de)magnification 
of the CR-fluxes, and not every incoming particle direction is allowed between a given 
source and the detector. This generates sky patches which are virtually unobservable 



from Earth ||294| , |295| . On the other hand, the enhancement of the probability to 
detect events from a given source in a narrow energy range near the caustic implies 
a concentration of events around the location where the new pair of images forms. 
Therefore, magnetic lensing around caustics becomes a potential source of clustering 



in the angular distribution of arrival directions |[296|| . Precise predictions depend upon 
the detailed structure of the magnetic field configuration, which is not so well known. 
However, new experiments such as PAO may make it possible to use CRs as probes of 



the magnetic field structure |[298| , [299|| 



At present, surprisingly little is known about extragalactic magnetic fields (EGMFs). 
There are some measurements of diffuse radio emission from the bridge area between the 
Coma and Abell superclusters ||300| , which, under assumptions of equipartition allows 



21 



At smaller energies the drift and diffusive regimes dominate [ 297 1 
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an estimate of 0(0.2 — 0.6) fiG for the magnetic field in this region. 22 Such a large field 



may possibly be understood if the bridge region lies along a filament or sheet ||302|| . In 
addition, the absence of a positive signal in Faraday rotation measurements on distant 
quasars ||303|| provides upper limits on magnetic fields of any origin as a function of 
reversal scale 



0041 



These bounds depend significantly on assumptions about the 
electron density profile as a function of the redshift. When electron densities follow 
that of the Lyman-a forest, the average magnitude of the magnetic field receives an 
upper limit of B ~ 10 -9 G for reversals on the scale of the horizon, and B ~ 10~ 8 G 

The latter upper bound is roughly 
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for reversal scales on the order of 1 Mpc 
coincident with the lower limit in the Galactic neighborhood [197|. 

If the Larmor radius of a particle is sizably larger than the coherence length i of 
the extragalactic magnetic field, the typical deflection angle from the direction of the 
source, located at a distance d, can be estimated assuming that the particle makes a 
random walk in the magnetic field ||306|| , 

d y /2 ( t ^ 1/2 

1 Mpc J 



9(E) ~ 0.54° Z 



1 Mpc 



B 



nG 



E- 



(90) 



20 



where B n Q = B/l nG, and E20 = E/10 eV. Note that deflection also implies an 
average time delay 

, „ s d $ 2 . 

relative to rectilinear propagation with the speed of light. With increasing magnetic 
field strength, the Larmor radius decreases and the hypothesis of straight line propaga- 
tion between scatterings, that Eq. ( p0[) relies upon, breaks down. A diffusive approach 
is more appropriate for this situation ||307[ [308| , |309| , |310| , |311| , p!2| , |313 |. 

The Local Supercluster is a flattened overdensity of galaxies extending for ~ 30 — 
40 Mpc with a width of ~ 10 Mpc ||272| . The Virgo cluster is approximately in 
the center of this distribution, about 17 Mpc away from the Milky Way, which is 
located at the edge. This large scale distribution defines the super-Galactic plane 
where it intersects almost perpendicular the Galactic plane. A broad ~ 50° region 
centered at the Local Supercluster seems to be endowed with a magnetic field with a 
uniform component of strength ~ 1.5 fiG 1314 ], and a random component of strength 
~ 1 /iG ||315|| . If this is indeed the case, all CRs traveling to Earth will propagate 
diffusively. 

The evolution of the nucleon spectrum is governed by the balance equation 



On 
~dt 



d[b(E)n] 
dE 



D(E) V 2 n + Q 



(92) 



22 Fields of O(fxG) are also indicated in a more extensive study of 16 low redshift clusters 1 301 
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which takes into account the conservation of the total number of particles in the spec- 
trum, i.e., J nd 3 f = Uq. In the first term on the right, b(E) = dE/dt is the mean 
rate at which particles lose energy, and the last term, Q is the number of nucleons per 
time generated by the source(s). The second term stands for the diffusion in the extra- 
galactic medium. Two extreme regimes can be distinguished depending on whether the 
particle remains trapped in magnetic subdomains or not, yielding different functional 
dependence on energy of the diffusion coefficient pi6 



1 B 2 

D(E) = -r,c- F55 — . (93) 

y J 3 L jTdkP(k) V ; 

Here, P(k) is the magnetic field power spectrum. As mentioned earlier, extragalactic 
magnetic field strengths and coherence lengths are not well established, but for r L < £, 



it may be plausible to assume a Kolmogorov [|188|| form for the turbulent magnetic field 
power spectrum, P(k) = Po(k / ko)~ 5 ^ 3 , with coherent directions on scales of 0.5 - 1 Mpc. 
Here, ko ~ l/£is the small wavenumber limit for the magnetic field and 3 Pq k /2 = B 2 . 
One can now easily estimate that protons with energies E c < 10 21 £mpc-B^g eV remain 
trapped inside magnetic subdomains of size £m P c = £/(l Mpc), attaining efficient dif- 
fusion when the wave number of the associated Alfven wave is equal to the Larmor 



radius of the particle [|189| , |316[ . The replacement of the Kolmogorov spectrum into 
Eq. © leads to 

(J£ fi2 \ V 3 
20 Mpc j Mpc 2 /Myr. (94) 

As the particle energy increases (E > E c ), there is a transition to Bohm diffusion 
when r L ~ £, and the diffusion coefficient is of the order of the Larmor radius times 
the velocity. Numerical simulations [ fjQ9|| show that for Bohm diffusion, D ~ 3r L c, or 
equivalents, D{E) ~ 0.1E 2Q /B IJlG Mpc 2 /Myr. 23 

Now we consider the propagation of a burst from a single source. Idealizing the 
emission to be uniform with a rate dN/dt = n /r, we have 

Q = — 5 3 (r) [9(f - to„) - 8(f - toff)], (95) 
r 

where J Qd 3 rdt' = N tot , and 9 is the Heaviside step function. Here t on (t oS ) is the 
time since the engine turned on (off) its CR production. If the energy loss term is 
neglected, the solution to Eq. ( j92f) reads, 

n(f,t)= [ dt' [ d 3 fG(f-r',t-t') Q(7,t f ), (96) 



23 Criticisms of these assumptions have been raised in 131 
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where 



G(r-f*,t-t r ) = [AnD(t-t')}- 3/2 Q(t - t') exp{-(r - f>) 2 /AD(t - t')} , (97) 



is the Green function ||318|| . Before integrating over the range of possible propagation 
times, we note that there is a negligible contribution from times prior to the arrival 
time of the diffusion front, 

r 2 

TD = AD 

For a bursting source, i.e, t Q g — t on ~ 5t, Eq. (R 



{9i 



n{r,t) = — 

T 

n 



dt' 



-r 2 /4,(t-t>) 



leads to [ pi9| 
1 



n 



[4tt D(t -t')} 3 / 2 rAirDr 



'AD(t - t' 



[47TD] 3 / 2 



T 



-1/2 



-t: 



-1/2 



+ 0{r 3/2 )] 



(99) 



where $>(x) is the error function, and we have set present observation time at t — 0. 
Note that for Kolmogorov diffusion D ~ i? 1 / 3 . Now, shock acceleration results in a 
spectrum at the source Q oc E~"< with 7 slightly grater than 2. Therefore, the purely 
diffusive solution results in a spectrum n(E) oc £' _7+1 / 3 , close to the observed spectrum 
E~ 2 - 7 BOa. The condition of diffuse propagation is 



To > r s 



(100) 



where t s = r/c is the time required for a straight line propagation. 

If photopion production becomes important, one has to resort to numerical Monte 
Carlo simulations. Detailed studies of CR propagation in turbulent plasmas have been 
reported ||309|, |311||. In these analyses the magnetic field is characterized by a Gaussian 



random field with zero mean and a power spectrum with (B 2 (k)) oc k nB for k < ko and 
(B 2 (k)) = otherwise, where ko = 2n /£ represents the numerical cutoff scale and the 
r.m.s. strength is B 2 = J Q °° dkk 2 (B 2 (k)). The field is then calculated on a grid in real 
space via Fourier transform. Figure |16| shows an example of the distribution of arrival 
time and energies obtained via simulation of a bursting source located at d = 10 Mpc 
and assuming n B = —11/3, i.e., a turbulent Kolmogorov type spectrum [[309i |311| . 
Note that for B = 0.3 fiG, £ = 1 Mpc, r = 10 Mpc, and E = 10 19 eV, Eq. fl94]) leads 
to to = 750 Myr, in the midregion of the time-delay given in Fig. 16. For E = 10 20 eV 
(where proton energy losses are still small), td = 350 Myr, not inconsistent with the 



results of Fig. [16], although at the margin of the distribution of time-delays for this 
energy. This marginal agreement is indicative of the transition to Bohm diffusion, 
where the validity of Eq. (|94]) begins to flag. 
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FIG. 16: The distribution of time delays t(E) and energies E for a burst with spectral index 
7 = 2.4 at a distance d = 10 Mpc and r.m.s. strength B = 3 x 10 -7 G. The inter-contour 
interval is 0.25 in the logarithm to base 10 of the distribution per logarithmic energy and 



time interval. Published in Ref. [309] 



V. OUR COSMIC BACKYARD 

In order to analyze the effect of energy losses on the observed spectrum, it is conve- 
nient to introduce the accumulation factor / acc , defined as the ratio of energy- weighted 
fluxes for "low" (10 18 ' 7 eV - 10 19 5 eV) and high (> 10 20 eV) energy CRs above the 
ankle. For ordinary baryonic CRs, if the Earth is located in a typical environment and 
all CR-sources have smooth emission spectra, the observed spectrum above the ankle 
should have an offset in normalization between low and high energy given by / acc . In 
the case where the cosmic rays are protons from a uniform distribution of sources active 
over cosmological times, the cutoff due to photopion processes relates the accumulation 



factor to a ratio of GZK distances [319] and leads to / acc ~ 100. A similar value for 



f acc is obtained for nuclei due to photodisintegration. The smoothness of the observed 



CR spectrum P0[ , i.e., / acc ~ 1, suggests that the power of nearby sources must be 
comparable to that of all other sources (redshift z > 0.5) added together. Otherwise, 
one would have to invoke an apparently miraculous matching of spectra to account 
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for the smoothness of the CR energy spectrum. Of course, the GZK cutoff can be 
overcome if nearby sources are significantly more concentrated. However, this does not 
seem to be the case: if one simply assumes that the distribution of CR sources follows 
the distribution of galaxies, the local overdensity is only a factor of two above the 
mean, and thus insufficient to explain the measured flux above 10 20 eV 
proposals with a smooth high energy spectrum due to nearby sources have been put 
forward. In this section we discuss these models. 



32C]. 24 Several 



A. Galactic sources 



Galactic sources would satisfy the power condition trivially. However, the observed 
approximate isotropy would be difficult to reconcile with a galactic explanation, unless 
CRs diffuse in the magnetic field. Specks of dust could constitute diffuse galactic CRs, 
though the depth of shower maxima recorded so far are not consistent with the exact 
picture of showers initiated by dust grains ||322| , fj23| 
in the Galactic magnetic field [ |324| 



Heavy nuclei will also diffuse 



Indeed, it was suggested by Blasi, Epstein, and 
Olinto ||175|| that iron ions from the surfaces of newly formed strongly magnetized 



pulsars may be accelerated to super-GZK energies through relativistic MHD winds. 
The iron ejected with energies ~ 10 20 eV will reach the Earth after some diffusion 
through the Galactic magnetic field. Note that the Larmor radius (Z 2 q = Z/26) 
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J 2f> 



B 



3fiG 



-1 



kpc. 



(101) 



is a few times the typical distance to a young neutron star, which is of order d ~ 8 kpc. 

Pulsars begin their life with very fast rotation, lu ~ 3000 s _1 , and very large magnetic 
fields, B s > 10 13 G. Inside the light cylinder, i.e 
density 

(102) 



r < c/uj, a magnetosphere of 



B(r)uj 

n„Ar) ~ 

GjV ; 4vrZec 



corotates with a dipole magnetic field component that scales as B{r) = B s (r ns /r) 3 . As 
the distance from the star increases, the dipole field structure cannot be maintained 
and beyond the light cylinder the field is mostly azimuthal. From the light cylinder 
a relativistic plasma with Alfven speed close to the speed of light expands as a MHD 
wind. 



24 We note that an earlier analysis 1 32 1 overestimated the local overdensity by an order of magnitude. 
This large discrepancy was caused by neglecting the necessary galaxy selection functions which 



account for the fact that nearby galaxies are easier to detect than far away galaxies J3 20| 
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The surface of a young neutron star is composed of elements near to the iron peak 
formed during the supernova event. Iron ions can be stripped off the hot surface of 
the neutron star by strong electric fields and can be present throughout much of the 
magnetosphere ||326| , [32711 . The typical CR energy can be estimated by considering the 



magnetic energy per ion at the light cylinder 



£cr = , (103) 



where B\ c ~ 10 10 Bi 3 u 3k G is the magnetic field strength at the light cylinder, ncj ~ 
1.7 x 10 11 B 13 cuj k /Z cm" 3 , u 3k = u;/3000 s" 1 and B 13 = B s /10 n G. Replacing by 
fiducial values, Eq. (|103|) leads to [[175(1 



E CR = 4 x 10 2t % 6 B l3 iu 2 3k eV , (104) 

less than the maximum energy of particles that can be contained in the wind near the 
light cylinder PH,|329|. 

The injection spectrum is determined by the evolution of the rotational frequency. As 
the pulsar spin slows down because of electromagnetic and gravitational radiation ||330 



the energy available for CR production decreases. For B s > 10 13 G, the rotation 



speed decreases mainly by magnetic dipole radiation [ 332 



^I ^/S^ > ( 105 ) 

L+ts #13 uj i3k 

where uj^k is the initial spin period and tg = t/10 8 s. With this in mind, the particle 
spectrum from each neutron star is given by [|175|1 

N(E) = 4^|^ GeV- 1 , (106) 

-t>13 -^20 ^26 

where e < 1 is the efficiency for accelerating particles at the light cylinder. 

Even though young neutron stars are surrounded by remnants of the presupernova 
star, the accelerated particles can easily escape the supernova remnant without sig- 
nificant degradation for a wide range of initial magnetic fields and spinning rates. 
Specifically, consider a supernova that imparts -Esn ~ 10 51 -Esi erg on the stellar enve- 
lope of mass M env ~ 10Mi solar masses, where E 5 i = _E/10 51 erg, and Mi is the mass 
of the envelope in units of 10 M Q . The column density of the envelope surrounding 
the neutron star is given by 

E ~ - Mcnv ^ (107) 

where v e ~ y^2E SN /M env is the envelope dispersion velocity, t is in seconds, and r ~ 
10 13 ri3 cm. The condition for iron nuclei to traverse the supernova envelope without 
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significant losses is X < 100 g cm 2 . At late times compared with t e = ro/v e , this 
transparency condition gives t > t tl ~ 1.3 x 10 7 M\E^ 2 s. 
The evolution of the maximum energy is thus given by 



Z 2& -&L3 ^| 



scrw ~ 4 x io 2 vrri ev . (los) 



The maximum energy decreases as the source evolves, and so the condition to produce 
ultrahigh energy CRs is that E CR exceeds the needed energy when the envelope becomes 



transparent 175 



E C K(t tI ) > E 20 , (109) 

or equivalently, 

3000 



Bll 2 [AZ W E^ - 0.13 M 1 B 13 E^ 2 }^ 



B -1 . (110) 



Eq. ( |110| ) translates into upper bounds on the surface magnetic field strength and the 
star initial spin period Pj = 2tt/uj^ 

31Z 26 E^ 2 

and 

P i <8nB 1 1 i 2 Z 26 E w 1 . (112) 

For Mi = 2 and E 20 = E 51 = Z 26 = 1, Eq (|ITT|) gives B 13 < 15.4, whereas Eq. (|TT2D 
leads to Pi ^ 10 ms, not very restrictive values for a young neutron star. Moreover, 
for a reasonable pulsar production rate, the CR-flux on Earth is easily consistent with 
observation, i.e., e^4x 10~ 7 ||175| . 



B. Centaur us A 

Centaurus A (Cen A) is the nearest active galaxy, at a distance of ~ 3.4 Mpc ||333 |. 
It is a complex FRI radio-loud source with galactic coordinates / ~ 310°, b pa 20°, 
and identified at optical frequencies with the galaxy NGC 5128. 25 Different multi- 
wavelength studies have revealed a rather complex morphology: it comprises a compact 



Galactic longitude, I, is measured eastward along the galactic equator from 0° (in Sagittarius) to 
360°. The galactic equator is defined by the central plane of the Milky Way. Galactic latitude, b, 
is the angular distance from the galactic equator (b — 0°) toward either north (b = 90°) or south 
(b = —90°) galactic pole. Note that before August 1958 a different, now obsolete, system was used. 
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core, a jet also visible at X-ray frequencies, a weak counterjet, two inner lobes, a kpc- 
scale middle lobe, and two giant outer lobes. The jet would be responsible for the 
formation of the northern inner and middle lobes when interacting with the interstellar 
and intergalactic media, respectively. There appears to be a compact structure in the 
northern lobe, at the extrapolated end of the jet. This structure resembles the hot 
spots such as those existing at the extremities of FRII galaxies. However, at Cen A 
it lies at the side of the lobe rather than at the most distant northern edge, and the 
brightness contrast (hot spot to lobe) is not as extreme |[334| . 



In order to ascertain the capability of Cen A to accelerate particles to ultrahigh 
energies, one first applies the Hillas criterion [|138|| (see Fig. [TO]). Low resolution po- 



larization measurements in the region of the suspected hot spot give fields as high 
as 25 fiG [|334|| . However, in certain regions where measurements at both high and 



low resolution are available, the S-field amplitude at high resolution can be seen to 
be twice that at low resolution. The higher resolution can reveal amplification in the 



post-shock region ||335|| , yielding 5-fields possibly as high as 50 — 60 /iG ||336| , |337 [ 



The radio-visible size of the hot spot can be directly measured from the large scale 



map ||338|| , giving i? H s — 2 kpc. The actual size can be larger by a factor ~ 2 because 
of uncertainties in the angular projection of this region along the line of sight. 26 All in 
all, if the magnetic field of the hot spot is confined to the visible region, the limiting 
energy imposed by Eq. (|34"1) is ~ 2 x 10 20 eV. Estimates of the radio spectral index 
of synchrotron emission in the hot spot and the observed degree of linear polarization 
in the same region suggests that the ratio of turbulent to ambient magnetic energy 



density in the region of the shock is u ~ 0.4 ||339||. The jet velocity is model dependent: 



possible values range from ~ 500 km s _1 to 0.99c [ [334j| . For FRI galaxies, the ratio 
of photon to magnetic energy densities, a, is expected to be C 1. Figure |T7| shows 
curves of constant ultrahigh energy in the f3- ]et - a plane, for B = 60 /zG, R = 4 kpc, 
according to Eq. (|4]). Since the range of values for a, and the jet velocity in units 
of c conform to expected values, it is plausible that Cen A can accelerate particles to 
energies > 10 20 eV. 

Recent observations of the gamma ray flux for energies > 100 MeV by EGRET [340 ] 27 



allow an estimate L 7 ~ 10 41 erg s 1 for the source. 28 This value of L 7 is consistent 
with an earlier observation of photons in the TeV-range during a period of elevated 



26 For example, an explanation of the apparent absence of a counterjet in Cen A via relativistic 
beaming suggests that the angle of the visible jet axis with respect to the line of sight is at most 



36° [334 1, which could lead to a doubling of the hot spot radius. It should be remarked that for a 



distance of 3.4 Mpc, the extent of the entire source has a reasonable size even with this small angle. 

27 The Energetic Gamma Ray Experiment Telescope (EGRET) was on board the CGRO mission. 

28 Note that the received radiation is negligibly affected by interactions with the various radiation 



backgrounds [194] 
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activity ||341| , and is considerably smaller than the estimated bolometric luminosity 
-^boi ~ 10 43 erg s~ 1 ||333|| . Data across the entire gamma ray bandwidth of Cen A is given 
in Ref . ||342|| , reaching energies as high as 150 TeV 1 343|| , though data at this energy await 
confirmation. For values of B in the /xG range, substantial proton synchrotron cooling 
is suppressed, allowing the production of high energy electrons through photomeson 
processes. The average energy of synchrotron photons scales as _E 7 ~ 0.29i? 7 [344]. 



With this in mind, it is straightforward to see that to account for TeV photons, Cen 
A should harbor a population of ultra-relativistic electrons with E ~ 6 x 10 18 eV. We 
further note that this would require the presence of protons with energies between one 
and two orders of magnitude larger, since the electrons are produced as secondaries. 29 
There are plausible physical arguments [|176|, |193j| as well as some observational 



29 Consecutive factors of ~ 2 energy loss occur in the processes pj — > Ntt°, 
Eq.(|39|) then implies proton energies of ~ 10 20 eV for 100 TeV photons. 



77, 7 
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reasons [|345|| to believe that when proton acceleration is being limited by energy losses, 
the CR luminosity Lqr ~ L 1 . Defining e, the efficiency of ultra high energy CR 
production compared to high energy 7 production - from the above, e ~ 1 - and using 
equal power per decade over the interval (E min , E max ), the source luminosity is found 
to be |3T9| 

E 2 dNl +n 6.3eL 4 il0 52 eV/s /11oN 
dEdt \n(E max / E min ) 

where L 41 = luminosity of Cen A/10 41 erg s -1 and the subscript "0" refers to quantities 
at the source. 

It is straightforward to see, using Eqs. ([)4]) and fl9"8|) that for E = 10 19 eV, B = 
0.5 nG, I = 0.5 Mpc, the diffusive distance traveled ctd = 50 Mpc 3> d = 3.4 Mpc 
is in very good agreement with the diffusion condition given in Eq. ( |10U| ). For higher 
energies, the validity of the diffusive approach must be checked on a case-by case 
basis. For these purposes, in the case of a continuously emitting source, the definition 
of a diffusion time is somewhat arbitrary. In the following discussion we adopt as 
the diffusion time r^, a choice partially consistent with simulations (see Sec.IV-E). 
Moreover, one can easily check that for 3.4 Mpc the diffusion time of any proton with 
energy above the photopion production threshold is always less than the GZK-time, 
and consequently energy losses can be safely neglected. The density of protons at 
the present time t of energy E at a distance r from Cen A, which is assumed to be 
continuously emitting at a constant spectral rate dNQ +n /dEdt from time t on until the 
present, is ||346 | 

dn(r t) dN p+n 1 f l e -r 2 /4D{t-t>) (jN p+n 1 
dn{r, t) _ dJ\ I dt , e = dp_ 1 ^ (m) 



dE dEdt [4:7rD(E)} 3 / 2 J ton (t-t'f/ 2 dEdt47iD(E)r 

where D(E) stands for the diffusion coefficient, x = 4DT on /r 2 = T on /r D , T on = t — t on , 
and 

1 f°° du 



I(x) = -= -= e-" . (115) 
V 71 " Ji/x Vu 

For T on — > 00, the density approaches its time- independent equilibrium value n eq , while 
for T on = t d , n/n eq = 0.16. 

To estimate the power of Cen A, one can evaluate the energy-weighted approximately 
isotropic proton flux at 1.5 x 10 19 eV, which lies in the center of the flat "low" energy 
region of the spectrum, 

Fr F 2 JN p+n 

E ^ - dElt " 7.6xl0- e£ »I.Vm-^«-. (116) 

In Eq. ( [L16| ) we have used the fiducial values of B and t as given in the previous 
paragraph, and set E m[n = 1 x 10 19 eV, E max = 4 x 10 20 eV, assuming (as discussed 
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in Sec.III-B.2) that the shock structure can be more extended than the visible region 



in the non-thermal radioemission. As noted by Farrar and Piran [ [319 



by stretching at most the source parameters, that the "low" energy flux from Cen A 
could be comparable to that of all other sources in the Universe. To this end, first fix 
e L41 1 = 0.40, after comparing Eq. ( gTg ) to the observed CR-flux: E 3 J ohs {E) = 10 24 ' 5 
eV 2 m~ 2 s~ 4 sr _1 |20|. Next, eL 41 ~ 1, determines / ~ 0.40, and consequently the 
required age of the source T on to be about 400 Myr. There is still a great uncertainty 
on jet evolutionary behavior. Experimental approaches and theoretical studies suggest 
ages that run from a few tens of Myr, up to half the Hubble time in some extreme 
cases ||186||. Moreover, T on ~ 400 Myr is consistent with FRI lobe ages ||345|| . To 



maintain flux at the "ankle" for the same T on , one requires an approximate doubling 

of L CR at 5 x 10 18 eV. Because of the larger diffusive time delay at this energy, this 

translates into an increased luminosity in the early phase of Cen A. From Eq. (j39|), 

the associated synchrotron photons are emitted at energies < 30 MeV. The increase 

in radiation luminosity in this region is not inconsistent with the flattening of the 

spectrum observed at lower energies | )347| , |348| |. 

Let us turn now to the discussion of anisotropy. This can be found by computing the 

incoming current flux density DVn as viewed by an observer on Earth, and one finds 

for a continuously-emitting source a distribution ~ (1 + a cos 9) about the direction of 

the source, where 9 is the angle to the zenith and 

2D(E) I' , . 

a = — ^ ■ -. 117 
cr I 

Here, 

1 f°° 

I\x) = —= \ du^i e~ u , (118) 
V 71 " Jl/x 

with x = T on /rr), and / as defined in Eq. ( |115| ) [|346f . For our choices of B and £, and 
T on = 400 Myr, we find for E = 10 19 eV (E = 10 20 eV) that a = 0.04 (a = 0.07). 

Neutrons at the highest energies could survive decay and produce a spike in the 
direction of the source (see Appendix E). However, those that are able to decay will 
beget secondary proton diffusion fronts with asymmetry parameters given by 

2D(E) I" . . 

a = ^^-— , 119 
cr I 



where 



(120) 

and k = X(E)/r, where X(E) = 0.9 E 2 o Mpc is the neutron decay length. In spite of 
the complicated nature of Eq. ( |120p , the results for a are very similar to the ones for 
the primary diffusion front given above. 
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FIG. 18: The angular image in equatorial coordinates, averaged over all 20 magnetic field re- 
alizations of 5000 trajectories each, for events above 40EeV (upper panel) and above 100 EeV 



(lower panel), as seen by a detector covering the whole sky, for the case suggested in Ref. |319(| 
corresponding to B = 0.3 /iG, and the source Cen A. The grey scale represents the integral 
flux per solid angle. The solid line marks the Super-Galactic plane. The pixel size is 1°; 
the image has been convolved to an angular resolution of 2.4° corresponding to AGASA. 
Published in Ref. |34 



Detailed Monte Carlo simulations 349, 350 



however, reveal that the predicted 
distributions of arrival directions are more anisotropic than current data even at 10 19 eV 
(see Fig. [18].). In particular, for a spectral index 7 = 2.4 and a maximum injection 
energy of 10 21 eV, the angular power spectrum shows a 3 a quadrupole deviation from 
observation, and the auto-correlation function is not consistent with the clustering 
at small scale observed by AGASA ||350|| . The discrepancy between phenomenology 
and Monte Carlo simulations could arise from a number of factors. The simulations 
are considerably more detailed than the phenomenological arguments; for instance 
effects of photopion production are taken into account. Furthermore, they simulate 
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the subtleties in transitions from Kolmogorov to Bohm diffusion. On the other hand, 
the Monte Carlo simulations rely on a statistical description of the magnetic field, 
in which the field is arranged into domains or "cells" with a characteristic coherence 
length. There is some uncertainty in the sensitivity to the geometry of the cell, which is 
reflected in the presence or absence of resonant modes in the spectrum. For instance, 
if one considers coherence domains generated in some Gaussian manner (say with 
spatial spread R) rather than as small cubes, then fc = 1/R instead of k = 27r /£, 
and for R ~ 0.7 Mpc, the diffusive resonant condition is satisfied even for the highest 
observed energies. Whichever point of view one may find more convincing, it seems 
most conservative at this point to depend on experiment (if possible) to resolve the 
issue. In particular, PAO will be in a unique position to search the sky in the direction 
of Cen A. 



M87 



20 Mpc, M87 is the dominant radio galaxy in the Virgo cluster 
||351|| . This powerful FRI has been under suspicion to be the 



At a distance of ~ 
(I = 282°, b = 74°) 

primary source of ultrahigh energy cosmic rays for a long time [ [307| , |352jl . The emission 
of synchrotron radiation with a steep cutoff at frequencies about 3 x 10 14 Hz from its 
radiojets and hot spots [ 353 . 354|| implies an initial turbulence injection scale having 
the Larmor radius of protons at 10 21 eV (see Sec.III-B.2). Hence, because of its relative 
proximity to Earth, M87 becomes a potential candidate to account for the observed 
events above 10 20 eV. The major difficulty with this idea is the observation of the nearly 
isotropic distribution of the CR arrival directions. One can again argue that the orbits 
are bent. However, the bending cannot add substantially to the travel time, otherwise 
the energy would be GZK-degraded. An interesting explanation to overcome this 
difficulty relies on a Galactic wind, akin the solar wind, that would bend all the orbits 
of the highest energy CRs towards M87 [|355| , [356|| . Indeed, it has long been expected 
that such a kind of wind is active in our Galaxy ||357| , [358| , [359|| . In the analysis of ||355|| , 
it was assumed that the magnetic field in the Galactic wind has a dominant azimuthal 
component, with the same sign everywhere. This is because in a spherical wind the 
polar component of the magnetic field becomes negligible rather quickly, decaying like 
1/r 2 , and thus the azimuthal part of the magnetic field quickly becomes dominant, 



with ~ sin#/r i n polar coordinates [ 360 1. Under these considerations one is left 



with two degrees of freedom: the strength of the azimuthal component at the location 
of the Sun, and the distance to which this wind extends. Recent estimates suggest that 
the magnetic field strength near the Sun is ~ 7 fiG ||361 . The second parameter is 



more uncertain. Our Galaxy dominates its near environment well past our neighbor, 
M31, the Andromeda galaxy, and might well extend its sphere of influence to half way 
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FIG. 19: Directions in polar projection of the highest energy CR events when they enter the 
halo of our Galaxy. The two highest energy events have two entries, under the assumptions 
that they are either protons or helium nuclei (filled black symbols). The cross indicates 
the position of M87 and the grey band stands for the super-Galactic plane. Published in 
Ref. |5|. 



to M81. This implies an outer halo wind of ~ 1.5 Mpc. With this in mind, the mean 
flight time of the protons in the Galaxy is ~ 5.05 x 10 6 yr <C t s , the time for straight 
line propagation from M87 [ 362| |. Figure [TP| shows the directions where the 13 highest 
energy CR events point towards when they leave the halo wind of our Galaxy. Note that 
except for the two highest energy events, all other events can be traced back to within 
less than about 20 degrees from Virgo. If one takes into account the uncertainty of the 
magnetic field distribution, all events are consistent with arising originally from Virgo. 
Moreover, if one assumes that the two highest event are helium nuclei, all 13 events 
point within 20° of Virgo. Note that a bending of 20° could be easily accommodated 
with a fine tuning of the magnetic field strength within the super-Galactic sheet from 
here to Virgo (more on this below). Criticisms of this model ||363|| have been addressed 
in [EM . 



We now focus attention on the power of M87. To this end, we analyze another 
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interesting limit of Eq. (|92|) in which diffusion is found to be extremely small. In such 
a case, the second term in the r.h.s. can be neglected, and the solution for a single 



point source reads ||228 

1 f°° 

n( - E ^ = We) J e ^ E ^ ) dE ^ ( m ) 

where 



Ie b{E) 

and Eq is the energy of the nucleon when emitted by the source. Using Eq. (p2|) we 
model the injection spectrum of a single source located at t from the observer by 

Q(E,t) = KE-~<5(t-t ), (123) 

where k is a normalization constant and for simplicity we consider the distance as 
measured from the source, i.e., t = 0. At very high energies, the total number of 
particles at a given distance from the source is 

n(EAa W)[ E ~ al5 { t -^) d ^ (124) 

or equivalently, 

n(E,t) ^ KE-^e- { ^ 1)ct , (125) 



where C is given in Eq. (|5lj|). Thus, the injection spectrum is uniformly damped by a 
factor that depends on the proximity of the source. 

At low energies, in the region dominated by baryon resonances, the parameterization 
of b(E) does not allow a complete analytical solution. However using the change of 
variables, 

r E 9 dE 

with Eq = £(E,t) and dt = dEQ/b(E ), one easily obtains, 

n(E,t) = — / aE,t)-~< 5(t-t) b[S(E,t)] dt, (127) 
b{E) Jq 



and then, the compact form, 



n(E,t) = ^E ^b(E ), (128) 



where the relation between E and E is given in Eq. (0). Now, to take into account the 
extension of the cluster, we assume that the concentration of potential sources at the 
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center is higher than that in the periphery, and adopt a spatial gaussian distribution. 
With this hypothesis, the particle injection rate into the intergalactic medium is given 
by 



E~ J f-(T-t 



DC 



Q{E,t) = Kj -^^S(t-T)exp\ K \ 2 0) \dT (129) 
'2 7r a I 20 



\2 



A delta function expansion around t , with derivatives denoted by lower case Roman 
superscripts, 

5(t — T) = 5{t - t ) + S {i) (t - to)(T - t ) + ^ u \t - t )(T - t ) 2 + ... (130) 
leads to a convenient form for the injection spectrum, which is given by, 

Q(E, t') = kE^ [ 8{t' - t ) + ^ S^(t' - t ) + ^ 6W(t' -t ) + ...]. (131) 

By replacing the above expression into Eq. ( |121| ) one obtains ||234|| 

E^b(E Q ) f e 2 A 2 e- 2B / E ° , . a 4 A 4 e- 4B /^ . . . S\ , . 
n = k ° b{ ^ 0> |l + 21 + 41 F2 ^°) + ' ( 132 ) 

where 

F^Eo) = 2B 2 E 2 + (2 - St)^ 1 + (1 - 7 ) 2 , (133) 

F 2 (E ) = 24B 4 £ " 4 + (4 - 50^)B 3 E 3 + (35 7 2 - 25 7 + S)^ 2 ^- 2 

+ (-IO7 3 + 20 7 2 - 15 7 + A)BE- X + (1 - 7 ) 4 , (134) 

and A and B as given in Eq. (|55|) . Figure ^0| shows the evolved energy spectrum 
of nucleons assuming a cosmologically homogeneous population of sources (universal 
hypothesis) |230|1 together with the modified spectrum as given by Eq. (|132[) for an 



extended source modeled by a Gaussian distribution of width 2 Mpc at a distance of 
18.3 Mpc. The extended source spectrum was normalized in the "low" energy region 
with the spectral index, 7 = 3.27, obtained from a maximum likelihood fit to the Fly's 



Eye data [p65|| . From Fig. £20] one can see that the proton component of the Virgo cluster 
partially reproduces AGASA data, but apparently cannot account for the super-GZK 
Fly's Eye event. However, lensing effects may come to the rescue [ [366|| . The magnetic 
field in the Galactic wind strongly magnifies the flux, reducing the power requirements 
of the source. Specifically, with fine-tuning of the source direction relative to the 
symmetry axis of the wind, M87 could be as high as > 10 2 times more powerful than if 
unlensed at energies below £/Z ~ 1.3 x 10 20 eV [ |366|| . Moreover, as noted in ||367|| , the 
super-Galactic plane sheet can focus ultrahigh energy CRs along the sheet. Hence, the 
particles would arrive at the boundary of our Galactic wind with the arrival directions 
described by an elongated ellipse along the super-Galactic plane sheet. This in turn 
would allow magnetic lensing to produce many hot spots, just as sunlight dapples the 
bottom of a swimming pool, and so easily explain the apparent CR-clustering. 
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FIG. 20: The CR flux as observed by AGASA (square) and Fly's Eye (triangle) experiments. 
The solid line indicates the CR spectrum assuming the universal hypothesis (UH), i.e., ho- 
mogeneous source distribution, whereas the dotted line indicates the evolved spectrum from 
the Virgo cluster (VC). Published in Ref. p64|] . 



D. Starbursts 



If the super-GZK particles are heavy nuclei from outside our Galaxy, then the 
nearby (~ 3 Mpc p8| ) starburst galaxies M82 (I = 141°, b = 41°) and NGC 253 
(/ = 89°, b = —88°) would probably be the sources of most ultrahigh energy CRs ob- 
served on Earth. Starbursts are galaxies undergoing a massive and large-scale star 
formation episode. Their characteristic signatures are strong infrared emission (orig- 
inating in the high levels of interstellar extinction), a very strong HH-region-type 
emission-line spectrum (due to a large number of O and B-type stars), and consid- 
erable radio emission produced by recent supernova remnants (SNRs). Typically, the 
starburst region is confined to the central few hundreds of parsecs of the galaxy, a region 
that can be easily 10 or more times brighter than the center of normal spiral galaxies. 
In the light of such a concentrated activity, the existence of galactic superwinds is not 
surprising ||368| . 

Galactic-scale superwinds are driven by the collective effect of supernovae and mas- 
sive star winds. The high supernovae rate creates a cavity of hot gas (~ 10 s K) whose 
cooling time is much greater than the expansion time scale. Since the wind is suf- 
ficiently powerful, it can blow out the interstellar medium of the galaxy preventing 
it from remaining trapped as a hot bubble. As the cavity expands a strong shock 
front is formed on the contact surface with the cool interstellar medium. The shock 
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velocity can reach several thousands of kilometers per second and ions like iron nu- 
clei can be efficiently accelerated in this scenario, up to ultrahigh energies, by Fermi's 
mechanism ||173|| . 



In a first stage, ions are diffusively accelerated at single supernova shock waves 
within the nuclear region of the galaxy. Energies up to ~ 10 15 eV can be achieved in 
this step ||f 47|| . Heavy nuclei are not photodissociated in the process despite the large 



photon energy densities (mostly in the far infrared) measured in the central region 
of the starburst. The escape of the CR outflow is convection dominated. In fact, 
the presence of several tens of young SNRs with very high expansion velocities and 
thousands of massive O stars, with stellar winds of terminal velocities up to 3000 km 
s -1 , must generate collective plasma motions of several thousands of km per second. 
Then, due to the coupling of the magnetic field to the hot plasma, the magnetic field is 
also lifted outwards and forces the CR gas to stream along from the starburst region. 
Most of the nuclei escape in opposite directions along the symmetry axis of the system, 
as the total path traveled is substantially shorter than the mean free path [[173 . 



Once the nuclei escape from the central region of the galaxy, with energies of 
~ 10 15 eV, they are injected into the galactic-scale wind and experience further accel- 
eration at its terminal shock up to 10 20 eV. For this second step in the acceleration 
process, the photon field energy density drops to values of the order of the microwave 
background radiation, as we are now far from the starburst region, and consequently, 
iron nuclei are safe from photodissociation. In terms of parameters that can be deter- 
mined from observations, the nucleus maximum energy is given by ||1 73 



E max ^\zeB^T on , (135) 

where E sw ~ 2.7x 10 42 erg s -1 is the superwind kinetic energy flux and M = 1.2M© yr _1 
is the mass flux generated by the starburst |[368|| . The age T on can be estimated from 
numerical models that use theoretical evolutionary tracks for individual stars and make 
sums over the entire stellar population at each time in order to produce the galaxy 



luminosity as a function of time [|369|[ . Fitting the observational data, these models 
provide a range of suitable ages for the starburst phase that goes from 50 Myr to 
160 Myr ||369 |. These models must assume a given initial mass function (IMF), which 



usually is taken to be a power-law with a variety of slopes. Recent studies have shown 



that the same IMF can account for the properties of both NGC 253 and M82 ||370 
Besides, a region referred to as M82 "B" near the galactic center of M82, has been 
under suspicion as a fossil starburst site in which an intense episode of star formation 



occurred over 100 Myr ago P7TI |372|[. The derived age distribution suggests steady, 



continuing cluster formation at a modest rate at early times (> 2 Gyr ago), followed by 
a concentrated formation episode 600 Myr ago and more recent suppression of cluster 
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formation. In order to get some estimates on the maximum energy, let us assume 



B ~ 50 fiG, a choice consistent with observation ||373|| . Inserting all of these figures in 



Eq. ( |135| ), already for T on = 50 Myr one obtains for iron nuclei 

^ e ax >10 2 ° eV. (136) 

Now, we can use the rates at which starbursts inject mass, metals and energy into 
superwinds to get an estimate on the CR-injection spectra. We once again use e 
to denote the efficiency of ultrahigh energy CR production by the superwind kinetic 
energy flux. Using equal power per decade over the interval 10 18,5 eV < E < 10 20,6 eV, 
we obtain a source CR-luminosity 

E 2 dN _. or „ 



dEdt 



3.5elO M eV/s. (137) 



The assumption that the giant air showers with E > 10 20 eV were triggered by heavy 
nuclei implies ordered extragalactic magnetic fields B n Q < 15 (at least in the outskirts 
of the Galaxy), or else nuclei would be captured in magnetic subdomains suffering 
catastrophic spallations. Moreover, even for fields O (nG), CR nuclei with energies 
E c < 10 18 ^Mpc ZB n Q eV remain trapped inside cells of size ^Mpc, attaining efficient 
diffusion, with 



D{E) « 0.048 ° mpc MpcVMyr. (138) 
V Z B nG J 

The power of the starbursts can then be estimated by straightforward generalization 
of the procedure discussed in Sec.V-B. Specifically, first evaluate the energy-weighted 
approximately isotropic nucleus flux at 10 19 eV, 

E * J W = T^mwS'' " 23 * ">^eVV- W , (139) 

where h = Ims2 + hcc 253, B nG = 15, £ Mpc = 0.5, and (Z) = 20. Then, fix 

eh = 0.013, (140) 

after comparing Eq.( |139|) to the observed CR-flux. Note that the contribution of 7m82 
and /ngc 253 to I* critically depends on the age of the starburst. Figure ^l] shows the 
relation "starburst-age/superwind-efficiency" derived from Eq. ( |140| ), assuming that 
both M82 and NGC 253 were active for 115 Myr (e w 10%). Beyond this epoch, 
CR-emission must be associated with M82 "B" . 

Above > 10 20 ' 2 eV iron nuclei do not propagate diffusively. Moreover, the CR- 
energies get attenuated by photodisintegration on the microwave background radiation 
and the intergalactic infrared background photons. In the non-diffusive regime, the 
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FIG. 21: Age of the starbursts as a function of the efficiency of CR-production, e. Published 
in 



EH 



accumulated deflection angle 9(E) from the direction of the source can be estimated 
using Eq. (|90f) . Therefore, if B ~ 15 nG all directionality is lost. The resulting time 
delay with respect to linear propagation is given by Eq. (ffil), and the total travel time 
is 

t « - c (1 + ±0 2 ) . (141) 

As an interesting exercise one can apply these considerations to the highest energy 
Fly's Eye event. Including statistical and systematic uncertainties, the energy of this 
event is 3.2 ± 0.9 x 10 20 eV. Eqs. (|63f) and (^) relate the uncertainty in energy to the 
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uncertainty in the attenuation time 



£ V -^20 / \ E J 

From these considerations, we find that the upper limit on the transit time for a nuclear 
candidate for the highest energy Fly's Eye event is ~ 6 x 10 14 s. The arrival direction 
of the highest energy Fly's Eye event is 37° from M82 [ |172j] . With d ~ 3 Mpc and 
9 = 37°, one finds from Eq. ( |141| ) a transit time t ~ 3.4 x 10 14 s, well within the stated 
upper limit. 

For average deflections of 60°, the time of flight is ~ 3.9 x 10 14 s, and consequently the 
CR-spectrum falls off sharply in the neighborhood of the GZK limit. A slight anisotropy 
should arise just before this cutoff. As the candidate sources in the northern and 
southern hemispheres are presumably at different distances, 30 a north-south asymmetry 
should also eventually emerge. It is rather difficult to assess whether events with 
energies > 10 20 ' 5 eV are plausible, because the maximum energy depends strongly on 

''"delay • 

The photodisintegration process results in the production of nucleons of ultrahigh 
energies with the same Lorentz factor as the parent nucleus. As a consequence, the 
total number of particles is not conserved during propagation. However, to a very 
good approximation one can treat separately the evolution of the heaviest fragment 
and those fragments corresponding to nucleons emitted from the traveling nuclei. With 
this in mind, the evolution of the differential spectrum of the surviving fragments is 
given by Eq. fl9"2"|). Now, generalizing the formalism discussed in the previous section 
for the case of a single source located at to from the observer, one obtains a slightly 



modified expression for the evolution of the spectrum [[25S 



n(E,t)~^^ } (143) 



where the relation between Eq and E is given in Eq. (|66|). Putting all this together, 
the energy-weighted flux beyond the GZK-energy due to a single M82 flare 

E *J(E) = . E n9 y% ~*' /56 ~ 2.7 x lO^oee^^eV^^s^sr- 1 , (144) 
[Andy dE dt 



is easily consistent with observation [ 20[ . It is important to stress that these phe- 



nomenological arguments are in agreement with Monte Carlo simulations ||375|| . In 5 



years of operation PAO will provide enough statistics to test the starburst anisotropy 
predictions ||299|| . 



The location of NGC 253 is still subject to large uncertainties [374] 



80 



VI. GZK EVADING MESSENGERS 



Several attempts have been made to explain the high end of the spectrum as a man- 
ifestation of some kind of physics beyond the SM. These suggestions became really 
interesting after Farrar and Biermann [ 376 ] reported an intriguing correlation between 
the arrival directions of 5 extremely high energy CRs and compact radio quasars at 
high redshifts (z = 0.2 — 0.3). However, with the present data, such "evidence" for 
directional correlation is still a subject of debate ||377| , |378| , |379| , |380|| . The most eco- 
nomical proposal to sneak away the GZK-cutoff involves a familiar extension of the 
SM, namely, neutrino masses. As noted by Weiler [351, 382|| , neutrinos can travel over 
cosmological distances with negligible energy loss and could, in principle, produce Z 
bosons on resonance through annihilation on the relic neutrino background within a 
GZK distance of Earth. If this were the case, the highly boosted decay products of the 
Z could be observed as super-GZK primaries, since they do not have to travel cosmo- 
logical distances to reach us, and would be pointing directly back to the source. This 
proposal requires very luminous sources of extremely high energy neutrinos through- 
out the Universe. Another possibility in which the ultrahigh energy CR should point 
directly to its source assumes that the CR is a new, neutral, stable or very long lived 
super symmetric hadron of mass a few GeV |383| . From all the candidates en vogue, 
the light glueballino seems the most plausible, because it can be efficiently produced 
in pp collisions, it is not strongly absorbed by the CMB, and produces extensive air 
showers with longitudinal development very similar to those observed |384|. An even 



more radical proposal postulates a tiny violation of Lorentz invariance, such that some 
processes become kinematically forbidden [385, 386 , 387| , 388 , 389 1. In particular, 



photon-photon pair production and photopion production may be affected by Lorentz 
invariance violation. Hence, the absence of the GZK-cutoff would result from the fact 
that the threshold for photopion production "disappears" and the process becomes 
kinematically not allowed. 

For models that rely on GZK-evading messengers, the accumulation factor f acc de- 
pends on the details of the model. In the case of messengers which can induce showers 
across the entire energy spectrum, one expects enhancement on the low energy side 
only from the baryonic component, and f acc depends on the interaction length in the 
CMB and on the energy spectra at the different sources. For messengers whose at- 
tenuation length is comparable to the horizon, and which do not shower at the lower 
energies, f acc -C 1. In what follows, we give a somewhat more detailed account of these 
general considerations. 
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A. Super-GZK CRs from the edge of the Universe 



The neutrino is the only known stable particle that can propagate through the CMB 
essentially uninhibited even at the highest CR energies. Specifically, the corresponding 
vv annihilation mean free path on the cosmic neutrino background, 

\ v = (n„ ow)" 1 « 4 x 10 28 cm , (145) 

is just above the present size of the horizon (recall that H^ 1 ~ 10 28 cm) ||390|| . One 
may then entertain the notion that neutrinos are indeed the super-GZK primaries. 
However, in the SM model a neutrino incident vertically on the atmosphere would pass 
through it uninhibited, never initiating an extensive air shower. Consequently one has 
to postulate new interactions so that these neutrinos acquire a strong cross section 
above 10 20 eV. 

The idea that neutrinos are the super-GZK primaries was introduced by Beresinsky 



and Zatsepin in 1969 |[249| . Since then, explicit models have been devised in which the 
neutrino-nucleon cross section is enhanced by some new physics beyond the electroweak 
scale. One interesting proposal is that leptons are bound states of dual QCD gluons, 
and therefore can interact "strongly" with all partons in the nucleon attaining an 
effective cross section comparable to the geometric nucleon cross section ||391| , |392 



More recently, low scale gravity models have provided new impetus to the 
hypothesis of strongly interacting neutrinos. In particular, the neutrino nucleon cross 
section can be made large enough due to a quick rise of the density of states, possibly 
increasing exponentially, p ~ eV^*, above a characteristic energy scale s | ]395| j. A 
specific implementation of this idea is given in theories with n additional large compact 



dimensions and TeV scale quantum gravity [ 396 . A detailed discussion of this last 
proposal is left until Sec.VII-A. For now, we will focus on the Z-burst model, in which 
everyday weekly interacting neutrinos explain the super-GZK primaries. 



B. Z-burst 



To ascertain whether the Z-burst hypothesis is viable, we would like to know what 
neutrino energy and flux would actually be required to generate observed highest energy 
CRs. This in turn depends on the neutrino mass, the distribution of thermal neutrinos 
and the kinematics of the Z decay products. These three points are discussed in the 
following paragraphs. At the end of the section, we comment on the compatibility of 
the Z-burst hypothesis with current EAS data and plausible astrophysical sources. 

In recent years, stronger and stronger experimental evidence for neutrino oscillations 
has been accumulating. Certainly, this evidence - which is compelling for atmospheric 
neutrinos |97], |9| g9|, [400|, gOT), strong for solar neutrinos g0|, |03|, [404], |0|, |406 
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407|| , and so far unconfirmed for neutrinos produced in the laboratory and studied by 
various groups |08|, g0|, gTO], gll], gT|, III' B HI HI ~ would extend the SM 
by requiring non-vanishing neutrino masses and mixing. A precise estimate of the 
^-masses and mixing parameters would be an important clue to any physics beyond 
the SM. However, neutrino oscillations are sensitive to the mass (squared) splittings 
Am% = m?, — vr?,, , and so only a lower bound on the mass of the heaviest neutrino can 

ij Ui Vj i <J 

be obtained. 31 Strictly speaking, using the atmospheric mass splitting in a 3 neutrino 
flavor scenario one obtains 



ni 



> VAmL >0.04eV. 



(146) 



whereas by using the LSND mass splitting in a 4 flavor scenario the limit is shifted an 
order of magnitude, i.e., 



rn 



;/ 4 



> 



AKsnd ^0.4eV 



(147) 



For recent surveys see |[417| , [418| , [419|| . The search for mass imprints in the endpoint 
spectrum of tritium f3 decay combined with experimental constraints from oscillations 
yield upper bounds on the mass of the heaviest neutrino, 



m 



< 



m 2 p + Ami 



< 9 5pV 

atm ~ c v J 



in a 3 flavor, and 



< 



S + ^ m LSND <3.8eV. 



(148) 



(149) 



in a 4 flavor scenario [|420| . Additional information on the neutrino mass scale can be 



derived through cosmological and astrophysical observations ||421 |. For instance, anal- 
yses of galaxy clustering, including recent CMB measurements and other cosmological 
constraints, give an upper bound 



E 



m n < 1. 



- 4.4 eV 



(150) 



on the sum of the neutrino masses 142 



423] P|, p| (4261, p? 



Big bang cosmology predicts that the Universe is filled with a shallow degenerate 
Fermi sea of neutrinos | 428 |. Neutrino decoupling occurred at a temperature of 1.9 K, 
just one second after the bang. This thermal neutrino background (TNB) has an 
average number density of 



(n Ui ) = (n Pi ) = 22 ( n ^° ~ 56 cm ~ 



(151) 



CMB 



i, j run over the possible neutrino species. 
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per neutrino flavor. The dominant interaction mode of extremely high energy neutrinos 
with the TNB is the exchange of a W boson in the t-channel [yiVj — > klj), or of a 
Z° boson in either the s-channel [yiVj — » //) or the t-channel [y%Vj — > i^Pj). Here, Z 



denotes a charged lepton, / a charged fermion, and i ^ j for the first reaction | 390 |. 
Therefore, one expects neutrinos within a few Z widths of the right energy, 

E^ = ^L = A (—)x 10 21 eV , (152) 



2m v . \ m 



to annihilate with the TNB into hadrons at the Z-pole with large cross section, 

/ds 
— 2 < 7 ann (s) = 27rV2G F ~40.4nb, (153) 
z 



producing a "local" flux of nucleons and photons [[42 9| , |430|| . 32 Remarkably, the energy 



of the neutrino annihilating at the peak of the Z-pole has to be well above the GZK 
limit. At this stage, it is worthwhile to point out that in view of the expected rapid 
decrease of the ultrahigh energy v flux with rising energy, the t-channel W- and Z- 
exchange annihilation processes can be safely neglected. On resonance, the s-channel 
Z-exchange interactions completely overwhelm them. 

The mean energies of the ~ 2 nucleons and ~ 20 7-rays in each process can be 
estimated by distributing the resonant energy among the mean multiplicity of 30 sec- 
ondaries. The proton energy is given by 

(E p ) - - 1.3 ( —) x 10 20 eV, (154) 



60 m,,. V m v 



whereas the 7-ray energy is given by 



~ i - 7 (& ) X 1020 eV ' < 155 > 



3 



The latter is a factor of 2 smaller to account for the photon origin in two body tt° 
decay. 

If the neutrino sources are randomly distributed in space, then the total rate of 
super-GZK events induced by vv annihilation at the Z pole within a distance D of the 



Earth is 429 



F z ~ E?* F v (a ann ) z Bi J d 3 x £^ = E^ F v (a ann ) z B z h J n v dr . 



(156) 



Gf = 1.16639(1) x 10 5 GcV 2 is the Fermi coupling constant. 
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where F v is the incident neutrino flux evaluated at the resonant energy, B% ~ 0.70 
is the hadronic branching fraction of the Z, and n u is the column number density of 
relic neutrinos. In deriving Eq. (|156[) one assumes that (a ann ) dr n u <C 1. The Z- 
burst rate can be amplified if neutrinos are clustered rather than distributed uniformly 



throughout the Universe [^29, 430 ] . In such a case the probability of neutrinos to 
annihilate within the GZK zone is generally on the order of 1%. 

An exhaustive analysis of the Z-burst parameter space has recently been reported 



by Fodor, Katz and Ringwald (FKR) ||431| , 432] . The analysis includes two possibilities 
for a diffuse background of protons, as distinct from protons resulting from the Z- 
burst itself. On the one hand they analyze the case where all events above 4 x 10 19 eV 
consist of protons which are produced within our Galactic halo or at least within the 
GZK zone. For this case, hereafter referred to as halo background (Halo bk'd), no 
GZK attenuation is considered. On the other hand, they analyze the situation where 
CRs are protons which originate from uniformly distributed extragalactic sources, or 
so-called "extragalactic background" (EG bk'd). In view of the observed distributions 
of arrival directions, the EG bk'd seems to be phenomenologically more realistic. In 
such FKR maximum likelihood fit yields 0.08 eV < m v . < 1.3 eV, while for 

a Halo bk'd 2.1 eV < m n < 6.7 eV, both at the 68 % CL. 

The FKR analysis seems to be pretty insensitive to the precise values of the cosmo- 
logical parameters: h, fl^, and Q\, the normalized Hubble expansion rate, and the 
matter and vacuum energy densities, respectively. The high energy z/-flux is described 
by 

F u (E u ,z)=F u (E u ) (l + z) a , (157) 

where z is the redshift and a characterizes the source evolution. Since neutrinos are 
produced as secondaries in hadronic astrophysical sources, the flux at zero redshift is 
expected to fall off like, F V {E V ) oc E^ 1 , with 7 > 1. The FKR analysis takes into 
account two extreme scenarii: (i) strong 7-ray attenuation, where ultrahigh energy 
photons from the Z-bursts do not contribute to the observed flux [this is certainly the 
case if the radio background is sufficiently large and/or the EGMF is 0(nG)], and (ii) 
minimal universal radio background with vanishing EGMF. 

EGRET 1 4=33 1 measurements of the diffuse 7-background in the energy range between 
30 MeV and 100 GeV strongly constrain the source evolution, yielding a ^ quite 
independently of different assumptions about the radio background. It should be noted 
that a correlation between ultrahigh energy CRs and BL Lac objects at redshifts z > 0.1 



has been claimed [|1 52f| . There is evidence that such a sub-class of AGN has zero or 



negative cosmological evolution. 

The required neutrino fluxes for the Z-burst hypothesis are given in Fig. together 
with existing upper limits and projected sensitivities of present and near future exper- 
iments. Power-law extrapolation of the fluxes F v oc E^ 1 below the resonance energy 
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FIG. 22: Neutrino fluxes, F = g X)f = i (i* 1 ^ + ^H); required by the Z-burst hypothesis for 
the case of a halo and an extragalactic background of ordinary cosmic rays, respectively 
(a = 0, h = 0.71,0m = 0.3, f^A = 0.7, z max = 2). Shown, as two crosses, are the necessary 
fluxes obtained for the case of a strong ultrahigh energy 7 attenuation. The horizontal errors 
indicate the 1 a (solid) and 2 a (dotted) uncertainties of the mass determination and the 
vertical errors also include the uncertainty of the Hubble expansion rate. Also shown are 



upper limits from Fly's Eye [125] on F Ve + Fp e and the Goldstone Lunar Ultrahigh energy 
neutrino Experiment GLUE [136] on Yl a =e ui^a as wen as projected sensitivities of 



AMANDA pH on F v + F Pp and PAO on F Ve + F Ue . Published in Ref. plj . 



with spectral indices 7 ^ 1.5 are excluded by the Fly's Eye data. Furthermore, in the 
Z-burst hypothesis, the flux required for power-law extrapolations with indices 7 > 1 to 
lower energies is larger than the theoretical upper limit from "hidden" hadronic astro- 
physical sources, E 2 F V ~2x 10 7 eV m~ 2 s" 1 sr _1 [ 435 1 . Hidden hadronic astrophysical 
sources are those from which only photons and neutrinos can escape. Therefore, one 
has to invoke sources that have to be opaque to their primary protons, and should 
emit secondary photons (unavoidably produced together with the neutrinos) only in 
the sub-MeV region to avoid conflict with the diffuse 7-ray background measured by 
the EGRET experiment g36|, @37| . 



It is fair to say that nowadays no convincing astrophysical sources are known which 



N6 



can accelerate protons at least up to 10 23 eV, have zero or negative cosmological evo- 
lution, are opaque to nucleons, and emit photons only in the sub-MeV region. 33 It is 
still an open question whether such challenging conditions can be realized in Nature. 
As one can see in Fig |22|, the PAO should give a definite answer to this question. 



C. SUSY U 

Certainly, a novel beyond-SM-model explanation to the high end of the spectrum 
is to assume that ultrahigh energy CRs are not known particles but a new species of 
particle, generally referred to as the uhecron, U [ (440| , [44l|| . The meager information we 
have about super-GZK particles allows a naive description of the properties of the U. 
The muonic content of EASs suggests U should interact strongly. At the same time, 
if uhecrons are produced at cosmological distances, they must be stable, or at least 
remarkably long lived, with mean-lifetime 

where d is the distance to the source and mjj, the uhecron's mass. In addition, to avoid 
photopion production on the CMB, mu ^1.5 GeV, because the threshold energy for 
the reaction Uj — ► Utt, increases linearly with mu, 

E th = m 7T . (159) 

w 

Moreover, to avoid deflections on EGMFs and the consequent energy loss due to pair 
production and other mechanisms, the hadron has to be neutral. Note that the latter 
may not be an essential requirement, depending on the distance to the accelerator. 

It has been known since the early days of Supersymmetry (SUSY) that the gluino 
g, the spin ~ SUSY partner of a gluon, could be the lightest supersymmetric particle 
(LSP), with the possible exception of the gravitino #3/2 [ |442| , [443| , |444j , [445|| . 34 There- 



fore, in a SUSY background the plausible candidates satisfying the U requirements are 
gluino containing hadrons (g-hadrons). 35 QCD sum rules suggest that the glueballino 



:i:S 



The authors of [438] have pointed out that an extra U(l) boson (the Z') with sufficiently low mass 



could reduce the extreme proton energies required in the Z-burst hypothesis. Speculations about a 



Z' of a few GeV have been entertained as one possible explanation of the NuTeV anomaly |43S , 439 1. 
If the LSP is the gravitino, the gluino decays gravitationally into a gluon and gravitino g — ► 533/2, 
and its lifetime can be long enough. 
35 We note that an alternative proposal (where CRs above the GZK-energy originate from sources at 
cosmological distances) that relies on a supersymmetric extension of the SM, assumes the existence 



of a new exotic pseudo-scalar axion-like particle [ 44G 
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G = gg is the lightest hadron. In the 1980s, G was suggested as the possible primary 
particle in attempts to explain the hadronic nature of CRs from Cygnus X-3 |447| , |448 



There are, however, arguments against a light quasi-stable gluino ||449|| . The main 



problem follows from the fact that along with a neutral hadron G there must exist 
a charged hadron of the type gqqq having baryon number (gluebarino) . The lightest 
gluebarino must be a hadron (guud) in which the triplet of light quarks has total 
spin | and is in a color octet state. The gluebarino must be stable to the same 
extent as the gluino, since the lightest two particle state with the same quantum 
numbers (gg) + (uud) must be heavier than (guud) by the mass of a constituent gluon 
(> 0.6 GeV). Now, production of charged gluebarinos in the Earth atmosphere by 
CRs and their accumulation in the oceans would result in a too high abundance of 
anomalous heavy isotopes of hydrogen and oxygen in contradiction with observational 
data 1 449 1 . However, it may well be that the lightest baryonic state is the neutral flavor 



singlet S° = guds ||383| , [450|| , due to strong quark attraction in this state. We stress 
that even in this case conflict with observational data emerges if gu<is-gluebarino and 
proton are bound into anomalous deuterium. 

In addition, direct searches for glueballino decays |f451| , [452| , |453|| as well as for 
decays of other unstable g-hadrons [|454|1 have severely eroded the attractiveness of 



the light gluino scenario. A practically model-independent analysis, which takes into 
account contributions of the light gluino to the running of a s and to QCD colour 
coefficients excludes light gluinos with m g = 3(5) GeV with 93% (91%) CL [|455|| . By 
combining these analysis with the determination of QCD colour coefficients from the 
analysis of multi-jets events [ 156[| , the conclusion of |[455|| becomes much stronger: light 
gluinos with mass ^ 5 GeV are excluded with at least 99.89% CL. Recent analyses 
based on currently available OPAL and CDF data conclude that the range 3 GeV < 
nig < 130 — 150 GeV can be excluded at the 95% CL ||457|| . For certain choices of the 
parameters, a window in the intermediate mass region 25 GeV ^ m g < 35 GeV remains 



open |[458|, [459|| . This leaves two narrow windows for allowed masses of g-hadrons 



1.5 GeV ^ mg_h ^ 3 GeV and 25 GeV ^ m g -h ^ 35 GeV. Gluino-containing hadrons 
corresponding to the second window would produce EASs very different from those 
detected by current CR-experiments [|384j| . However, light (/-hadrons corresponding 
to the first gluino window produce EASs very similar to those initiated by protons. 
Furthermore, as shown in Fig. ^3| light glueballinos, ttlq <^ 1.5 GeV, have a spectrum 



G 

with GZK-cutoff beyond the currently observed energy range |384 



D. Lorentz symmetry violations with anomalous kinematics 

At present, there is no reason to anticipate the existence of a universal scale below 
which our present notion of flat spacetime geometry is not valid. However, local Lorentz 
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FIG. 23: The three solid curves that accurately fit the observational data from AGASA above 
10 20 eV indicate the diffuse glueballino flux from uniformly distributed sources with injection 
spectra dN$/dE oc E~ 2 ' 7 and intrinsic cutoff E = 10 21 eV for tuq = 1.5, 2 and 3 GeV (from 



below). The proton flux is also shown for comparison. Published in [384] 



invariance (LLI) should not be accepted on faith but rather as a plausible hypothesis 
subject to experimental test. It is possible to introduce the notion of LLI violation ei- 
ther with or without accompanying anomalous kinematics. If no anomalous kinematics 
is involved, any search for LLI non-conservation effects will require testing length scales 
below 10~ 16 cm or less [460]. However, introducing anomalous kinematical constraints 
allows tiny departures from LLI, which would be undetectable at the electroweak scale, 
to be magnified rapidly with rising energy. LLI-violation with anomalous kinematics 



could then help to break the GZK barrier [385, 386, 387 , 388 , 382 . 

For instance, it was noted by Coleman and Glashow ||388|[ that renormalizable and 
gauge invariant perturbations to the SM Lagrangian that are rotationally invariant in a 
preferred frame, but not Lorentz invariant, lead to species-specific maximum attainable 
velocities q for different particles. In such a case, the possible departure from LLI 
can be phrased in terms of the difference between the particle maximum attainable 
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velocities 



5i 



-3 ■ 



Within this framework the energy conservation can be expressed by [|389|| 

4 w > 5 Ap E 



E 



(160) 



(161) 



If Lorentz symmetry is unbroken, Sa p = and Eq. ( |161| ) leads to the conventional 
threshold for a head-on collision. Otherwise, Eq. (|161|) is a quadratic form in E with 



discriminant Aw 2 



ml)5 Ap . For 



5a p > 8(w) = 



Aw 2 



m a — mi 



3.5 x 10 



-25 



W 
W 



(162) 



the discriminant < 0, and therefore the reaction py — > A — >• pre is kinematically 
forbidden for all E. Now, recalling that the relic photons follow a Planck distribution, 
it is easily seen that for 5& P comparable to 5(wq), the GZK effect would be relaxed. 

There are no direct experimental constraints on the parameter 5& p . However, if we 
compare the speed of photons to that of high energy CRs, it is possible to obtain a 
very stringent bound on violations of LLI. This bound follows from the emission of 
Cerenkov radiation and consequent loss of energy by charged CRs, a process which is 
Primary protons (electrons) with energies up to 10 20 eV (1 TeV) 



allowed if c 7 < c CR . 



have been seen, thus 5 P1 < 3.0 x 10" 23 |fH] (<5 7e < 1.3 x 10~ 13 f462j ). Moreover, if 



LLI is broken and c e > c 7 the threshold energy for pair production is altered, yielding 
5 7e < 2m 2 /-E 2 . Multi-TeV 7-ray observations from Markarian 501 then lead to 5 ie < 



1.3 x 10 15 ||462|| . On the other hand, if c 7 > c e , the decay of a photon into an 



electron positron pair is kinematically allowed for photons with energies exceeding 
E ma , x = m e ^2/\5 ej \. The decay would take place rapidly, so that photons with energies 
above E max would not be observable. Detection of primary photons with energies 
£ 7 > 50 TeV leads to 5 ei < 1.3 x 10" 15 " 



462| . The failure to detect velocity oscillations 



of neutrinos translates into \6 VV >\ < 2 x 10 21 , whereas the detection of neutrinos from 
SN 1987a yields |5„ 7 | < 10 -8 ||463| . An additional constraint \8 mi \ < 6 x 10 -22 results 



from a Hughes-Drever type experiment [464]. In the latter, the maximum attainable 
velocity of material matter c m was taken to be the same for all massive particles. It 
is interesting to remark that none of these constraints reaches the level of sensitivity 
to Lorentz violation needed to overcome the GZK paradox. However, the future PAO 
observations of faraway sources could provide constraints on, or even a measurement 
of, the violation of Lorentz symmetry, yielding essential insights into the nature of 
gravity- induced wave dispersion in the vacuum ||465| , |466| , [467] , |468| , \i69[ |470j . 
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VII. SPACETIME'S UNSEEN DIMENSIONS 



In recent years, it has become evident that a promising route towards reconciling the 
apparent mismatch of the fundamental scales of particle physics and gravity is to mod- 
ify the short distance behavior of gravity at scales much larger than the Planck length. 
Such modification can be most simply achieved by introducing extra dimensions, gen- 



erally thought to be curled- up, in the sub-millimeter range |[471] , [472| , [4T3| , |474|| . Within 
this framework the fundamental scale of gravity can be lowered all the way to O (TeV) , 
and the observed Planck scale turns out to be just an effective scale valid for energies 
below the mass of Kaluza-Klein (KK) excitations. 36 Clearly, while the gravitational 
force has not been directly measured far below the millimeter range IfFTSU , SM interac- 



tions have been investigated well below this scale. Therefore, if large extra dimensions 
really exist, one needs some mechanism to prevent SM particles from feeling them. Re- 
markably, there are several possibilities to confine SM fields to a 4 dimensional subspace 
(referred to as a brane- world) within the (4 + n) dimensional spacetime 1 479| , |480| . 



Theories with TeV scale gravity can be classified into two broad categories according 
to whether they do or do not assume that the full spacetime manifold is of factorized 
form. In these theories, in addition to the 4 dimensions that we see (with coordinates 
x^) there are n unseen dimensions with coordinates y m . In the case of non-factorizable 
geometries the metric takes the form 

ds 2 = dx" dx u + g mn {y) dy m dy n , (163) 

where the characteristic size of the large extra dimensions is of the order of the funda- 
mental Planck length, explaining their invisibility. In the case of warped extra dimen- 
sions, the scale of the 4-dimensional metric may vary depending on the location in the 
extra dimension, due to the warp factor e 2A( - y) that can be thought of as the carrier of 
a position-dependent redshift. The metric in this case reads, 

ds 2 = e 2Aiy) dx» dx v + g mn {y) dy m dy n . (164) 

To illustrate the two cases we treat the following representative examples. 

Large extra dimensions: Arkani-Hamed, Dimopoulos, and Dvali (ADD) ||471|| imagine 



the spacetime as a direct product of ordinary four- dimensional spacetime and a (flat) 
spatial n-torus with circumferences Li = 2irri (i = l,...,n), generally of common 
linear size = r c . As mentioned above, SM fields cannot propagate freely in the extra 



36 Theories with small periodic internal dimensions where the low energy degrees of freedom are 
restricted to zero modes and separated by a large mass gap from the massive modes were first 



proposed by Kaluza Jf75| and Klein |47§ |77 . 
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dimensions and remain consistent with observations. This is avoided by trapping the 
fields to a thin shell of thickness 5 ~ M" 1 . Assuming that the higher dimensional 
theory at short distance is a string theory, one expects that the fundamental string 
scale M s and the higher dimensional Planck scale M* are not too different. The only 
particles propagating in the (4+n) dimensional bulk are the (4+n) gravitons. Because 
of the compactification, the extra n components of the graviton momenta are quantized 

h= r ^j 2 - = —-, i = l,...,n. (165) 

Thus, taking into account the degeneracy on ti G Z, the graviton looks like a massive 
KK state with mass 

/ n \ 1/2 

^,..A=(E^) r ?- ( 166 ) 

It is important to stress that the graviton's self-interactions must conserve both ordi- 
nary 4-momenta and KK momentum components, whereas SM fields (that break trans- 
lational invariance) do not have well defined KK momenta in the bulk for £/r c < M s . 37 
Therefore, interactions of gravitons with SM particles do not conserve KK momentum 
components. Applying Gauss' law at r C r c and r r c , it is easily seen that the 
Planck scale of the four dimensional world is related to that of a higher dimensional 
space-time simply by a volume factor, 



M P1 \ 2/n 1 iq /TeV\ /M P1 \ 2/n , x 

2.0 x 1(T 19 — — -rr cm > ( 16T ) 



M* J M* ' \ J \ M, / 

so that M* can range from ~ TeV to M-p\ ~ 10 19 GeV, for r c < 1 mm and n > 2. For 
n < 6, the mass splitting, 

1 / M \ 2/n f M \ n+2/2 

is so small that the sum over the tower of KK states can be replaced by a continuous 
integration. Then the number of modes between \£\ and \£\ + d£ reads, 

dN = dhd£ 2 ...d£ n = S n _i \l\ n - l d£, (169) 

where 

2 7T"/ 2 



37 For ljr c > M s we would expect higher order quantum gravity and string effects to become impor- 
tant, providing some sort of a natural cut-off in the field theory. 
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is the surface of a unit-radius sphere in n dimensions. Now using Eqs. ( |166| ) and ( |167|) , 
Eq. (|169|) can be re-written as 



— ) m n - l dm. (171) 

From the 4-dimensional viewpoint the graviton interaction vertex is suppressed by Mpi. 
Roughly speaking, a m oc Mp 2 . Now, introducing da m /dt, the differential cross section 
for producing a single mode of mass m, one can write down the differential cross section 
for inclusive graviton production 

- 5 7^A 2 i (172) 



dt dm V M * J M * dt 

or else, the branching ratio for emitting any one of the available gravitons 

s n/2 

r »~^ < 173 > 

where s 1 / 2 is the cm. energy available for graviton-KK emission. As one can see by 
inspection of Eq. ( |173| ), the enormous number of accessible KK-states can compensate 



for the Mpj factor in the scattering amplitude. 

Warped extra dimensions: Randall and Sundrum (RS) | }474| | have proposed a simple 
and attractive scenario to ameliorate the hierarchy problem in which our Universe plus 
a hidden brane are embedded in a 5-dimensional bulk with a negative cosmological 
constant. The set-up is in the shape of a gravitational capacitor: two 3-branes with 
equal and opposite tensions rigidly reside at Si/Z 2 orbifold fixed points at the bound- 
aries (y = and y = 7ir c ) of a slab of anti-de Sitter (AdS) space of radius I. The line 
element satisfying this Ansatz in horospherical coordinates reads 

ds 2 = e~ 2lyl/£ Vfiu dx»dx u + dy 2 , (174) 

where rj^ is the metric of Minkowski space. Examination of the action in the 4- 



dimensional effective theory leads to ||474 

M P1 = MH (l - e - 2lTrc/e ) , (175) 

where Mpi is the reduced effective 4-dimensional Planck scale, and M the fundamen- 
tal scale of gravity. The AdS warp factor is exponential in the ^/-coordinate, and so 
the energy scales on the negative tension brane become exponentially redshifted. In 
other words, a field with the fundamental mass parameter m will appear to have 
the physical mass m = e~ nTc ^m on the the 3-brane living at y = 7rr c . There- 
fore, the weak scale is generated from fundamental scales of order Mpi through the 
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exponential hierarchy, requiring r c ji « 12. Additionally, in the "single brane" limit 
(r c — > oo) ||481|| , the bulk continuum can be described by a Conformal Field Theory 
(CFT) |[482| , [483] , |484|| , reproducing the duals discussed in ||485|| . For finite but large r c , 
deviations from conformality are exponentially damped in the infrared ||486|, |487|| . In the 



RS model there exists a KK tower of gravitons with essentially electroweak couplings 
and masses m n = t~ x x n e _7rrc//f , where x n is the n th root of the Bessel function J x 11481 ]. 
This implies that the tower mass spectrum, m n = m\x n /x\, is completely determined 
by the lowest lying excitation. Note that while the zero mode graviton couples with a 
strength M P1 , all the remaining states couple as CMp\ e~ nVc ^) . The phenomenology 
of this scenario (recall that only gravity spills into the extra dimension) is governed by 
two parameters: mi, and c = (£M P i) _1 which is expected to be near, though somewhat 
less than, unity gS8|. The analysis |S8| of Tevatron data QiSiJ g5g, g9TJ, [492], [493], g53 



for anomalous Drell-Yan and dijet production as well as the calculation of indirect con- 
tributions to electroweak observables | [495|| yield m x > 500 GeV, whereas AdS/CFT 
considerations suggest c < 0.1 ||496|| . 

Clearly, the ADD and the RS model would have different manifestations in scatter- 
ing processes. In the case of the product spacetime, each excited state couples with 
gravitational strength, and the key to observing KK-modes in particle collisions is the 
large multiplicity of states, due to their fine splittings. However, in the RS geometry, 
instead of gravitational strength coupling ~ energy/M P1 , each excited state coupling 
is of order energy/TeV, and hence each resonance can be individually detected via its 
decay products. In the first part of this section we discuss how virtual graviton ex- 
change would disturb high energy neutrino interactions ||396||. Under some extremely 



speculative hypotheses this phenomenon may allow neutrinos to interact strongly in 
the atmosphere |[497| , |498|| . More realistic assumptions would lead to a lower cross 
section which would predict an increase in the event rate of quasi-horizontal deeply 
developing showers. In the last part of the section, we discuss the possibility that ul- 
trahigh energy neutrino interactions on the TNB produce gravitons of weak scale mass 
and coupling, resulting in "gravi-burst" fragmentation jets that can contribute to the 



super-GZK spectrum in a way that is similar to the Z-burst mechanism |[499 



A. Influence of KK-modes on the development of extensive air showers 

A novel feature of the contribution of KK graviton exchange to high energy scatter- 
ing cross sections is the fact that it is projectile independent, and thus the same for 
scattering of pp, 7p, up, etc. We focus here on up interactions where the contribution 
from KK-modes leads to a radical departure from SM particle physics. 

A simple Born approximation to the elastic z/-parton cross section ||503|| (which un- 
derlies the total z/-proton cross section) leads, without modification, to a*"' ~ s 2 . 
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Unmodified, this behavior by itself eventually violates unitarity. This may be seen 
either by examining the partial waves of this amplitude, or by noting the high energy 
Regge behavior of an amplitude with exchange of the graviton spin-2 Regge pole: with 
intercept a(0) = 2, the elastic cross section 



da e i \A R (s,t)\ 



2 



dt S 2 

whereas 



„ s 2a(0)-2 „ g2j (176) 



lm[A R (0)] (oj.! 

crtot ~ ~ s w ~ s, (177; 

s 

so that eventually o e \ > a tot- Eikonal unitarization schemes modify these behaviors: in 
the case of the tree amplitudes ||396|| the resulting (unitarized) cross section (n = 2) 



a up ~ s ' the all-order loop resummation yields 0"*°* ~ s 2 ' n ||500|| , whereas for the single 



Regge pole exchange amplitude, ~ In (s/so) ||501 



The other relevant parameter which has direct influence on the shower profile is the 
inelasticity The KK-modes couple to neutral currents, and thus the scattered neutrino 
carries away 90% of the incident energy per interaction. Therefore, the cross section 
has to be at least a few 10 -26 cm 2 to be consistent with observed showers which start 
within the first 50 g/cm 2 of the atmosphere. Specifically, the survival probability iV at 
atmospheric depth X of a particle a with mean free path 

\ a = — , (178) 

®a— air 

is given by 

N(X) = e~ x/Xa , (179) 

where m air ~ 2.43 x 10 -23 g is the mass of an average atom of air, and cr a _ air the cross 
section on air. For a proton energy ~ 10 20 eV, the mean free path is X p ~ 40 g/cm 2 . 
Therefore, a proton air shower is initiated at the top of the atmosphere. The key 
feature in the evolution of the shower is the ratio of decay to interaction of secondary 
hadrons along their path in the atmosphere. The latter strongly depends both on 
particle energy and target density. For protons of this energy, (X max ) ~ 800 g/cm 2 , 
corresponding to 20 hadronic interaction lengths. The neutrino nucleon cross section 
estimates discussed above indicate that X u > 10A P . Due to the larger mean free path, 
X miB is shifted 360 g/cm 2 downwards in the atmosphere. Moreover, in a neutral 
current interaction the neutrino only transfers 10% of its energy to the shower and 



38 Note that the KK contribution to the cross section at s ~ 2 x 10 4 TeV 2 gives ctrk ^ 10 34 
whereas the usual hadronic pp cross section at this energy is a p 1 ^ v 1 n ~ 25 



'VP 



1.5 x 10 J cm 2 , so the impact 



on the pp scattering process is almost innocuous [ 502 1 . 
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consequently elongates the longitudinal development even further ||501| , |504|| . Events 
consistent with these features have not been observed. However, the model may still be 
salvaged if the neutrino-nucleon cross section can be enlarged enough (o vN > 20 mb), 
then multiple scattering within the nucleus may provide a sufficiently large energy 
transfer to reproduce the super-GZK events [|497| , |498| . 



The question of whether the interaction cross section of neutrinos with matter could 
reach typical hadronic values at high energies is yet undecided. The Regge picture of 
graviton exchange is not yet entirely established, the apparently increasing dominance 
assumed by successive Regge cuts due to multiple Regge pole exchange, as well as 
the presence of the zero mass graviton can introduce considerable uncertainty in the 



eventual energy behavior of the cross section [|505| , |506jl . However, it is fair to say that 
recent calculations for the rate of rise of the cross section in the context of string theory 
tend not to support hadronic-like cross sections in ^-interactions [|507| , [508|1 . 



Now we turn to the possibility of experimental input into the issue. In the follow- 
ing we discuss a method to test scenarii in which neutrinos are super-GZK primaries 
which relies directly on EAS observables and is independent of the type of interaction 
enhancing physics. Any physics beyond the SM that increases the neutrino-nucleon 
cross section to typical hadronic values at E v > 10 20 eV should also affect standard 
shower observables at lower energies, where the cross section attains sub-hadronic- 
sizes. In particular, for o vN > 1CT 27 cm 2 , one expects neutrinos to trigger moderately 
penetrating showers, where 1000 g/cm < A max < 2500 g/cm . For the specific case 



discuss in Ref. ||497| , |498|| the cross section is likely to be sub-hadronic near the energy 
at which the cosmogenic neutrino flux peaks (see Fig. and so this kind of shower 
should be copiously produced ||504j| . 39 Clearly, the absence of moderately penetrating 



showers in the CR data sample should be understood as a serious objection to the 
hypothesis of neutrino progenitors of the super-GZK events. Even though neutrino- 
nucleon interactions in TeV-gravity models seems insufficient to explain the showers 
above the GZK-limit, virtual graviton exchanged can still lead to interesting new phe- 
nomena which may be observed in upcoming CR and neutrino experiments [|510| , |511| . 



Recall that the atmosphere provides a detector medium with a column depth ~ 36000 
g/cm 2 for horizontal arrival directions, and so can probe cross sections in the range 
~ 10~ 29 — 10~ 27 cm 2 . Due to the increase column depth, water/ice detectors would 
probe cross sections in the range ~ 10 -31 — 10 -29 cm 2 . 



:!!) 



There is no significant signal of showers with X max > 1000 g/cm 2 in the Fly's Eye data [509 
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FIG. 24: Energy weighted total cross section for hadron production in units of that for the 
Z pole in the Weiler Z-burst model for 7 = as a function of cm. energy. The relatively 
flat lower curve shows the SM, and the upper curve shows the RS model with c = 0.1 and 
mi = 600 GeV. Published in Ref. pSfl . 



B. Gravi-burst 



The RS [474] model of localized gravity would open up new channels for high en- 
ergy neutrinos to annihilate with the TNB to produce a single graviton KK state 
on resonance which subsequently decays hadronically |[499 |. For neutrino masses 
m Uj ~ 1CT 2 — 10 1 eV, and graviton resonance of order a TeV, super-GZK events can be 
produced. Therefore, if one assumes that the incoming neutrino spectrum extends in 
energy with a reasonably slow fall-off, the existence of a series of s-channel KK graviton 
resonances may lighten the requirements on neutrino fluxes given in the Z-burst model 



shown in Fig. 22 



To estimate the event rate of a vv graviton mediated process, consider a point in 
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the parameter space allowed by all current constraints: m\ = 600 GeV and c = 0.1. 
For these parameters, Davoudiasl, Hewett, and Rizzo [f499|| have estimated the cross 



section for vv annihilation into hadrons. Such a cross section has a number of distinct 
contributions, because gravitons not only lead to pairs of quarks, and W and Z final 
states, but also to pairs of gluons and Higgs bosons. By combining all these individual 
process cross sections and assuming the neutrino spectrum above the Z pole energy 
falls like ~ E' 1 , it is straightforward to compute the full energy dependence of the 
total vv — > hadrons cross section, and afterwards the ratio of expected CR rates in 
units of the Z-pole induced rate given in Eq. (|156|) . 



«(^) = F Z = Ml(a m YBl ' (180) 



To get an idea of what this ratio looks like, Fig. |24| shows the case of 7 = 0. Integration 
of R over a range of \fs leads to the relative rate of events expected in the RS model 
to those originating from the Z-burst model. Note that the integral of R under the Z 
pole gives the value unity as it should to reproduce the Z-burst results. If the ultrahigh 
energy neutrino spectrum above the Z-pole falls with spectral index 7 ~ 0.5, it is easily 



seen from Eq. ( 180|) that correction from the first RS excitation to the z/-fluxes required 



to fit the observed spectrum (see Fig. ^2|) is ^ 5%, quickly falling below 1% for 7^1. 



VIII. EXOTICA 

The difficulties so far encountered in modeling the production of ultrahigh energy 
CRs arise from the need to identify a source capable of launching particles to extreme 
energy. In contrast to the "bottom-up" acceleration of charged particles, the "top- 
down" scenario avoids the acceleration problem by assuming that charged and neutral 
primaries simply arise in the quantum mechanical decay of supermassive elementary 
X particles. Sources of these exotic particles could be: 



Topological defects (TDs) left over from early Universe phase transitions associ- 
ated with the spontaneous symmetry breaking that underlies unified models of 
high energy interactions ||512| , |513| , |514| , |515| , |516j , |517|| . 



Some long-lived metastable super-heavy {mx ~ 10 12 GeV) relic particles pro- 
duced through vacuum fluctuations during the inflationary stage of the Uni- 
verse fy||, HI HI H3 • 



40 



For numerical purposes the mass of the Higgs is set to 120 GeV. 
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Due to their topological stability, the TDs (magnetic monopoles, cosmic strings, do- 
main walls, etc.) can survive forever with X particles (mx ~ 10 16 — 10 19 GeV) trapped 
inside them. Nevertheless, from time to time, TDs can be destroyed through collapse, 
annihilation, or other processes, and the energy stored would be released in the form of 
massive quanta that would typically decay into quarks and leptons. Similarly, super- 
heavy relics would also have quarks and leptons as the ultimate decay products. The 
strongly interacting quarks fragment into jets of hadrons resulting in typically 10 4 - 
10 5 mesons and baryons. In this way, very energetic CRs, with energies up to mx, 
can be produced directly without any acceleration mechanism. Another exotic expla- 
nation of the highest energy CRs postulates that relic TDs themselves constitute the 
primaries [|522| , |523f . General features of the exotic scenario have been discussed in sev- 
eral comprehensive reviews [ 170 , |524| , |52l 



overview with emphasis on recent developments. 



528 . Here we give just a general 



A. Top-down origin 



In an epic paper, well ahead of its time, Lemaitre ||529|| introduced what may well 
be the ultimate top-down model. According to Lemaitre the entire material filling the 
Universe, as well as the Universe's expansion, originated in the super-radioactive dis- 
integration of a "Primeval Atom", which progressively decayed into atoms of smaller 
and smaller atomic weight. The CRs were introduced as the energetic particles emitted 
in intermediate stages of the decay-chain. As a matter of fact, Lemaitre regarded the 
CRs as vestigial evidence of the primeval fireworks. Of course, we now know that the 
existence of the CMB precludes the origin of super-GZK CRs in a very early cosmo- 
logical epoch. However, it is amusing that one of the earliest explanations considered 
for the CR origin was a top-down mechanism. 

Topological defects: A more contemporary top down model proposes the vestiges 
of phase transitions in the early Universe as responsible for the highest energy CRs. 
According to current unified models of high energy interactions, the Universe may 
have experienced several spontaneous symmetry breakings, where some scalar field, 
generally referred to as the Higgs field, acquired a non-vanishing expectation value 
in the new vacuum (ground) state. Quanta associated with these fields are typi- 
cally of the order of the symmetry-breaking scale, which in Grand Unified Theories 
(GUTs) can be ~ 10 16 — 10 19 GeV. During a phase transition, non-causal regions may 
evolve towards different states, so that in different domain borders the Higgs field may 
keep a null expectation value. Energy is then stored in a TD whose characteristics 
depend on the topology of the manifold where the Higgs potential reaches its mini- 
mum [|530| , [53 1| , |532| , |533| . The relic defects such as magnetic monopoles [|512| , |517| , 



cosmic strings ||514| , |515|| , superconducting cosmic strings ||513|| , vortons (superconduct- 
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ing string loops stabilized by the angular momentum of the charge carriers) ||534 , |535 
and cosmic necklaces [ |536|| 41 are all relatively topologically stable, but can release part 
of their energy (through radiation, annihilation, or collapse) in the form of X particles 
that typically decay to quarks and leptons. The quarks hadronize producing jets of 
hadrons containing mainly pions together with a 3% admixture of nucleons. Note that, 
in contrast to the case of bottom-up acceleration models, top-down models predict a 
flux above 10 20 eV which is dominated by gamma rays and neutrinos produced via 
pion decay Therefore, the photon/proton ratio can be used as a diagnostic tool in 
determining the CR origin ||538|| . In light of the mounting evidence that ultrahigh 
energy CRs are not gamma rays ffl]| , one may try to force a proton dominance at 
ultrahigh energies by postulating efficient absorption of the dominant ultrahigh energy 
photon flux on the universal and/or galactic radio background. However, the neutrino 
flux accompanying a normalized proton flux is inevitably increased to a level where it 
should be within reach of operating experiments such as AGASA |[539|| . 

Certainly, the precise decay modes of the X particles and the detailed dynamics of 
the first secondary particles depend on the exact nature of the particles under consid- 
eration. However, one expects the bulk flow of outgoing particles to be almost indepen- 
dent of such details. Moreover, the gross features of the hadronic jet systems can be 
reasonably well described by using Local Parton-Hadron Duality (LPHD) [f540| , |541 



In this approach, the primary hadron spectrum is taken to be the same, up to an overall 
normalization constant, as the spectrum of partons in the parton cascade after evolving 
the latter all the way down to a cutoff transverse momentum (&j_)^ ofT ~ ^ ew hundred 
MeV. A common picture for the parton cascade evolution is provided by the so-called 
Modified Leading Logarithmic Approximation (MLLA) of QCD ||542| , |543j| . Within 
this approximation the energy spectrum can be expressed analytically, and simplifies 
considerably when the QCD scale A e g is equal to the transverse momentum cutoff, Qq. 
In such a limiting case the energy distribution of the partons is given by [|540| , |541| , |544j 
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where dN part is the number of partons with a fraction (x, x + dx) of the energy E- ]et of 
the original jet-initiating quark q, £ = ln(l/x), Y = \n(E^ ct / A eS ), and a = [tanh _1 (l — 
2£/y) + i£]. Ib is the modified Bessel function of order B, where B = a/b with 



41 A cosmic necklace is a possible hybrid TD consisting of a closed loop of cosmic string with monopole 
"beads" on it [031 
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a = [lljV c /3 + 2n / /(3iV c 2 )] and b = (11N C - 2n/)/3, n f is the number of flavors of 
quarks, N c = 3 is the number of colors, and Cp = (N% — 1)/2N C = 4/3. Now, using 
LPHD one obtains the fragmentation function due to the hadronization of a quark q, 

x^ = K(Y)x^, (182) 
dx dx 

where x = Eh/ E- ]et . The overall normalization constant K(Y), which takes into account 
the effect of conversion of partons into hadrons, is fixed from the conservation of energy 

x dN h (Y, X ) dx=1 
dx 

Of course, there is a great uncertainty in the extrapolation of the QCD (MLLA + 
LPHD) spectra, which has been tested so far only at collider energies, up to super- 
ultrahigh energies > 10 23 eV. In particular, new processes may alter energy thresholds 
as well as the content of particles in the jets. For instance, if SUSY turns on at an 
energy scale Mgusy of O (TeV), the shower development is expected to tie up not 
only quarks and gluons but also their supersymmetric partners with equal probability, 
provided the 4 momentum transfer Q is above the SUSY scale. Once Q < Msusy, 
the SUSY particles in the cascade would decouple from the cascade, and eventually 
decay into the stable LSPs. In such a case, the final state may contain ultrahigh energy 
LSPs [j545|l , thus modifying the shape of the fragmentation spectrum [|546| . 



It is clear that the wide variety of TDs produce X particles at different rates. How- 



ever, as noted in [516], the X particle production rate, fix, may be parametrized on 



dimensional grounds in a very general way, 



-*W = — (184) 



m x \t / 

where t is the Hubble time, t denotes the present age of the Universe, and Q = 
fix{to) mx is the rate of energy injected in the form of X particles of mass mx per unit 
volume in the present epoch. In most cases p = 1. Exceptions are superconducting 
string models where p < 0, or decaying vortons which have p = 2. Certainly, the 
evolutionary properties of the TD system are unknown, and so one is not able to infer 
the value Qo a priori in a parameter-free manner. However, if a TD scenario is to 
explain the origin of ultrahigh energy CRs, Eq. (|184|) can be normalized to account for 



the super-GZK events without violating any observational flux measurements or limits 



at higher or lower energies, as shown in Fig. [25|. The top-down neutrino and gamma ray 



fluxes depend on the energy released integrated over redshift, and thus on the specific 
TD model. Note that the electromagnetic energy injected into the Universe above the 
pair production threshold on the CMB is recycled into a generic cascade spectrum below 
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FIG. 25: Energy-weighted spectra of nucleons, 7-rays and neutrinos for the TD model with 
mx = 10 16 GeV, p = 1, and decay mode X — ► g + q, assuming an EGMF of 10 -10 G. Thick 
and thin lines represent the SUSY and no-SUSY fragmentation functions, respectively. The 
1 sigma error bar crosses represent combined data from the Haverah Park, Fly's Eye, and 
AGASA experiments above 10 19 eV. Also shown are piecewise power law fits to the observed 
charged CR flux below 10 19 eV, the EGRET measurement of the diffuse 7-ray flux between 
30 MeV and 100 MeV, experimental neutrino flux limits from Frejus and Fly's Eye, as well as 
projected neutrino sensitivities of PAO and NASA's OWL project. Published in Ref. [547]. 



this threshold on a time scale short compared with the Hubble time. Therefore, it can 
have several potential observable effects, such as modified light element abundances 
due to 4 He photodisintegration, or induce spectral distortions of universal gamma-ray 
and neutrino backgrounds ||548 , |549|| . In particular, measurements of the diffuse gamma 
ray background in the 100 MeV region, to which the generic cascade spectrum would 
contribute directly, limit significantly the parameter space in which TDs can generate 
the flux of the highest energy CRs |23|, ggg, ggj, g>g ■ 



Super-heavy relics: The highest energy CRs may also be produced from the decay of 
some metastable superheavy relic particle with mass ^ 10 12 GeV and lifetime exceed- 
ing the age of the Universe ||518| , |519| , |520| , |521]| . Of course, there are no metastable 



superheavy relics within the SM. However, a number of candidate metastable relics are 
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predicted in theories beyond the SM [§5j, [553, EDD1 HI HI S S HI B 
|563| , |564| , |565| , |566|| . Here, the predicted CR flux is driven by the ratio of the density 
of the relics to their lifetime. Since the decay should occur within the GZK-sphere, a 
reasonable parametrization of the decay rate is 



n x 



n x 



;i85) 



where t x is the relic's lifetime, n x = p c Q X h 2 /m x , is the relic density, Q x is the 
cosmic average mass density contributed by the superheavy relics in units of the crit- 
ical density p c ~ 1.05 x 10~ 4 /i 2 GeV cm" 3 , and h is the present value of the Hubble 
constant in units of 100 km sec" 1 Mpc -1 . Clearly, neither Q x nor t x is known with 
any degree of confidence. Additionally, as in the case of TDs, the details of the CR 
spectrum depend on the fragmentation function. Several models have been used to ap- 
proximate the fragmentation in superheavy relic decay, including the HERWIG Monte 
Carlo program | 567 |, MLLA ||568| , |569| , numerical integration of the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equations ||570| , |571 , 572 ], as well as a combination 



of MLLA (for small x) and DGLAP (for large x) ||573|| . In general, these approaches 
all yield similar results in the region 10 19 -10 20 eV ||568|| . As an example, Fig. |26] shows 
the resulting ultrahigh energy CR flux from the decay of relics (m x = 5 x 10 12 GeV) 
clustered in the halo of our Galaxy (more on this below) obtained by numerical so- 
lution of the DGLAP equations, including effects of SUSY ||570|| . The hard injection 
spectra in top-down scenarii (oc E~ l ) have roughly the inverse structure to that of 
the attenuation length of baryonic CRs and 7-rays. Therefore, over the limited en- 
ergy range ~ i0 19,5-20 e y, where the attenuation length drops quickly, the attenuated 
spectra can be easily fitted to the observed spectrum. However, the spectra of such 
models have entirely the wrong shape below the photopion production threshold, and 
at 10 18 ' 7 eV, where the flux is well measured, the predicted flux is more than an order 
of magnitude too small (see Figs. [25] and |26|) . Thus, to be consistent with observations, 
top-down models need to be smoothly matched to a bottom-up model which describes 
the spectrum below the GZK energy ||574|| . An elegant explanation which can account 
for all the events with energy > 10 190 eV can be concocted by mixing the neutrino 

0.07 eV) |75J, |76 



fluxes of top-down models with the Z-burst mechanism 
This hybrid scenario predicts that most primaries above the ankle should be nucleons 
up to 10 20 eV (with a slight accumulation at the GZK energy) and photons at higher 
energies. In addition, the model predicts a new break in the observed spectrum above 
10 20 eV, reflecting the hard top-down spectra. 

Under certain circumstances the superheavy relics could also constitute the dark 
matter (DM) of the Universe. 42 Even though the nature of the DM is still unknown, 



■12 



Nowadays there is mounting evidence that about 90% of the matter in our Universe may be dark. 
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FIG. 26: The best SUSY evolution fit to the CR data with a decaying particle mass of 

. The dotted line indicates the extrapolation of the power law 



5 x 10 12 GeV (c/. Fig. |2_ 
component from lower energies, the dashed line shows the decay spectrum, and the solid line 
is their sum. It is important to stress that this model would have an anomalous accumulation 
factor (see Sec.V). Published in |57Cj]. 



the DM hunt traditionally concentrates on particles with mass of the order of the 
weak scale and with interaction with ordinary matter on the scale of the weak force. 
Importantly, big bang nucleosynthesis constraints imply that such weakly interacting 
massive particles (WIMPs) cannot be baryonic, and must therefore comprise non-SM 
particles. If the WIMP is a thermal relic, then it was once in local thermodynamic 
equilibrium in the early Universe, and its present abundance is determined by its 
self-annihilation cross section. The largest annihilation cross section at early times 
is expected to be oc m^ 2 - This implies that heavy WIMPs would have such a small 
annihilation cross section that their present abundance would be too large. Therefore, 
the mass of a thermal WIMP is found to be less than about 500 TeV [|584]| . Nevertheless, 
it was recently put forward that DM particles may have never experienced local thermal 
equilibrium, and so their masses may be as high as 10 12 -10 19 GeV | [559| , |560| , [56 1| , |562 



Current observations supporting the DM hypothesis include gravitational lensing |577, 5781, peculiar 



velocities of large scale structures [579, 580 1 , CMB anisotropy 581, 582], and recession velocities of 
high redshift supernovae [p83|. 
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563], |564| , [565]. Such monsters have been christened "wimpzillas" |[585]| . One of the 



more promising candidates is the "crypton", the analog of the hadron in the hidden 
sector of supersymmetry breaking string theories ||560|| . Cryptons naturally emerge 
from the theory with about the right mass, they are cosmologically stable, and their 
interaction rate is sufficiently weak such that thermal equilibrium with the primordial 
plasma was never obtained. Furthermore, a sufficient abundance of these particles 
could have been produced near the end of the inflationary epoch. 

If superheavy relics play the role of DM, then irrespective of the abundance fix the 
ratio of the mass density contribution from X particles to that from DM particles should 
be roughly the same everywhere in the Universe, because both X and DM particles 
respond to gravity in the same way. Therefore, since DM particles (by definition) 
cluster on galactic halos, so will the X particles. This suggests that a clean signature of 
the superheavy relic X hypothesis is the anisotropy imposed by the asymmetric position 
of the sun in the Galactic halo ||586| , |587| , |588|| . The exact value of this asymmetry is, 
however, model dependent, because the distribution of dark matter inside the halo is 
by no means certain. The study in ||589|| , following the cusped Navarro- Frenk- White 
(NFW) DM density distribution [f590|| , gives a substantial anisotropy signal in the 
arrival direction of ultrahigh energy CRs. Indeed, the analysis seems to indicate that 
less than 10% of the ultrahigh energy CRs could come from relic particles in the 
Galactic halo. However, recent observations suggest that the NFW profile, for which 
the inner regions are dominated by DM, does not give an accurate description of the 
Galaxy [091 1 . Specifically, the mass density implied by the luminous disk is already 
sufficient to account for the rotation curve in the inner Galaxy without any contribution 
from DM [092 ]. In addition, the microlensing optical depth to the red clump stars 
already shows that almost all the density in the inner Galaxy must be in the form of 
compact objects that yield microlensing, and therefore cannot be DM ||593| , 594|| . 

Recently, Evans, Ferrer and Sarkar (EFS) [|595|1 have analyzed in detail the depen- 
dence of the expected CR anisotropy with the DM density distribution in the Galaxy. 
They studied 4 typical models of the dark matter halo: cusped, isothermal, triaxial, and 
tilted (see Appendix F). In these models, the amplitude of the anisotropy is controlled 
by the extent of the halo, whereas the phase is controlled by its shape. The results of 
EFS, given in Fig. |27|, show that the amplitude, which is ~ 0.5 for a cusped halo, falls 
to 0.3 for an isothermal halo with realistic core radii. The phase points in the direction 
of the Galactic center, with deviation of up to 30° when considering triaxial and tilted 
haloes. These amplitudes and phases are very similar to those obtained by Medina- 
Tanco and Watson in a separate analysis ||597|1 - Unfortunately, the current CR data are 



too sparse above 10 19 eV to have statistically significant discriminators between any 
dark halo model density profiles, or to either confirm or refute any correlation with 
the Galactic halo. Moreover, due to the limited size of the present ultrahigh energy 
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FIG. 27: Contour plots in equatorial coordinates (TZA, S) of the predicted ultrahigh energy 
cosmic ray sky for 4 different dark halo models (from left to right, downwards: cusped, 
isothermal, triaxial and tilted). The effect of the halo of M31 is seen in the upper right of 
each plot: a hotspot at TZA « 00 h 43 m , b ra 41° (referred to the J2000.0 epoch). Published 



in [596]. 



CR sample, nothing can yet be said about the hypothesized existence of an ultrahigh 
energy CR contribution originating in the dark halo of Andromeda (M31) ||595| , 597 1. 
Certainly these issues will be resolved by PAO after a few years of operation. 



B. Monopoles and other topological defects as primaries 



It has long been known that any early Universe phase transition occurring af- 
ter inflation which leaves an unbroken £7(1) symmetry group may produce magnetic 
monopoles ||598| , 599| . For instance, minimal SU(5) breaking may lead to "baryonic 
monopoles" of mass M ~ T c /a, with magnetic charge ?7(1)em and chromomagnetic 
SU(3)c ||600 |. Here a stands for the fine structure constant at symmetry breaking 
temperature T c . These monopoles easily pick up energy from the magnetic fields per- 
meating the Universe and can traverse unscathed through the primeval radiation. Thus, 
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they are likely to generate extensive air showers ||522| , |601| , |602|| . 43 Before proceeding 



further, it is important to point out that if the monopoles are formed at the usual 
GUT scale ~ 10 15 GeV, the energy density overdoses the Universe. Thus, to avoid 
this effect, the symmetry breaking scale associated with the production of monopoles 
has to be shifted to lower energies. Remarkably, if the GUT scale is at ~ 10 9 GeV, one 
would end up with an abundance of relativistic monopoles well below the closure limit, 
and yet potentially able to explain the tail of the CR-spectrum. In addition, for such a 
critical temperature the observed flux of ultrahigh energy CRs is below the flux allowed 
by the Parker limit. 44 Unfortunately, contrary to the observed CR arrival directions, 
the expected flux of relativistic monopoles would be highly anisotropic, roughly aligned 
with the magnetic lines near the Earth f60~5j . 



In models with large extra dimensions, the low-scale unification enables the produc- 
tion of light-mass monopoles, say M ~ 100 TeV. Furthermore, the physical embodi- 
ment of these theories would allow a natural generalization of the 't Hooft-Polyakov 
monopole providing a convenient set of representations for Dl-branes ending on D3- 
branes, and consequently even lighter monopoles [ |606|| . The light-mass monopoles 
could lose and gain energy as they random-walk towards the Earth. The maximum 
energy attainable before hitting the atmosphere is roughly 10 25 eV [ |607|1 . Therefore, 



these "particles" would be ultra-relativistic, and the expected flux would have no im- 
print of correlation with the local magnetic field. Note, however, that direct searches at 
accelerators pretty much exclude masses below a few hundreds of GeV, whereas bounds 



stemming from quantum effects on current observables turn out to be ~ 1 TeV ||608 



The value for the lower limit of the mass of the monopole is still under debate ||609| . 



To mimic a shower initiated by a proton the monopole must transfer nearly all 
of its energy to the atmospheric cascade in a very small distance. The large inertia 
of a massive monopole makes this impossible if the cross-section is typically strong, 
say ~ 100 mb. Wick, Kephart, Weiler and Biermann (WKWB) |607|1 have recently 



pointed out that this problem can be avoided in models in which the baryonic monopole 
consists of g-monopoles confined by strings of chromomagnetic flux. To describe the 
interactions of such a monopole in air, WKWB have developed a model based on the 
four following axioms: i) before hitting the atmosphere the monopole-nucleus cross 
section (unstretched state) is roughly hadronic (cr ~ Aqcd)' attaining a geometric 
growth after the impact; ii) in each interaction almost all of the exchanged energy goes 
into stretching the chromomagnetic strings of the monopole; iii) the chromomagnetic 



43 The idea of monopoles as constituents of primary cosmic radiation is actually quite old, it can be 
traced back at least as far as 1960 [j603| . 

44 This bound requires that there not be so many monopoles around as to effectively "short out" the 



Galactic magnetic held 1 604 ] 
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strings (of tension T ~ Aqq D ) can only be broken to create monopole-antimonopole 
pairs, a process highly suppressed and consequently ignored; iv) the average fraction 
of energy transferred to the shower in each interaction is soft AE/E = r)~ Kqqd/M. 
Generally speaking, in this set-up the monopole will penetrate into the atmosphere, 
because the cross section is comparable to that of a high energy proton. However, 
since the geometrical cross- sect ion grows proportionally with the Lorentz factor T, 
the interaction length after the impact shrinks to a small fraction of the depth of 
the first interaction. Stated mathematically, the unstretched monopole's string length, 
L ~ Aqq D , increases by 5L = AE/T. Recalling that nearly all of the exchanged energy 
goes into stretching the color magnetic strings, the fractional increase in the length is 
SL/L = T, yielding o\ ~ (1 + r)/AQ CD . Now, the total mean free path after the iV-th 
interaction reads, 



\ at rs_/ r^j 

N ^ ^ i-x , v^W 



A 2 A 2 

iv QCD ^QCD 
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where we have assumed a constant density of nucleons (n nuc i (4/3) n A i?jj) and we 
have used the approximation ~ (1 — Aqcb/M) n T ~ T. Here A stands for the mass 
number of an atmospheric nucleus, and Ro ~ 1.2 — 1.5 fm. It should also be stressed 
that for N = r/ _1 the approximation has an error bounded by limjv__ +00 (l— iV^ 1 )^ = e~ l . 
For rj^ 1 1, the total distance traveled between the first interaction and the r/^-th 
interaction is then 



AX 



A 2 v 1 A 2 



yl^^i nr] -\ ( i 8 7) 



r n nuc fi^N Tn 



Note that the mean free path for all secondary interactions is 0(1/T) of the first one. 
Thus, a baryonic monopole encountering the atmosphere will diffuse like a proton, pro- 
ducing a composite heavy-particle-like cascade after the first interaction. A distinctive 
feature of the monopole shower would be the great number of muons among all the 
charged particles [|610|| . Although this feature was observed in a poorly understood 
super-GZK event ||611 |, it seems unlikely that a complete explanation for the ultrahigh 



energy CR data sample would be in terms of magnetic monopoles alone. Moreover, 
any confirmed directional pairing of events would appear difficult to achieve with the 
monopole hypothesis. 

An alternative TD that can be easily accelerated to ultrahigh energies is the vor- 
ton ||523|| . However, the interaction properties of vortons with ordinary matter are 



completely uncertain, and so no clear predictions can be made. 
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IX. THE ATMOSPHERE AS A BLACK HOLE FACTORY 



In this chapter we shift the focus of our discussion from what CRs tell us about astro- 
physical phenomena to what they tell us about fundamental interactions. In particular, 
ultrahigh energy CR neutrinos provide a means to probe the fundamental Planck scale 
in a regime not currently accessible to man-made accelerators. If Planck scale is of 
0(TeV), then cosmic neutrinos could produce microscopic BHs when they interact in 
the atmosphere. The BHs would decay promptly, initiating deeply developing showers 
far above the SM rate [|612|| , and with very distinctive characteristics when the BH 



entropy ^> 1 ||613|| . In the subsequent sections, we discuss first the phenomenology of 
BH production, and then experimental signatures for BHs created in CR events ||253 |. 
For dessert, we examine the parameter space left for p-brane production pTjl |615|. 



A. Black hole production in particle collisions 



Black holes (BHs) are among the most fascinating and inaccessible phenomena in 
nature. It has been known for quite a long time that microscopic BHs can be pro- 
duced in particle collisions with cm. energies above the fundamental scale of grav- 
ity |p!6| , |617| , |618| , |619|| , where they should be well described semi-classically and ther- 
modynamically ||620|| . In the ordinary 4-dimensional scenario, where the fundamental 
scale of gravity is ~ 10 19 GeV, the study of such BHs is hopelessly beyond the realm 
of experimental particle physics. However, if TeV scale gravity is realized in nature, 
production of BHs should be observed in particle collisions with ^Js ^> 1 TeV and 
sufficiently small impact parameter ||621| , |622| , |623| , |624j , |625 . 



The ensuing discussion will be framed in the context of large extra dimensions | 471 |, 
and will be valid at distances that are small compared to the compactification radius 
r c . The requirements for validity of the picture in the warped scenario [[474 



626 



are discussed in Appendix G. For purposes of normalization we use the notation of 
|627|| , where the fundamental mass scale reads 
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The radius of a Schwarzschild BH of mass Mbh = in (4 + n) dimensions is [|628 
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If one envisions a head-on collision involving partons i and j with cm. energy \/I = 
Mbh and impact parameter less than r s , semiclassical reasoning suggests that a BH 
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would be formed [|629|| . 45 The total cross section of the process can be estimated from 
geometrical arguments f|62^ 



625|| and is of order 



a 



(190) 



Criticisms ||63CI| , |631|| of this cross section, which center on the exponential suppression of 
transitions involving a (few-particle) quantum state to a (many-particle) semiclassical 
state, have been addressed in ||632| , |633| , |634| , |635| , |636|| . However, it is worthwhile to 
point out that the heuristic arguments supporting Eq. ( |190| ) only determine a up to an 
overall factor of order one [ 253|| . Uncertainties in this factor are due, for instance, to 
the inclusion of angular momentum when considering BH creation for collisions with 
non-zero impact parameter. 

Note that even though the colliding particles are confined to the brane, if r s <C r c 
then the BH should be treated as a fully (4+n) dimensional object in an asymptotically 
Minkowskian spacetime. In addition, for yfs 3> Mp, BH production should dominate 
over all other SM processes because of the rapidly increasing cross section a oc s 1 /^ 4 " 1 ). 
To calculate the total production cross section, we have to take into account that only a 
fraction of the total cm. energy in a collision is available to the parton participating. In 
our discussion of BH production in the atmosphere, we are most interested in collisions 



of ultrahigh energy neutrinos with the nucleons in the air molecules, where [|612 



dx <Ji(\fxs) fi(x, Q) 



(191) 



Here, s = 2m^E y) fa are parton distribution functions, M^ n is the minimum BH mass 
for which the parton cross section into BHs is applicable, and the sum is carried out 
over all partons in the nucleon. The momentum transfer Q is set to min{MBH, 10 TeV}, 
where the upper limit comes from the CTEQ5M1 distribution functions [ |637j ] . The cross 
section ova^bh is highly insensitive to the details of this choice ||612 |. Once produced, 
the BH will Hawking evaporate with a temperature proportional to the inverse radius 



T, 



n 



n + 1 
47rr s 



(192) 



Note that the wavelength A = 2tt/Th corresponding to this temperature is larger than 
the BH size. Hence, to first approximation the BH behaves like a point-radiator with 
entropy |[253j 

47r M BH r s 

S = . 193 

n + 2 K J 



45 



Here, s = xs, where x is the parton momentum fraction and s is the square of the cm. energy. 
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The BH lifetime estimated from ||638 



dE 
~dt 



A T* +n 



is found to be 



' BH 
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BH 



3+re 
1 + n 



(194) 



(195) 



M D \ M D 

where A is the horizon area. For Mbh ^ M^, the BH is a well defined resonance and 
may be thought of as an intermediate state in the s channel. This means that the 
decay of a BH is only sensitive to the radial coordinate and does not make use of the 
angular modes available in the bulk. Therefore, BHs decay with equal probability to 
a particle on the brane and in the bulk [|622| , |624j , |625|| . 46 Since there are many more 
particles on the brane, we expect the BH to decay visibly into SM particles giving 
rise to events with large multiplicity, (N) ~ M B n/ (2T H ), large total transverse energy, 
and a characteristic ratio of hadronic to leptonic activity of roughly 5:1. As Mbh 
approaches Mp, semiclassical arguments are no longer valid, because the BHs become 
"stringy" and their properties rather complex. 

The magnitude of the entropy determines the validity of the semiclassical approx- 
imation. Thermal fluctuations due to particle emission are small when S ^> 1 ||640|| . 



and statistical fluctuations in the microcanonical ensemble are small for \/S ^> 1 11624 



Searches for BH production at colliders are viable when x m j n = M^/Md is high 
enough that the decay branching ratios predicted by the semiclassical picture of BH 
evaporation are reliable. The QCD background is large, and therefore the extraction of 
signal from background at hadron colliders depends on knowing the BH decay branch- 
ing ratios reliably [§24], §25], §41], §42], §43], §44], §4|, §46], §47], §4|, §4|, §50], §51]. 
Therefore, in searches at colliders a cutoff of x m j n = 5.5 or more may be appropriate. 
In contrast, the search for deeply penetrating quasi- horizontal showers initiated by BH 
decay products can afford to be much less concerned with the details of the final state, 
since the background from hadronic showers is filtered out by the atmosphere. As a 
result, the signal relies only on the existence of visible decay products, which, in this 
context, includes all particles other than neutrinos, muons, and gravitons. Indeed, 
there is very little about the final state, other than its total energy and to some degree 
its multiplicity and electromagnetic component |613| , that we can reasonably expect to 
observe, since detailed reconstruction of the primary BH decay process is not possible 
at cosmic ray detectors. With this in mind, one can choose a significantly lower value of 
than the one needed for collider searches. Although BHs of mass around Mo will 
be outside the semiclassical regime, it seems quite reasonable to expect that they, or 
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One caveat is that BHs may jump off the brane, as suggested in a recent article [ S3S ] . 
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FIG. 28: Cross sections 0Vjv->bh for n = 1, . . . , 7 (bottom to top) for Mr, = 1 TeV, x m - m = 1. 



The SM cross section | 124 | is indicated with a dotted line. Published in Ref. [253] 



their stringy progenitors, will nevertheless decay visibly, whatever stringy or quantum 
gravitational description applies. 

Cross sections for BH production by cosmic neutrinos are given in Fig. ^S[ 47 They 
scale as 



Ml 



(196) 



Despite the fact that the BH production cross section reduces the neutrino interaction 
length to 



X v = 1.7 x 10 7 kmw.e. ( 

a 



(197) 



the Earth's atmospheric depth, which is only 0.36 km w.e. even when traversed hor- 
izontally, is still smaller than X u . Neutrinos therefore produce BHs uniformly at all 



47 It should be stressed that the neutrino-gluon scattering may dominate neutrino-nucleon collisions 
below 10 8 GeV |52 |. 
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atmospheric depths. 



B. Probes of TeV scale gravity in CR experiments 

Given the provocative new features of models with large-compact-dimensions 
(LCDs), there is a strong motivation for phenomenological studies to assess their ex- 
perimental viability. The most obvious consequence of the existence of LCDs is the 
violation of Newton's Gravitational Law at distances of order r c . Applying Eq. (|167 ) 



for n = 1 one obtains r c ~ 10 13 cm, which immediately suggests this case is excluded 
because Newtonian gravity would be modified at the scale of our solar system. For 
n > 2, r c is sufficiently small so that these scenarii are not yet excluded. For n = 2, sub- 



millimeter tests of the gravitational inverse- square law constrain Md > 1.6 TeV |[478 
For n > 3 LCDs become microscopic and therefore elude the search for deviations in 
gravitational measurements. In the presence of LCDs, however, the effects of gravity 
are enhanced at high energies due to the large multiplicity of KK-modes, yielding as- 
trophysical bounds on r c . For instance, the requirement that the neutrino signal of SN 
1987A not be unduly shortened by the emission of KK modes into the bulk leads to 
M d > 50 TeV, for n = 2, and M D > 5 TeV, for n = 3 |5|, |54|. These limits could be 
even more restrictive when taking into account KK graviton decay in typical astrophys- 
ical environments, yielding Md > 600 — 1800 TeV for n = 2 and Md > 10 — 100 TeV 
forn = 3 |>5§ |656], |B5E|| . 

For n > 4, only high energy colliders and CR-experiments provide reasonably sen- 
sitive probes of LCDs. The effects of direct graviton emission, including production 



of single photons or Z's, were sought at LEP ||659| , |660|| . The resulting bounds are 
fairly model- independent, as the relatively low energies at LEP imply a negligible de- 
pendence on any soft-brane damping factor. For n — 4 (6), these null results imply 
Md ~ 870 (610) GeV ||661|| . The effects of low-scale gravity can also be seen through 
virtual graviton effects. These are most stringently bounded by the D0 Collaboration, 
which recently reported [|662j| the first results for virtual graviton effects at a hadron 
collider. The data collected at y/s = 1.8 TeV for dielectron and diphoton production 
at the Tevatron agree well with the SM predictions and provide restrictive limits for 
n > 4. Specifically, if one adopts a Gaussian cutoff on the transverse momentum of 



the graviton due to brane excitations, Md ^ 1 TeV [253] . 



Limits on the fundamental Planck scale derived from CR-observations are mainly due 
to searches for BH production in the atmosphere by ultrahigh energy neutrinos ||253 



363| , |664|1 . The number of BHs to be detected by a given CR-experiment is 



N = I dE u N A — a uN ^ BU (E u ) A{E V ) T , (198) 
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FIG. 29: Bounds on the fundamental Planck scale Md from tests of Newton's law on sub- 
millimeter scales, bounds on supernova cooling and neutron star heating, dielectron and 
diphoton production at the Tevatron, and non-observation of BH production at AGASA. 
Future limits from PAO ground array, assuming 5 years of data and no excess above the 
SM neutrino background, are also shown. The range in Tevatron bounds corresponds to the 
range of a brane-softening parameter which sets the Gaussian cutoff (see text) . The range in 
cosmic ray bounds is for x m \ 



1 - 3. Published in Ref. [252]. 



where A(E U ) is the experiment's acceptance in cm 3 w.e. sr, Na = 6.022 x 10 23 is 
Avogadro's number, d$/dE v is the source flux of neutrinos, and T is the running time 
of the detector. As discussed earlier, for neutrinos with large incident zenith angles 
the likelihood of interaction is maximized. Observed EASs with 9 > 70° typically 
must traverse > 2000 g/cm 2 of atmosphere before interacting. Since the interaction 
length of the atmosphere for protons is ~ 45 g/cm 2 one expects the background from 
hadronic cosmic rays to be eliminated. In 1710.5 days of data collected from De- 
cember 1995 to November 2000, the AGASA Collaboration found 6 candidate events 
with X max > 2500 g/cm 2 | 665 |. At AGASA, the location of the shower maximum 
is determined through correlation to two measurable quantities: r], which parame- 
terizes the lateral distribution of charged particles at ground level, and 5, which pa- 
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rameterizes the curvature of the shower front. The candidate neutrino events must 
satisfy X^ ax ,X^ ax > 2500 g/cm 2 . The expected background from hadronic showers 
is l-72i^o7io!4i) where the first uncertainty is from Monte Carlo statistics, and the 
second is systematic. Of the 6 candidate events, however, 5 have values of X max that 
barely exceed 2500 g/cm 2 , and are well within AX max of this value, where AX max is 
the estimated precision with which X max can be reconstructed. The AGASA Collab- 
oration thus concludes that there is no significant enhancement of deeply penetrating 
shower rates given the detector's resolution. The AGASA results imply lower bounds 
on the scale of low-scale gravity. Of course, the bounds are subject to both small un- 
certainties inherent in the parton level cross section and uncertainty in the cosmogenic 
neutrino flux. For the conservative cosmogenic neutrino fluxes given in Fig. |12| the ex- 
pected rate for deeply penetrating showers at AGASA from SM neutrino interactions 
is about 0.02 events per year, which is negligible. Given 1 event that unambigu- 
ously passes all cuts, and the central value of 1.72 background events, the AGASA 
results imply an upper bound of 3.5 black hole events at 95% CL [|666|| . Therefore, 
the absence of deeply penetrating showers in the AGASA data implies limits on the 
size of LCDs which are summarized in Fig. ^9] | 253|| . These bounds are conserva- 
tive. Note that larger cosmogenic neutrino fluxes, as predicted by some models [B67|, 
will strengthen them, possibly dramatically. Note also that cosmic neutrino interac- 
tions from sub-Planckian extra-dimensional physics can only serve to strengthen these 
bounds. The absence of a deeply-penetrating signal in the Fly's Eye data also implies 
lower bounds on Md- However, these are consistently weaker. The Fly's Eye group 
searched for deeply developing showers in the sample recorded with the Fly's Eye I 



during an observation time of 6 x 10 s [|125|| . Although they found some events with 



large incident zenith angle (6 > 75°) ||668| , the distribution of X 0} the depth of the 



first observed interaction, is consistent with Xq < 2000 g/cm 2 , and there are no events 
with X max > 2500 g/cm 2 ||125|| . This result leads to limits on the size of LCDs. For 



example, if n = 6, Md > 900 GeV [ |663|| . The bounds ||126|| derived with the combined 
exposure of the AGASA and Fly's Eye experiments (the latter integrated over all its 
operating epochs) extend up to 2 TeV for n = 7 and x m { n = 1. Moreover, assuming 
the conservative value x m i n = 3, for which the entropy S > 10, the bounds derived 
with the combined exposure, for n = 5, 6, 7, are Md > 1-26 TeV, 1.30 TeV, 1.40 TeV, 
respectively. All of these exceed the Tevatron bounds ||662 , and represent the best 
existing limits on the scale of TeV-gravity for n > 5 extra spatial dimensions. 

We now discuss the sensitivity of future CR-experiments. PAO is expected to be- 
come fully operational in 2003, and thus would have a running time of roughly 5 years 
before the LHC begins operation. If no enhancement of quasi-horizontal showers is 
seen during the pre-LHC epoch, PAO will probe fundamental Planck scales as large 
as Md = 4 TeV [253 , |612| , |669]| . Moreover, given the prospects for fairly high statis- 
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tics in the region of high entropy, detailed studies of BH shower profiles [|613|| are in 
principle possible. In addition, the projected sensitivity of the OWL satellite project 
will substantially extend the region of Md probed in CR-observations before the first 
collisions at LHC ||670|| . It should also be stressed that neutrinos that traverse the atmo- 
sphere unscathed may produce BHs through interactions in the ice or water. Detailed 
simulations ||671] , |672| , |673| , |674|| also find observable BH rates at neutrino telescopes. 

In summary, in searches for evidence of TeV scale quantum gravity, CR-experiments 
are competitive with colliders, and in fact, until the LHC-era, they will provide the 
most stringent limits. 



C. Prospects for distinguishing BH production in neutrino showers 



Up to now we have only discussed how to set bounds on physics beyond the SM. 
An actual discovery of new physics in cosmic rays is a tall order because of large un- 
certainties associated with the depth of the first interaction in the atmosphere, and 
experimental challenges of reconstructing cosmic air showers from partial informa- 
tion. However, in the following paragraphs we will discuss an observable which skirts 
these uncertainties, and may therefore provide a technique for actually discovering BHs 
in cosmic air showers if TeV-scale quantum gravity exists. If an anomalously large 
quasi-horizontal shower rate is found, it may be ascribed to either an enhancement on 
the incoming neutrino flux, or an enhancement in the neutrino-nucleon cross section. 
However, these two possibilities may be distinguished by separately binning events 
which arrive at very small angles to the horizontal, the so-called "Earth-skimming" 
events ||127| , |128| , |129| , 130 , |131| , |132| , |133|| . An enhanced flux will increase both quasi- 
horizontal and Earth-skimming event rates, whereas a large BH cross section suppresses 
the latter, because the hadronic decay products of BH evaporation do not escape the 
Earth's crust flBl^l . 

Consider a flux of r-neutrinos in the energy decade 10 18 to 10 19 eV. During propa- 
gation inside the Earth, a r-neutrino can produce a tau lepton via a charged current 
interaction. The interaction length for a neutrino near the Earth's surface, which is 
roughly homogeneous with density p s = 2.15 g/cm 3 , is given by 
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is the charged current cross section accurate to within 10% as given by CTEQ4-DIS 
parton distribution [[124|j. Here we neglect neutral current interactions, which at these 
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energies serve only to reduce the neutrino energy by approximately 20%. Therefore, 
for E v ~ 10 19 eV, we have A^ c ~ 250 km. However, if one takes into account the 
possibility of BH production, the mean free path is reduced to 
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The energy of the outgoing r is degraded in Earth by Bremsstrahlung, pair produc- 
tion and deep inelastic scattering, and can be parametrized by 



dE T 

dz 



-(a T + (3 T E T )p s 



(202) 



where a T is negligible at energies of interest, and j3 T ~ 0.8 x 10 6 cm 2 /g [|675|| . The 
maximal path length for a detectable r is, then, 
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where E mSuX « E v is the energy at which the tau is created, and E^m is the minimal 
energy at which a r can be detected. For cosmogenic neutrino fluxes and other rea- 
sonable sources, and the acceptances of typical cosmic ray detectors, taus cannot lose 
much energy and still be detected. For E max / E m i n = 10, L T = 11 km. 

For a given isotropic v T + v T flux, the number of taus that emerge from the Earth 



with enough energy to be detected is proportional to the "effective solid angle" [132 
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(205) 



is the probability for a neutrino with incident nadir angle 9 and azimuthal angle <fi to 
emerge as a detectable r 48 . Here I = 2R^cos6 is the chord length of the intersection 
of the neutrino's trajectory with the Earth, where _R e ~ 6371 km is the Earth's radius. 
Evaluating the integrals, one obtains [ 133 1 
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48 In Eq. ( p05[ ) we have assumed that BHs produced in the Earth will range out before they can 
produce a detectable signal. 
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For E v e (10 18 eV, 10 19 eV), A* ot < R e . Additionally, if the BH cross section is not 
very large, A* ot ^> L T , and thus Eq. (|206|) simplifies to 

\ tot 2 jjr 

«- - 2 'iipF • (207) 

Now, since the r production is proportional to the incoming neutrino flux the 
number of Earth-skimming events detected in 5 years reads 

Nes - C ES - , , (208) 

where Ces is the number of Earth-skimming events expected for the standard cosmo- 
genic flux $g in the absence of BH production. Assuming maximal neutrino mixing 



and the f3 T value given above, one obtains Ces ~ 3.0 for the PAO ground array [131 



The HiRes fluorescence detectors provide additional sensitivity ||132|| . In the following 
discussion we make a conservative assumption of Ces — 3.0 for the SM combined rate 
over 5 years. The presence of a significant cross section for BH production, as well as 
their prompt decay and the rapid absorption of the decay products in the Earth, gives 
rise to a substantial depletion of the ES event rate. 

In contrast to Eq. ( |208| ), the rate for quasi-horizontal showers follows simply from 
Eq. (|198j), and has the form 



QH - G QH-^ Tec ' y IUJ > 



where Cqh = 2.5 for the PAO ground array ||105 



Given a flux and BH cross section <j u n^bk, both N ES and Nqu are determined. 



Contours for these two rates are shown in Fig. |3CJ. As can be seen from the figure, 
it is impossible to differentiate between an enhancement from large BH cross section 
and large flux, given a quasi-horizontal event rate -/Vqh- However, in the region where 
significant event rates are expected, the TVqh and A^es contours are more or less orthog- 
onal, and thus provide complementary information. With measurements of A^qh and 
Nes, both u y jv^bh arid $ y may be determined independently, and neutrino interactions 
beyond the SM may be unambiguously identified ||253|| . 



As an example, consider the case in which <j u n^bk/<?vn = 3, and Q u /&q = 1, shown 
as a dot in Fig. ^0[ On average, one would then observe a total of -/Vqh = 10 deep 
quasi-horizontal showers, an excess of 8 above SM expectations. On average, one also 
expects TVes ~ 0.2 Earth-skimming events. A SM explanation (i.e., ct u n-,byl = 0) of 
the deeply penetrating event rate would require $ v /$q = 4 and predict 12 Earth- 
skimming events, a possibility that would be clearly excluded at high confidence level. 
More generally, one might try to salvage a SM explanation by attributing the observed 
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FIG. 30: Contours of constant number of quasi-horizontal showers Aqh (dashed) and Earth- 
skimming neutrino events N^s (dotted) as functions of source flux Q v and BH production 
cross section ovjv— »bh- 5 year running times for PAO and HiRes are assumed. The figure also 
shows confidence level contours, assuming = 
(A QH , A E s) ~ (10,0.2). Published in Ref. p5|. 



3>q and (Tt/jv— >bh = 3o"^y, corresponding to 



rates to statistical fluctuations in both ]Vqh and N ES . Figure |30] shows contours of 
constant x 2 using a maximum likelihood method for Poisson-distributed data [676]. It 
is easily seen that the possibility of a SM interpretation along the ovtv^bh = axis 
would be excluded at greater than 99.9% CL for any assumed flux. 

BH production will most likely be accompanied by more model-independent sub- 
Planckian effects. In particular, the virtual exchange of bulk gravitons (KK modes) 
leads to extra contribution to the neutrino nucleon cross section. As a consequence, 
the quasi-horizontal event rate will be increased. However, the exchange of spin 2 bulk 
gravitons has very little effect on the Earth-skimming event rate, because neutrinos 
suffer very little energy loss during this process [|500|| . Then one expects KK modes 
to further enhance the ratio Nq H /N ES , making it still easier to distinguish such effects 
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from SM expectations. 

It is important to stress that during this section we only consider the conservative es- 
timates of the cosmogenic neutrino flux. However, some models predict fluxes far above 
these conservative values [[435|| . 49 In such a case the atmosphere would be effectively a 
BH factory, and 100's of BH events would be observed each year at PAO. 



D. p-branes from Heaven 

In this penultimate section we come back to a main topic of this review and discuss 
yet another proposal ||677|| to sneak away from the GZK cutoff. This recent suggestion 



is a variant of those discussed in Sec.VI-A where neutrinos interact strongly beyond 
the GZK-energy. Therefore, the validity of this model is contingent on the observation 
of moderately penetrating showers. 

The basic scaling relation of Eq. ( |167| ) is the simplest one possible. Once one en- 
visages the notion of asymmetric compactification radii, a plethora of new phenomena 
may arise as one can consider the possibility of several compactification scales. Of 
particular interest here is the next simplest example, where there is a hierarchy of two 
sets of compactified extra dimensions, with m small dimensions of length L < 1 TeV -1 , 
and n — m large extra dimensions of length V . The relation between M* and the low 
energy 4-dimensional Planck mass Mpi ~ 1.2 x 10 19 GeV is [ |678|| 

Mp! = M 2 + n L m L' n - m . (210) 

For simple toroidal compactifications, L = 2rcr c , and L' = 2Tcr' c . Scenarii with low 
values of n — m are already tightly constrained from table-top gravity experiments, 
as well as from astrophysical and cosmological considerations ||478| , |653| , |654| , |655| , |656 



16571 BSH1. 



In spacetimes with asymmetric compactifications, extended higher dimensional ob- 
jects wrapped around small extra dimensions, so-called "p-branes", can be produced 



in super-Planckian scattering processes ||614j , |379| , |680|| . In general, the formation of 



higher- dimensional branes dominates the formation of lower-dimensional branes and 
spherically symmetric BHs (0-branes). The decay of p-branes is not well-understood. 
One possibility is that they may decay into lower dimensional brane-antibrane pairs, 
leading to a cascade of branes. In any case, there is no reason to expect them to decay 
only to invisible particles, and it is reasonable to expect their decays, as with BH de- 
cays, to be dominated by visible quanta (those on the brane) that could be detected 



by cosmic ray observatories ||614 . 



49 Clearly, any top-down origin of ultrahigh energy CRs is inconsistent with TeV scale gravity. The 
arguments in this section then do not apply to neutrino fluxes associated with a top-down origin. 
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p-brane production is negligible relative to BH production unless the p-brane wraps 
only Planck-sized dimensions. The enhancement in the p-brane cross section results 
from wrapping on small dimensions and is a consequence of the dependence of the 
p-brane radius ||61 i\ . 



r >=7^ 7(n ' p) Ut) • (211) 

solely on the density of the p-brane. Here, V p = (L/L^y is the volume wrapped by the 
p-brane in fundamental Planck units, 



7(71, p) 



8T 



3 + n — p \ / 1 + p 



(2 + n){2 + n - p) 



p 

(212) 



and L* is the fundamental Planck length. 50 p-brane production dominates over BH 
production only if L < L#. It is worthwhile to point out that in the string-based low 
energy Lagrangian, the gauge coupling is inversely proportional to the compactification 
volume, and a small volume corresponds to strong coupling. In certain explicit models, 
these small volumes can be removed from the gauge sector via a T-duality transforma- 
tion ||S81 |. In what follows, we avoid reference to specific models, and for illustrative 



purposes discuss p-brane production over an inclusive range 0.1 < L/L* < 10. 
The total cross section for brane production is given by 



"VAT— >brane ^ @vN— >p-brane i (^1*^) 
V 

where ^ 

fi/AT^p-brane = V] / dx &i(y/xs) fi(x,Q) , (214) 

is the cross section for p-brane production from neutrino nucleon scattering, 
<3"j :; ^p_brane(V / S) = 7r7-p, is the parton-parton cross section, the sum is over all partons in 
the nucleon, the fi are parton distribution functions, and M™ m is the minimum mass 
required for p-brane production which is set equal to M D . 

It was recently proposed [J677U that ultrahigh energy neutrinos interacting via p- 
brane production may trigger vertical EASs. The required cross sections for hadronic- 
like interactions (~ 100 mb at 10 20 eV), as given by Eq. (|213|) , are only attained for 



n — m = 1, 2 and of a few TeV, a region of the parameter space excluded by the sub- 
mm gravity experiments and astrophysical constraints. Therefore, cosmic neutrinos 



50 



Note that for p — 0, Eq. (211) reduces to the metric of a (4+n)-dimensional BH and r p becomes 



the Schwarzschild radius ]628| . 
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FIG. 31: Contours of fr^ brane (in mb) at E v = 10 11 GeV in the (M D ,L/L*) plane for 
(n = 6,m = 5). The right vertical axis shows the corresponding values M*. The shaded 
region is excluded by the non-observation of deeply penetrating showers at AGASA. The 
hatched region is excluded by requiring no large corrections to standard model physics at 
lower energies (see text). Published in Ref. [615]. 



with interaction strengths enhanced by p-brane production cannot initiate super-GZK 
showers in factorizable spacetimes. Although no explicit models for warped scenarii 
are available, in the spirit of [|626|[ one may consider the possibility that Eq. (pi 3) is 
still valid for M* ~ 1 TeV, and n — m = 1. Under these speculative assumptions, 
hadronic-like neutrino interactions could dominate the scattering to such an extent 
that the non-observation of deeply developing showers is satisfied because the cross 
section is too large — the showers would develop too high in the atmosphere. 



Figure |3l] displays contours in the (Md, L/L*) plane labeled by the value of the cross 
section attained at 10 20 eV (in mb), for n = 6, m = 5 (a similar result can be obtained 
for n = 7, m = 5) ||615|| . The shaded area is excluded by the absence of a significant 
signal of deeply developing showers reported by the AGASA Collaboration [ |665|| . It is 
important to stress that the AGASA data serve to restrict the parameter space even in 
the case where the total cross section at 10 20 eV is of (9(100) mb, because some among 
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the tower of p-branes are produced deep in the atmosphere. We note that the upper 
boundary of this shaded region is in agreement with existing limits ||510|| . In the lower 
left of Fig. [31] is an unshaded region with large cross section which ostensibly can evade 
the AGASA bound because of neutrinos showering high in the atmosphere. However, a 
new consideration enters in this situation: cross sections which are too large can lead, 



via a dispersion relation, to large deviations in SM predictions at lower energies ||682 
With the cross section behavior given by Eq. ( ^13| ), it is a straightforward exercise 



using the result in ||682|| to show that low energy corrections of (9(100%) can arise such 
that 

a„ N (W 20 eV) > 300 mb . (215) 



In Fig. [31] this region is cross-hatched, to indicate that it is problematic for SM physics 
at much lower energies (~ 100 GeV). As can be seen from Fig. 31, this leaves very little 
room for explaining the super-GZK events using p-brane physics: part of the potential 
parameter space is ruled out by AGASA data, most of the rest by requiring small 
feed-down of these high energy contributions to low energy neutrino physics. Only 
a tiny window is available in a region where the small extra dimension L < 0.2L*. 
Furthermore, in this allowed region one has strong coupling effects in the underlying 
stringy regime. 



X. LOOKING FORWARD 

It is difficult - and perhaps hazardous - to speculate where CR physics will go in the 
next ten years. One can never foretell the serendipitous discoveries that will transform 
our understanding. However, on-going and future experiments will surely provide us 
with the statistics to begin discriminating among the many promising ideas so far 
proposed to explain the origin and nature of CRs above the GZK energy limits. The 
superior angular and energy resolution of the Pierre Auger Observatory will allow the 
high end of the energy spectrum and the CR arrival directions to be measured with 
unprecedented precision. These observables provide the most powerful discriminators 
among candidate sources and primary species. 

Distinguishing among CR sources that produce anisotropic distributions is compar- 
atively straightforward. For example, CRs from radio galaxies would appear clustered 
in the super-Galactic plane, whereas those from super heavy relics are expected to 
cluster in the halo. It is more challenging to distinguish among the various models 
that reproduce the observed isotropy below 10 20 eV. In this case, there are two broad 
categories: scenarii where isotropy is realized by diffusive CR propagation, and those 
where the distribution of the sources themselves is isotropic. 

First we comment on the case of an isotropic CR distribution resulting from diffusion. 
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Diffuse propagation of CRs must be associated with relatively nearby sources, and the 
most likely candidates are characterized by distinctive signatures. Galactic sources, 
for example, will emit mainly nuclei which will begin to exhibit a correlation with 
the Galactic plane for energies in excess of 10 20 eV. If Centaurus A is the culprit, 
one can expect to observe neutrons above 10 20 eV which of course would point back 
to the source. In addition, magnetic fields of O(fiG) are required to produce enough 
diffusion to render an isotropic distribution for the lower energy particles. A third 
possibility is that CRs coming from M87 are diffused by a Galactic wind, in which 
case a north-south asymmetry should emerge in the galactic coordinates. Finally, if 
nearby starburst galaxies are the ultrahigh energy CR generators, one should observe 
primarily nuclei which, at high energy, are correlated with the sources. In this case, 
fields of O(nG) would provide enough diffusion to render the distribution isotropic 
below about 10 20 eV. 

An isotropic CR distribution could also arise from isotropic sources, such as highly 
red-shifted AGNs, GRB-fireballs, or decaying TDs. To ascertain whether CRs share 
a common origin with GRBs, we can make use of the data recorded by the HETE 
satellite and directly search for correlations between CR and GRB samples. To fill in 
the picture further, we have to be able to discriminate between the signatures of AGNs 
and TDs. TD scenarii predict that 7-rays should dominate the super-GZK energy 
spectrum. In contrast, AGNs could produce GZK-evading messengers like uhecrons 
or neutrinos that attain hadronic cross sections, both of which would generate typical 
hadronic EASs. These two types of showers can be distinguished by the rate of vertical 
compared to inclined showers. Concerning the uhecron option, we note that the al- 
lowed parameter space is narrowing and collider experiments will have the final word. 
Theories which postulate strong neutrino interactions at super-GZK energies also pre- 
dict that moderately penetrating showers should be produced at lower energies, where 
the neutrino-nucleon cross section reaches a sub-hadronic size. In TeV scale gravity 
models the neutrino-nucleon cross section is likely to be sub-hadronic near the energy 
at which the cosmogenic neutrino flux peaks, and so moderately penetrating showers 
should be copiously produced. Finally, the Z-burst model eludes a definitive confirma- 
tion based on the prediction of the arrival directions, since the distribution will depend 
on whether relic neutrinos are clustered in the halo or not. However, the sensitivity 
of the Pierre Auger Observatory will be sufficiently high to test if the neutrino flux is 
big enough to render the model viable. Interestingly, if the Z-burst model could be 
verified it would provide experimental evidence for the thermal neutrino background. 
Moreover, in this case the CR arrival directions would encode some information on 
the distribution of relic neutrinos in the Universe. Table [I] summarizes the signatures 
corresponding to some of the models in the discussion above. 

Future CR data will not only provide clues to the origin of the super-GZK events, 
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TABLE I: Possible signatures for different ultrahigh energy cosmic ray sources. 





Siprial 


M87 


Asymmetry in flux with respect to Galactic plane 


Pulsars 


Arrival directions correlated with Galactic plane 


Starbursts 


AnmntroTiv towattIk snnrrps ahrvvp 0. x 1 f)^ pV 

iilllDU Ul U U y uWVVcHiJ.HO OUU1 ^VU Oj kJ KJ V Kj _■ /x X \J vT< V 


Cen A 


Spike in source-direction due to n-channel 


Z-burst 


Can only check for required v flux 


Top Down 


Superheavy relics 


Anisotropy towards Galactic halo 




Flux dominated by v's and 7's 


Topological defects 


Isotropy 



but could enhance our understanding of fundamental particle physics. For example, 
if CR primaries are found to have a significant photon component above 10 20 eV, 
this could suggest an exotic ingredient in CRs, such as decay products of TDs, and 
thus could provide insight into the description of the early Universe as well as particle 
physics beyond the SM. On the other hand, an absence of ultrahigh energy photons 
would imply hadronic primaries interacting at cm. energies well above that achievable 
at any current or foreseeable collider, thus providing a unique opportunity to evaluate 
hadronic interaction models at ultrahigh energies. Moreover, a similar technique to that 
employed in discriminating between photons and hadrons can be applied to search 
for signatures of extra-dimensions. In this case, comparison of event rates of quasi- 
horizontal deeply developing showers and Earth-skimming neutrinos will allow better 
limits to be placed on low scale gravity. An optimist might even imagine the discovery of 
microscopic BHs, the telltale signature of the Universe's unseen dimensions. The puzzle 
of ultrahigh energy CRs may even have something to say about issues as fundamental 
as local Lorentz invariance. 

All in all, after 40 years of careful work by many research groups around the world, 
we are in possession of a tantalizing body of data, more than sufficient to stimulate our 
curiosity but not yet sufficient to unravel the mystery of the highest energy CRs. The 
upcoming high quality observations promise to make the next 10 years an extremely 
exciting time for CR physics. 



Note added: After this review was completed, the HiRes Collaboration ||683| , |684 



re-examined their sample and some of the highest energy events where discarded. The 
small modification induced on the energy spectrum is, however, of little help in attempts 
to understand the tail of the CR spectrum, leaving the nature and origin of the highest 
energy events still a mystery. The AGASA Collaboration [|685|1 have also re-evaluated 
their highest energy events and conclude that there are "surely events above 10 20 eV 
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and the energy spectrum extends up to a few times 10 20 eV." Certainly, the conclusion 
of this review remains the same: more data is needed to understand the shape of the 
energy spectrum above the GZK-limit 
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APPENDIX A: HADRONIC INTERACTION MODELS IN THE TERRA 
INCOGNITA 

The analysis of ultrahigh energy cosmic rays requires the extrapolation of hadronic 
interaction models more than 2 orders of magnitude in energy beyond the highest 
accelerator energies {^/s = 1.8 TeV). Actually, the required extrapolation is much 
greater than this because the showers involve nuclei as well as single hadrons both as 
targets and projectiles. 

Soft multi-particle production with small transverse momenta with respect to the 
collision axis is a dominant feature of most events in high energy hadronic collisions. 
Though strict calculations based on ordinary QCD perturbation theory are not feasible, 
some phenomenological approaches successfully describe the main properties of soft 
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diffractive processes (for reviews see Refs. [|686| , |687|| ). 

The most theoretically advanced type of model (like VENUS ||688j| , QGS JET ||689|| and 
dpmjet ||690| , |691|| ) is based on Gribov-Regge theory ||692| , |693|| . The interactions are 
not described by single particle exchange, but by the exchange of highly complicated 
collective modes known as reggeons. The slow growth of the cross section with for 
many different hadronic reactions measured at colliders can be well described by the 
form (jj = Yj s~ v + Xj s e , with the universal parameters i] ~ 0.47 and e ~ 0.08 | |245 



Here, Xj and Yj are relative amplitudes differing for each reaction j and they must 
be determined experimentally. The two terms represent the effect of the exchange of 
reggeons and pomerons, 51 respectively, between the scattering objects, the latter dom- 
inating at high energies. Elastic amplitudes and total cross-sections can be calculated 
on the basis of multi-pomeron exchange, whereas inelastic reactions are introduced by 
cut pomerons ||688[| . The transition from the number of cut pomerons to the secondaries 
produced is modeled by string theory ||694| , |695|| . 
On a different track, in minijet models (e.g. SIBYLL 



696 ) the rise of the cross section 



with energy is driven by a growth of the fraction of hard interactions. The probability 
distribution for obtaining N jet pairs (with p 3 ^ > Pj m , where p™ m is a sharp threshold 
in the transverse momentum below which hard interactions are neglected) in a collision 
at energy a/s, is computed regarding elastic pp or pp scattering as a diffractive shadow 
scattering associated with inelastic processes. The energy of each produced particle 
is generated using Lund techniques ||698|| . The algorithms are tuned to reproduce the 
central and fragmentation region data up to pp collider energies, and with no further 
adjustments they are then extrapolated several orders of magnitude. 

The most difficult point in connection with the air shower development is certainly 
the treatment of nuclear reactions. According to the traditional notion, which treats 



the projectile simply as a superposition of free nucleons ||699|| , fluctuations in a nucleus 
initiated shower should be smaller. This notion still stands if one uses the quark-gluon 
picture of nucleus-nucleus interaction, but the resulting fluctuations in EASs initiated 
by nuclei become nearly twice as large as semi-superposition model predictions ||700 l. 
Such an enhancement is due to the more adequate description of the stochastic behavior 
of projectile nucleons. Moreover, the different approaches used to model the underlying 
physics of pp collisions show clear differences in multiplicity predictions which increase 
with rising energy [701, |702| , 703 , 704]. Therefore, distinguishing between a proton 
and a nucleus shower is extremely difficult at the highest energies |705| . Recently, 
a complete analysis on the hadronic core of showers around the knee done by the 
KASCADE Collaboration ||706| provides important information about the quality of 



51 



This leading Regge trajectory with vacuum quantum numbers was originally proposed by Chew 



and Frautschi [ 397 1 
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the hadronic event generators. The study seems to indicate that QGSJET is the model 
which best reproduces the spectrum below 5 x 10 15 eV. However, above the knee 
deviations in several observables are seen, suggesting that further investigations (and 
data) are needed. 52 



APPENDIX B: REDSHIFTING 

There is now a substantial body of observations that support, both directly and 
indirectly, the idea that our Universe expanded from a super-dense hot phase roughly 
13 billion years ago ||708|| . In particular, in the late 1920's Hubble established that the 
spectra of galaxies at greater distances are systematically shifted to longer wavelengths. 
The change in wavelength of a spectral line is expressed as the redshift of the observed 
feature, 

-, . -^observed /t->-i\ 

l + z = - . (Bl) 

^emitted 

Interpreting the redshift as a Doppler velocity, for z<Cl, Hubble's relationship can be 
re-written as z ~ Hod/c, where Ho is the expansion rate at the present epoch. 

The adiabatic energy losses suffered by cosmic rays due to the expansion of the 
Universe can be associated with the internal energy loss of a relativistic gas with 
particle density n and temperature T doing work to expand its volume V, 

dU=-PdV, (B2) 

where PV = N kT or P = nkT (i.e., nV — N) is the gas pressure and k is the 
Boltzmann constant. Now, since the mean energy of each particle is 3kT, one obtains 
dU = n V dE = —n E dV/3. For a fixed number of particles 

(dE\ (B3) 

WWadiabatie 3 N dt ' 1 ] 

where dV/dt is the expansion rate of the region with field velocity v(r). The change in 
the volume dx dy dz moving with the flux is given by, 

^ = (v x+dx - v x )dy dz + {v y+dy - v y )dx dz + (v z+dz - v z )dx dy, (B4) 
or, using Taylor's expansion 

dV rdV *+^ + ^)dxdydz = (V.v)V. (B5) 



dt \ dx dy dz 



r .2 



For a more extensive discussion on hadronic interaction models the reader is referred to [707 
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Now, the adiabatic fractional energy loss is obtained by replacing Eq. ( p3|) into Eq. ( |B5| ) 

(§) =-\(V.v)E. (B6) 

\ al / adiabatic 6 

In standard Friedman-Robertson- Walker cosmology, with vanishing cosmological con- 
stant and scale factor R which expands at a velocity Vq, v = Vo(r/R), and 



V.i? = 



9 ( ^ 
or 



3|. (B7) 



Hence, 



*f) = ^E = ^E = H E. (B8) 

^/adiabatic R Rdt { ) 



APPENDIX C: KINEMATICS OF PHOTOMESON PRODUCTION 

The inelasticity K n depends not only on the outgoing particles but also on the 
kinematics of the final state. Nevertheless, averaging over final state kinematics leads 
to a good approximation of K n . The cm. system quantities are determined from the 
relativistic invariance of the square of the total 4-momentum of the photon-proton 
system. This invariance leads to the relation. 

8 = ( w c - m - + E c - m f = m 2 p + 2m p w r . (CI) 

The cm. system energies of the particles are uniquely determined by conservation 
of energy and momentum. For reactions mediated by resonances one can assume a 
decay, which in the cm. frame is symmetric in the forward and backward directions 
with respect to the collision axis (given by the incoming particles). For instance, we 
consider single pion production via the reaction p7 — > A — > pir. Here, 

Thus, the mean energy of the outgoing proton is 

{Ep } ~ 2^sm A 2m~ A ' M) 



or in the lab frame 



E (s - ml + m\) (m| - m 2 n + m p ) 
1 ~ ° 2m A 2 m A 



The mean inelasticity = 1 — ((E l ) / E) of a reaction that provides a proton after 
n resonance decays can be obtained by straightforward generalization of Eq. ( |U4jj , and 
is given in Eq. fl53l). 



129 



APPENDIX D: IRON PHOTODISINTEGRATION 



The photo-absorption cross section roughly obeys the electric dipole sum rule 



a{w r )dw r 



60^ MeVmb 
A 



(Dl) 



where n = A — Z is the number of neutrons. 53 These cross sections contain essentially 
two regimes. At w r < 30 MeV there is the domain of the giant dipole resonance where 
disintegration proceeds mainly by the emission of one or two nucleons. A Gaussian 



distribution in this energy range is found to adequately fit the cross section data [255 



At higher energies, the cross section is dominated by multinucleon emission and is 
approximately flat up to w r ~ 150 MeV. Specifically, 



^9(^-2)0(30-^)6-^-^ /iE d e(tw r -30) 

&Ai = 777^ r 



for 2=1, 2, and 



WAi 



&Ai 
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/A6K-30) 
120 



(D2) 



(D3) 



for i > 2 ||255|| . Here, W is a normalization factor given by 



W 



71 



$(^(30 - eoi)/Ai) + HV2(e 0l - 2)/A t 



is the error function, and Q(x) the Heaviside step function. The dependence 
of the width Aj, the peak energy eo?,, the branching ratio fi, and the dimensionless 
integrated cross section £j are given in Ref. ||255|| for isotopes up to 56 Fe. 

The photon background relevant for nucleus disintegration consists essentially of 
photons of the 2.7 K CMB. The background of optical radiation turns out to be of 
no relevance for ultrahigh energy cosmic ray propagation. The infrared background 



radiation | 710 | 



^ = i.ixio- 4 r^r 2 - 5 

dw \eVJ 



cm 3 eV 



(D4) 



only leads to sizeable effects far below 10 20 eV and for time-scales O (10 17 s) ||711[| . 



By substituting Eqs. ( p2| ) and ( p3| ) into Eq. ( pl|) the photodisintegration rates on 
the CMB can be expressed as integrals of two basic forms. The first one is 



h 



A 



2TH 2 h 3 c 2 



15/r 



dw(e 



w/kT 



J + 



i/r 



dw(e 



w/kT 



i-l n-l 



J' 



15/r 



(D5) 



53 Indeed, this integral is experimentally ~ 20 — 30% larger [^0' 
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TABLE II: Series and functions of Eq. Q 



A 




Si 




s 2 


£°° =1 ^J'" 1 exp{-i/rA;r}[$(x/2(2 - e 0i )/A i ) - $(^2(30 - e^/Ai)] 


S 3 


£°° =1 exp[e 2 ] + 15^8/A,) - $(B + v^/A,)} 


S 4 


£°° =1 exp{-15j/rfcT}[(£;T/j)(15/r) 2 + (A;T/j) 2 (15/r) + (fcT/j) 3 ] 


B 


jA. l /TkTV32-2e 0i /V2A i 


K. 


e i Ai $(V2(e i - 2)/Ai) + (A;/2) 2 exp{-2(e 0i - 2) 2 /A?} 



where the functions J and J' are given by the expressions, 



J 



A, 



<S>(V2(2Tw - eoiJ/Ai) + $(V2(e 0i - 2)/A; 



_|_ ( t±L ) J g-2((e 0l -2)/A l ) 2 _ e -2((2r«,-e 0l )/A I ) 5 



and 



J' 



$(^(30 - eoiVAi) + ^(^(eo, - 2)/Ai 



+ 



A, 



e -2((e 0s -2)/A I ) 2 _ e -2((30-e i)/A l ) 2 



The second basic integral is of the form 

!V w 2 dw 



I 2 = (^c 2 Y l a A 



15/r 



e w/kT _ I 



15/r 



r / ./ 15 , r e^-l 



With this in mind, Eq. fluTf ) can be re-written as ||712 
1 



RAi 
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2 h 3 c 2 T 2 2 



A /7T\ 1/2 



e -24/Af (Sl+( 5 2 



(D6) 



(D7) 



(D8) 



^:nrin(i-e- 15 / rA;T ) 



.4 

8 



— r e ~ 2 ^ ( ^) 1/2 + £ OT [ln(l - e~^ VkT ) - ln(l - e" 1 /^) 



A 

2 



+ ^ [ r 2 <S 4 + 15 2 A;Tln(l - e~ 15 / rA;T 



(D9) 



with A, Si, and /C as given in Table I. 
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APPENDIX E: THE NEUTRON CHANNEL 



As remarked in Sec.III-B.2, if the radiation energy density of the source is sufficiently 
high, photopion production leads to copious neutron flux (that can readily escape the 
system) and associated degradation of the proton spectrum. This occurs only near the 
maximum proton energy ||187|| . It is reasonable to assume that the ambient photon 
density of Cen A is sufficiently high so that near the end of the spectrum the efficiency 
of neutron production e n becomes comparable to the proton channel e p . Because of the 
leading particle effect, one expects a cutoff in the neutron spectrum at approximately 
half the maximum injection energy. In what follows, we adopt an energy of 1 x 10 20 
eV as a lower cutoff on the neutron spectrum, and simplify the discussion by assuming 
that in the narrow interval E 2 q G [1, 2] e n ~ e p . The neutron beam observed at Earth is 
further narrowed because of both decay en route and interactions with the CMB. The 
interaction length is approximately 6 Mpc, so that after 3.4 Mpc, 45% of the neutrons 
interact. Half go back into neutrons, so this 22% depletion effect falls within errors. 
However, because of the exponential depletion, about 2% of the neutrons survive the 
trip at 10 20 eV, and about 15% at 2 x 10 20 eV. Note that the increasing survival of 
neutrons at energies above 1.5 x 10 20 eV has, as a consequence of the Cen A model, 
that the observed diffuse flux E 3 J ohs (E) should begin to decrease at these energies 
(unless other factors contribute to an increase). 

We may now estimate a signal-to-noise ratio for the detection of neutron CRs in the 
southern hemisphere. If we assume circular pixel sizes with 2° diameters, the neutron 
events from Cen A will be collected in a pixel representing a solid angle Af2(CenA) ~ 
10 -3 sr. For PAO, the event rate of (diffuse) protons coming from the direction of Cen 
A (say in a 2° angular cone) is found to be [[346 



^ = S Afi(CenA) E 3 J p (E) ^ < ^ events/yr, (El) 
at J El E 6 E 120 



where we have assumed E 3 J p (E) to be (approximately) constant up to at least E 
3 x 10 20 eV, in agreement with the observed isotropic flux in this region, E 3 Job s (E) 
10 24.5±o.2 e v 2 m - 2 s- 1 sr- 1 [f20l. The neutron rate fl34~6 



dNn - JL_ dN p ~d/X(E) _i lfi r f E2 - 2 ° ^20 p -d/\(E) 

dt ~4vrd 2 L dEdt 6 " n 41 L„„ E? 



e -a/x { E) = 116enL41 j ^ e -*W) events/yr, (E2) 



is potentially measurable. For E20 G [1,2] one expects 

—j— ~ 4 e n L 4 i events/yr (E3) 

(JjL 

arriving from the Cen A direction of the sky. With e n L^\ pa 1/2, this gives about 2 
direct events per year, against the negligible background of Eq. (|E1|) ||346|| . An increase 
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in the maximum energy attainable at the source will shift the end of the spectrum to 
higher energy. This will lead to significant enhancement of the neutron flux (at the 
higher energy) because of a greater survival probability. 



APPENDIX F: RECIPES FOR A DARK MATTER HALO 

There is a great uncertainty in the structure of the Galaxy's dark matter halo. 
Strictly speaking, the HI gas rotation of the Galaxy cannot be traced for distances 
<; 20 kpc, and therefore the best estimates of the total mass (~ 2 x 1O 1O M ) and 
the extent of the Galaxy (•> 200 kpc) arise from kinematical analyses of the distant 
satellite galaxies |[713|| . 

Numerical simulations of the structure formation in hierarchical merging cosmogo- 



nies seem to favor a cusped density profile for dark halos | 590 | 



n(r) oc t r-. rrr , (Fl) 

(r + r e ) (r + r s ) 2 v ' 

where r is the position with respect to the Galactic center, r = \r\ is the spherical 
polar radius, r s is the scale radius (which is approximately 10 kpc for the Galaxy), and 
r e ~ 0.5 is set by the resolution limit of the simulations. This profile implies that the 
inner regions are dominated by DM. 

An alternative description, in which the luminous matter dominates the central 
regions and DM the outer parts, assumes the Galaxy's dark halo has an isothermal 
distribution with large core radius (r c ~ 10 kpc) ||714| 



o r 2 _i_ Jl 



Actually, the profile of Eq. ( |F2"|) is the spherical limit of a more general family of triaxial 
distributions [fT15|1 where 



Ax 2 + By 2 + CZ 2 + D 
n[T) K (r 2 + x 2 + y 2 p~ 2 + z 2 q- 2 ^ 2 ' [ ° ' 



with f— (x, y, z), A = (p~ 2 + q~ 2 — 1),B = p~~ 2 (l —p~ 2 + q~ 2 ), C = q~ 2 (l +p~ 2 — q~ 2 ), 
and D = r 2 (l +p~ 2 + q~ 2 )- Here, p and q are axis ratios of the potential, lip = 1, the 
halo is oblate, whereas if p — q, the halo is prolate. In Fig. |27], p = 0.9 and q = 0.75. 
This means that the density contours have semi-axes in the ratio 1 : 0.788 : 0.428, 
yielding a highly flattened profile with an ellipticity 54 of roughly E6. 



54 The ellipticity is 10 x (1 — b/a), where b and a are the projected minor and major axes respectively. 
A spherical galaxy is E0, whilst the most flattened elliptical galaxies are E7. 
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It was also pointed out [|716|| that in some cases the outer parts of the halo may be 
misaligned with the disk, because of the warping of the neutral gas disk. This effect, 
which is known to be present in the Galaxy, leads to a new family of tilted halos with 
density 

n(R> A cc '■ c 2 P + T 2 ) + ^ ,2 + ;V 2 P-T 2 ) (F4) 

where R' 2 = x' 2 + y' 2 . The coordinates (x', y', z') are related to (x, y, z) by a rotation 
through an angle 9 about the x-axis, on which the Sun lies. The pretext for this 
transformation is that the Sun lies nearly on the line of nodes of the warp. In Fig. |27|, 
to obtain the main qualitative features of tilted haloes, the rotation angle was set to 
an extreme value: 6 = 30°. 



APPENDIX G: TUBULAR BLACK HOLES 



In what follows, we discuss the domain of the RS |474]| parameter space for which 
a high energy collision, as viewed from the SM brane, can result in the formation 
of a 5 dimensional flat space BH. Two different kinds of BHs can be produced in 
super-Planckian particle collisions within the RS setup. On the one hand, there is 
the AdS/Schwarzschild BH that can propagate freely into the bulk, and in general 
will fall towards the AdS horizon once produced. On the other hand, there is the 
tubular pancake shape BH that is bound to the brane ||623 |. As in the ADD ||471 



scenario one expects that this type of BH radiates mainly on the brane with interesting 
phenomenological consequences. To characterize the main features of tubular BHs it 
is convenient to re- write Eq. ( |1 74j ) by introducing the coordinate z = £e y ' i , so that the 
metric 

i 2 

ds 2 = — (dz 2 + rjij dx l dx J ) (Gl) 



is manifestly conformally flat ||717[| . Next, to describe the SM brane, define w = z — z c , 
where \w\ G (0, w c ), z c = £ e Tc ^, w c = £ (e Vc ^ — 1), and locate the TeV brane at w — 0. 
A calculation of metric perturbations due to a source on the w = brane can be made 
in the flat 5-dimensional space-time approximation, i.e., ignoring the effects of the AdS 



curvature, if w z c |718| . This implies that the flat space BH formulae discussed in 
Sec.IX-A are also valid in the RS scenario if the Schwarzschild radius as measured by 
an observer in the canonical frame -C z c . The upper limit on the BH horizon translates 
into an upper bound on the BH mass 

< 24 c" 4/3 , (G2) 



M; 



D 



where, M D = (4 vr) 1 / 3 MZ/z cl and c = (£M Pl ) When the energy exceeds this 



bound, the behaviour of the cross section may be analyzed within the AdS/CFT dual 
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picture, and may assume the \n 2 E behavior conforming to the Froissart bound 
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